ae sag ly 

















We Aeud you beat wishes 


for a Happy Holiday Season 
from our Houston Headquarters. 


SCHLUMBERGER 


Well Surveying Corporation 


THE EYES OF THE OIL INDUSTRY 
Sage 





LIBURTOI 





Now you can get OeCeT quality in 
wear-resistant alloy trimmed chokes... 

















In scores of field tests, OeC*T chokes have proved 
superior to any other choke in wear-resistance 


Now you can add greater wear-resistance to the 
other advantages . . . economy, speed, flexibility, 
strength and simplicity . . . that have made O-C-T 
Production Chokes the most popular in the field. 
Designed specifically to get the job done at extreme 
high pressures, these O-C-T Chokes will give you 
trouble-free, long life service under any field con- 
ditions you encounter. They’re available through 
more than 700 supply store outlets . . . because 
more supply stores carry more O-C-T Chokes and 
choke parts than any other make. 


P. O. Box 3091, Houston, Texas 
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CLARK BROS. CO., Olean, New 
York — Gas turbines; engines; 
reciprocating, centrifugal, axial 
flow compressors—gas, steam, 
electric and diesel driven. 
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IDECO, Dallas, Tex. — Hydrair® 
and Power Rigs; Full-View® 
Masts; substructures; single, 
dual and Drive-in Rambler 
Rigs; blocks, swivels, rotaries; 
mud-pumps; petroleum equip- 
ment and supplies. 
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PACIFIC PUMPS, INC., Hunt- 
ington Park, Calif.—Centrifu- 
gal pumps for refineries, 
power stations, pipelines, and 
chemical plants; plunger 
pumps for oilwells. 
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DRESSER MANUFACTURING 
DIVISION, Bradford, Pa.—Pipe 
line couplings, pipe repair 
sleeves and clamps, weld- 
ments and forgings, welding 
fittings, flanges, rings. 


LANE-WELLS CO., Los Angeles, 
California, Houston, Tex., 
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ers and bridging plugs. 
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BLOWER DIVISION, Con- 
nersville, Ind.—Rotary positive 
blowers, gas pumps, centrifugal 
blawers, exhausters, and posi- 
tive displacement meters. 
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DRESSER-IDECO COMPANY 
Columbus, Ohio — Radio and 
television broadcasting towers, 
steel buildings, aircraft hangars, 
mechanical parking garages, 
electric power substations 


MAGNET COVE BARIUM CORP., 
Houston, Tex.- Magcobar®and 
Magcogel*drilling muds and 
other specialized oil well 
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SECURITY ENGINEERING 
DIVISION, Dallas, Tex., 
Whittier, Calif. — Rock bits, 
reamer rock bits and hole open- 
ers; reamers, casing scrapers 
and Neo-Red rubber stabilizers 


In Texas or Trinidad, in faet, wherever the oil, gas and chemical industrie 
are flourishing, you'll find Dresser products and services put to good use. 
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operate individually the better to serve your specific needs. Yet each Dresse1 
company’s brand of dependability is fortified by a great plus factor 
common to all — the Dresser Plus >. It’s a bonus in extra value...an inta 
that’s meaningless unless you’ve experienced it... but added strength, 
knowledge and resourcefulness when you call upon any Dresser company 
for any service. The inter-company teamwork of all nine Dresser 
companies assures the customer of any one the research, design, developm« 
and manufacturing services of all. This is the Dresser Plus. ¥ 

from builders of oil, gas, chemical equipment and services 

that are the standard of comparison the world over. 
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continues to maintain 
its leadership in the 


LUFKIN 


pumping unit field with the | INE 


OF IMPROVED CRANK BALANCED UNITS 





STANDARD 
ASSEMBLIES 





C-640DB-144-30 
C-640DB-120-30 
C-640DB-108-30 
C-456DB-144-30 
C-456DB-120-30 
C-456DB-108-30 
#* C-320D-120-25 
C-320D-84-30 
C-320D-84-27 
C-320D-74-27 
C-320D-74-25 
C-228D-74-27 
C-228D-74-23 
C-228D-74-20 
C-228D-64-23 
C-228D-64-20 
C-160D-74-20 
C-160D-64-23 
C-160D-64-15 
C-160D-54-18 
C-160D-54-17 
C-114DA-64-15 
C-114DA-54-17 


C-114DA-54-16A 
C-114DA-54-15 
C-114DA-54-14 
C-114DA-54-13.5 
C-114DA-48-14 
C-114DA-42-10.5 
C-80DB-48-14 
C-80DB-42-10.5 
C-57D-48-10 
C-57D-42-10.5 
C-40D-40-7.4 
C-40D-34-8 
C-25D-28-6 
C-25D-24-6 

















The new “C” series offers a greater variety 
of combinations of torque, counter-balance, 
and polished rod capacity to suit any 
conceivable well condition. 


LET YOUR LUFKIN REPRESENTATIVE ASSIST 
IN SELECTING THE PROPER AND MOST 
ECONOMICAL UNIT FOR YOUR REQUIREMENTS. 


EXPLANATION OF UNIT DESIGNATION 
C —For Crank Counterbalance 
320D—Gear Reducer (320,000 Lb. In. Peak Torque 


Double Reduction) 
120 —Maximum Stroke (120 Inches) 
25 —Peak Polished Rod Capacity (25,000 Lb.) 


These Units Are Also Available 
With Single Reduction Gear Reducers 


4 


FOUNDRY & MACHINE COMPANY 


LUFKIN, TEXAS 


e: Houston @ Dallas © New York @ Tulsa @ Los Angeles © Seminole © Oklahoma City © Corpus Christi © Odessa 
Kilgore © Wichita Falls © Casper, Wyoming ® Great Bend, Kansas @ Effingham, Illinois © Sterling, Colorado 
Lafayette, Louisiana @ Bakersfield, California © El Dorado, Arkansas 


» equipment in Canada is handled by THE LUFKIN MACHINE CO., LTD., 14321 108th Avenue, EDMONTON, ALBERTA, CANADA 








Quantitative analysis for 
origin of gases in drilling 
AND RESEARCH CORPORATION 
ALAMO NATIONAL BUILDING 
PATENTS PENDING ON 
EQUIPMENT & PROCEDURES SAN ANTONIO, TEXAS 


muds by spectroscopy. 








ATTENTION: Geologists, Petroleum and Electrical Engineers, Physicists, and Electronic Technologists. Employment opportunities exist. Contact us 





HIGH-INJECTION-RATE 


SANDFRAC 


DOWELL treatment injects fluid at rates averaging 30 to 40 barrels per minute 


Here’s a big fracturing treatment designed to pay off big tioners make it possible for Dowell engineers to tailor a 
for you. It’s field-proven Sandfrac® by Dowell, plus the fast, accurate treatment. 

production-boosting force of extra-high injection rates. Dowell research men developed high-injection-rate Sand- 
Powerful Dowell pumping units drive sand-bearing oil into frac to answer problems encountered in many = his 
the pay zone . . . at rates averaging between 30 and 40 high-speed fracturing technique has shown improved 
barrels . . . sometimes as high as 60 barrels per minute. recoveries over long periods in hundreds of wells, 
But even at such high injection rates, the treatment is For service, or for more information, call any of the 165 


precision-controlled. Large-volume, skid-mounted mixing Dowell offices in the United States and Canada; in Venezuela, 


tanks with pressurized suctions, specially constructed bulk contact United Oilwell Service. Or write Dowell Incorpo- 


sand-handling equipment and mechanical sand-oil propor- rated, Tulsa 1, Oklahoma, Department L-15 


services for the oil industry 


A SERVICE SUBSIDIARY OF THE DOW CHEMICAL COMPANY 








This is the fourth of a 
series of editorials to be 
written on “What AIME 
Means to Me’’ by long- 
standing members of the 
Petroleum Branch, AIME. 











EDITORIAL 





WHAT AIME MEANS to ME 


L. E. ELKINS 
MEMBER AIME 


We are often benefactors of an organization cr a sys- 
tem without really appreciating its impact on our lives. 
In retrospect, I see how AIME has helped me profes- 
sionally far more than I had realized. Briefly, AIME, 
both directly and indirectly, has enhanced most of my 
technological environment since my student days. This 
environment produced the mental tools, stimulated their 
use, and always pointed up the challenge that “there is 
a better way.” AIME has also provided the opportunity 
to become acquainted with many leading petroleum 
engineers and geologists throughout the country. 

Upon reflection, it becomes obvious that the least 
a petroleum engineer owes his profession is his mem- 
bership in AIME. So much for generalities — here are 
some specifics. 


AIME—The Educator 


Remove from our engineering schools the literature 
published through AIME, blank out of the educator's 
mind the stimulation he has derived from AIME dis- 
cussions, and education will be seriously handicapped. 
Thus, AIME enters the petroleum engineer's career at 
its foundation. 

I remember very vividly an early experience in my 
engineering growth. A number of young engineers met 
weekly to discuss technical papers, chapters in petro- 
leum engineering books, and articles relative to differ- 
ent aspects of petroleum engineering. Most of these had 
their origin in AIME. During one of these stimulating 
“digging into it” reviews, I received my first clear con- 
ception of fluid displacement mechanics. A once very 
hazy picture came into clear focus. It wasn’t the particu- 
lar discussion, or the one evening that clarified my 
thinking; it was the cumulative effect of many. This is 
just one illustration of the many ways that AIME makes 
its impact on engineers. Local AIME study groups offer 
this same opportunity. This stimulation and cross fer- 
tilization of ideas very likely have had more bearing on 
technical progress than we often recognize. This illus- 
trates perhaps the principal reason for the existence of 
AIME. 


AIME—The Standard Bearer 


I often refer to AIME as the “standard-bearer of 
petroleum engineering technology.” By and large, most 
of the papers which lead to the greatest advance in 
applied technology have first appeared in AIME forums. 
Their selection for a place in the Transactions heralds 
their general good caliber. They therefore become the 
building stones of technology. Quite often a paper may 
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ose some of its direct value several years after its writ 
ing. In most cases, as it loses its direct value, it stimu 
lates a succession of improvements. Thus, the standard 
of technology is steadily advanced — the book is neve: 
closed. Here again, all benefactors of petroleum engi 
neering progress are very much indebted to AIME. 

AIME, locally, regionally, and nationally, has pro 
vided me the opportunity to become acquainted with 
hundreds of petroleum engineers and geologists. Through 
many types of meetings planned by all geographical 
units of AIME, both professional and social contacts are 
provided. AIME would lose its savor if either type o! 
contact were to suffer. Perhaps this professional-socia! 
bond contributes much to the high standard of ethics 
for which AIME members strive. Very few men wil! 
compromise ethics under the watchful eyes of “know- 
ing” friends. Probably, after another 20 years in AIME, 
I will proclaim my AIME friends as the most treasured 
product of my membership. This could not be so were 
it not for the common desire to improve and apply 
technology. 


AIME—The Spokesman 


There is another satisfaction to be gained from AIMI 
Our institute is an organization of mineral technologists 
that can and does speak with influence on subjects 
within its scope. While avoiding politics, AIME does 
attempt to cast true light on issues that bear on the wel- 
fare of its membership and the mineral industry. For 
example, AIME joined with other technical societies in 
advising the War Manpower Commission on how to deal 
with the shortage of engineers. AIME also works with 
other societies in developing an accredited standard for 
engineering schools throughout the country. When rep- 
resentatives of the entire mineral industry speak on 
matters of this nature, they are heard. Sometimes a 
major controversy arises in some field of the mineral 
industry. AIME doesn’t issue a proclamation as to what 
is right. Through the appropriate branch or division, it 
arranges a forum intended to bring out all sides of the 
question and let the members determine what they be 
lieve is correct and fair. 

In appraising what AIME means to me, I have be 
come aware of some deep satisfactions sometimes over 
looked. Truly, AIME is the standard-bearer of tech 
nology and ethics for our profession. AIME stimulates 
an intense drive for high achievement for petroleum 
technologists and industry leaders. 

It is sincerely hoped that others will take time occa- 
sionally to make a similar inventory. wee 
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OIL PRODUCTION — A CONTINUING CHALLENGE 


JOHN R. SUMAN 
MEMBER AIME 


Rock Oil 


In 1855 “rock oil” (petroleum seepages), which sold 
for 30 to 40 dollars a barrel, was used for medicinal 
purposes and widely advertised as “good for man or 
beast.” We are all familiar with Col. Drake’s “rock oil 
test” which was drilled near an “old oil spring” at Titus- 
ville, Pa., in 1859. After drilling into hard rock at 36 
ft, “Drake’s Folly” became the joke of the countryside. 
After many disappointments, the test struck oil at 6912 
ft to open the first oil boom in the United States. Many 
wells were drilled in Pennsylvania, West Virginia, and 
New York during the next few years forcing the price 
of oil down to 10 to 20 cents per barrel. 

Transportation and refining facilities, as well as mar- 
kets (mostly for lighting and lubrication), increased 
steadily during the next 40 years. Consequently, oil 
prices were generally sufficiently low to encourage new 
uses for oil, and they were sufficiently high to provide 
a reasonable profit for the new industry. 

The completion of the gusher at Spindletop in 1901 
by Capi. Anthony F. Lucas ushered in a new era in 
production and marketing, with the use of fuel oil be- 
ing expanded rapidly during the next 10 years. It is 
significant to note that the initial flow of the Lucas well 
was 75,000 to 100,000 B/D—equal to 50 per cent of 
the daily oil production in the United States during 
1900. In 1901 Capt. Lucas presented a paper on the 
Spindletop discovery before the AIME. One of the 
interesting points in his discussion concerned formation 
evaluation. 


Overproduction 


Almost before the excitement of the Spindletop dis- 
covery had subsided, overproduction became a pressing 
problem. During the next 25 years, markets as well as 
daily oil production increased steadily. The ability of 
the United States to supply the fuel needs for the 
machines of war during World War I was considered 
to be an important factor in winning the war. 

Following World War I, oil production continued to 
increase. In 1927, E. B. Reeser, president of Barnsdall 
Oil Co., wrote an article for the Oil and Gas Journal 
entitled, “A Challenge to the Oil Industry.”” The theme 


Paper presented at AIME nation Evaluation Syraposium in 


Houston on Oct. 27-28, 1955. 


CONSULTANT 
HOUSTON, TEX. 


was: “Overproduction must be regulated by cooperation 
or government supervision. Leaders should develop a 
plan.” 

During the great depression of 1929-33, overproduc- 
tion grew worse and prices dropped to ruinous levels. 
By Dec., 1931, the situation had reached a critical 
stage; enormous flush production of the East Texas and 
Oklahoma City fields could not be absorbed by an 
already oversaturated market. Governors of Texas and 
Oklahoma attempted to shut down and control produc- 
tion by military force and martial law. This chaotic 
situation was gradually ironed out somewhat by inter- 
ested industry groups. In 1935, the Interstate Oil Com- 
pact Commission was formed for the purpose of deter- 
mining means that should be taken to conserve oil and 
gas, and for promoting the enactment of legislation and 
enforcement by competing states 


World War Il 


Excess potential continued to be a problem until 
1941. when the demand for oil skyrocketed during 
World War II. In addition to the oil and regular gaso- 
line furnished for the vehicles of war, the industry met 
the challenge to furnish the tremendous demand for 
100 octane gasoline and synthetic rubber. Supplying of 
fuel oil to power the huge United States Navy was in 
itself a real accomplishment 

At the beginning of the war, the OPA froze the price 
of oil at prewar levels. This freeze caught oil industry 
prices at an abnormally low point and, coupled with 
the shortage of steel, resulted in a relatively small oil 
field development program during the war. Although 
the industry emerged from the war with its producing 
reserves low, its plant facilities scrambled, and its mar- 
kets disrupted, praise for oil industry contributions to 
the war effort was universal 

The transition from war to peace was made in a 
relatively short time without great difficulty. The post 
World War II period in the United States has been 
marked by greatly increased cost of finding oil and very 
large drilling programs which resulted in the finding of 
large reserves. 


Technology Moves to Meet the Challenge 


Probably the most serious and continuous problem 
which plagued the industry during the first quarter of 
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this century was the political hue and cry raised against 
“waste of natural resources.” This situation created a 
particular impetus for scientific study of producing 
problems, including technologic improvement in reser- 
voir control. It was one of the factors which motivated 
the formation of a number of technical organizations 
to study industry problems. 


Technical Organizations 

1. The Petroleum Div. of the Bureau of Mines, estab- 
lished in 1912, was the first organization in the history 
of the petroleum industry to make a concerted co- 
ordinated effort to understand the mechanics of oil 
production. From its inception it emphasized the im- 
portance of fundamental understanding of the physical 
processes involved in oil production and of the applica- 
tion of this knowledge to the development and exploita- 
tion of reservoirs. Near the end of World War I, field 
and laboratory studies were initiated by Bureau of 
Mines’ engineers and scientists which were the precur- 
sors of our modern industrial production research and 
reservoir evaluation programs. Unquestionably the Bu- 
reau of Mines held the dominant position in the scien- 
tific study of reservoir behavior and control prior to 
1927. 

2. In 1914 the AIME set up a “Committee on Petro- 
leum” with Capt. Anthony Lucas as chairman and 
began to hold annual forums for discussions on tech- 
nical problems of the oil industry. 

3. In Feb., 1917, the Southwestern Association of 
Petroleum Geologists (changed to American Associa- 
tion of Petroleum Geologists in 1918) was formed for 
the promotion of the science of petroleum and gas. 


4. In 1919 the American Petroleum Institute was 
organized and soon became an active force in the pro- 
motion of technological developments within the in- 
dustry. 


Federal Oil Conservation Board 

When the Federal Oil Conservation Board was ap- 
pointed by President Coolidge in Dec., 1924, the API 
was Officially invited to assist the Conservation Board 
in assembling data and reaching conclusions in regard 
to oil conservation measures. An initial API committee 
of 11 members was appointed to prepare the industry’s 
case; however, final recommendations of this committee 
were somewhat negative because of their suspicion of 
the motives of the Conservation Board. 

As the work of the Conservation Board continued 
along constructive lines, the API's attitude changed and 
a new committee with E. W. Marland, president of the 
Marland Oil Co., as chairman was appointed to work 
with the Conservation Board. Members of this com- 
mittee were Henry McGraw, William S. Farish, Thomas 
A. O'Donnell, Sid H. Keoughan, J. Edgar Pew, and 
L. P. St. Clair. The Marland Committee was composed 
largely of non-technical men, and it was soon found 
that a technical working group was needed to solve 
some of the more complex and difficult technical prob- 
lems facing the industry. Accordingly, regional com- 
mittees were organized with M. E. Lombardi, chairman 
for the Pacific Coast; Max W. Ball, chairman for the 
Rocky Mountain area; John R. Suman, chairman for 
Texas and Louisiana; W. P. Haseman, chairman for 
Kansas and Oklahoma; and Earl Oliver, general secre- 
tary. The Marland Committee reported its findings to 
the API in Dec., 1927. Following are some of the reso- 
lutions proposed by the Marland Committee: 


“Resolved that . . . the industry adopt every reason- 
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able and practical means to effect the conservation o! 
natural gas in oil bearing formations in the interest o! 
economies in production costs, increased ultimate re 
covery, and in the interest of conserving our supplies 
of petroleum, and to the end that proper recognition 
of the correlative rights of land owners and lease owne: 
in a given pool might be attained. 

“Resolved that the American Petroleum Institute en 
courage the cooperative development of oil pools whe: 
ever found under diversified ownership, as the bes 
means of effecting the greatest economies in production 
of oil, and the best means of conserving our petroleum 
resources; and that to this end the Institute should en 
courage passage of the necessary state legislation to 
permit cooperative development of oil pools .. . 

“Resolved that the Institute should impress upon its 
members in the several states, where oil is being pro 
duced, the need for state legislation or other state regu 
latory measures to permit and require the producers o! 
oil and gas . . . to adopt measures to prevent the waste 
of natural gas found in oil sands.” 

The findings of the Marland Committee were as 
sembled by H. C. Miller of the Bureau of Mines and 
published in 1929 in a book entitled Function of Nat 
ural Gas in the Production of Oil, under joint sponsor 
ship of the U. S. Bureau of Mines and the API. 


Some Significant Industry Developments 
1901 - 1955 


From 1901 until 1951 many significant developments 
took place in the oil industry. Some of them are listed 
below. 

In 1903 Union of California used cement in Cali 
fornia oil wells to shut off water. 

By 1920, a technical development of formation evalu 
ition included use of strip formation logs for correla 
tion, and peg models. Coring was becoming widespread 
One of the earliest types of coring equipment used in 
the Southwest was the “poor boy” or “Texas” core bar 
rel. This consisted of a short piece of pipe with 
V-shaped teeth cut on its lower end and run on bottom 
of drill pipe. After a core was cut, excess weight was 
set down on the pipe to bend the V-shaped teeth in 
ward, thus holding the core in place. 

By 1921 core analyses and drill-stem tests were in 
general use, and Elliot's core barrel was being em 
ployed. This core barrel, having a stationary inner bar 
rel which greatly increased core recovery, is considered 
to be one of the most important developments of th« 
period. 


Search for Oil 

Although studies of subsurface strata of the earth 
had been made by physical measurements on the sur- 
face much earlier, it was not until after the close of 
World War I that suitable apparatus was developed to 
carry on such studies in the search for oil. 

1. In 1921 or 1922 a group of Germans (L. Min 
trop, et al), having developed portable refraction seis- 
mic equipment, introduced seismograph crews into 
Mexico and the United States. During this same period 
Anglo-Persian and Dutch Shell groups were making 
trial surveys with the torsion balance. In Dec., 1922. 
torsion balance instruments were first introduced into 
the Gulf Coast region, the first work being carried out 
(according to DeGolyer) at Spindletop. 

2. Although earlier work had been done with the 
refraction seismograph in Mexico and the Mid-Con 





tinent area of the U.S., it was not until March, 1924, 
that work of this type was done in the Gulf Coast. The 
first crew was obtained from Seismos Co. (Mintrop) 
and was under the general direction of Alexander Deus- 
sen for Marland Oil Co. 

3. During the period 1924-1930 an intensive and 
highly successful campaign of searching for shallow 
salt domes, varying from a few hundred to 5,000 ft in 
depth, was carried on by use of refraction seismograph 
and torsion balance instruments. 

About the same time the seismic and torsion bal- 
ance systems were being introduced, the magnetometer 
(also a German development) was put into use. The 
results obtained by the magnetic method were not as 
revealing and not nearly as spectacular as those shown 
by the gravity and seismic methods. The magnetometer, 
however, proved to be a useful instrument for explora- 
tion in regional type studies, particularly to indicate spe- 
cial types of structure such as serpentine plugs. 

5. J. C. Karcher began work on the reflection seismic 
method as early as 1919, and some field work was done 
in Oklahoma in 1921. However, it was not until 1928 
that the method was sufficiently perfected to be put into 
commercial use in Oklahoma. The reflection seismic 
method spread rapidly, and by 1931 it had replaced the 
refraction method almost completely. A very active 
campaign of reflection work began in the Gulf Coast 
region in 1931. 

6. As a culmination of several years of research ef- 
fort, principally by O. H. Truman, the first successful 
gravity meter was built by Humble in early 1931. After 
a period of development to make the instrument more 
field worthy, the lighter and faster gravity meter re- 
placed the torsion balance in the conduct of gravity 
surveys. 


Other Significant Developments 

In 1929, Marcel and Conrad Schlumberger brought 
their new electric logging tool over from France. 

In the lean early 30’s there was an accelerated trend 
toward greater portability of rotary rig equipment. 

One of the earliest uses of directional drilling in the 
Southwest was Humble’s killing of a wild well in Conroe 
in Jan., 1934, by drilling within 30 ft of the bottom of 
the wild well and pumping in water to kill it. To the 
best of my knowledge, this was the first successful 
effort to ki! a cratered oil well by this procedure. 

Development in radioactivity logging began in the 
late 30’s and was in commercial use by 1940. 


Reservoir Engineering 

The growth of reservoir engineering during the past 
25 years has probably been the most important indus- 
try development during the past 100 years. This new 
science has caused a fundamental shift to formation or 
reservoir evaluation instead of thinking about individual 
well behavior. Some of the pioneers in forming cur- 
rently accepted basic reservoir concepts were: Stewart 
Coleman, H. D. Wilde, Jr., T. V. Moore, S. E. Buck- 
ley, M. C. Leverett, Morris Muskat, Ralph Schilthuis, 
A. F. Van Everdingen, H. G. Botset, and R. D. 
Wyckoff. 

For a number of years a fundamental mathematical 
challenge to reservoir engineers has been the derivation 
of equations that describe the flow of reservoir fluid 
systems and the development of practical solutions for 
such equations. Basic equations were developed to pre- 
dict the rate of production and ultimate recovery of 
each of the fluids from a given reservoir under natural 
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or ortifically imposed conditions. Undoubtedly the de- 
velopment of electric analog computing machines, in- 
cluding the oil pool electric analyzer by Carter Oil Co. 
in 1942, was a major reservoir engineering achieve- 
ment. Within recent years the use of digital computing 
machines has made it possible to perform mathematical 
operations at rates hundreds of times those possible with 
a desk calculator. These machines have enabled the res- 
ervoir engineer to effect more complete (and, therefore 
more reliable) solutions of problems previously con- 
sidered insoluble because of the magnitude of the com- 
putation work involved 


Brilliant Scientific Achievement 
Although the science of reservoir engineering is not 
complete today, the advances in knowledge of reservoir 
behavior and control which emanated from the many 
competing companies of the industry between 1927 and 
World War II constitute a record of brilliant scientific 
achievement. The increased recovery of oil which has 
resulted is an excellent example of the progress that 
can be made by private enterprise operating in a climate 
cf economic freedom 
It is worthy of note that this era included one of 
the most serious economic depressions in the history 
of the United States, a condition that might have teen 
expected to discourage industrial research. 


Shortage of Engineers 

Unfortunately, one of the principal problems facing 
the industry today is the shortage of engineers and 
experienced teachers to train additional engineers. Cur- 
rent registration of engineers in the various colleges 
indicates that the number of engineering graduates will 
increase considerably during the next few years. In 
the meantime, maximum benefits from available talent 
can be realized by assigning engineers to strictly engi- 
neering duties, rather than requiring them to perform 
semiclerical and other non-technical tasks. 

There does not appear to be an immediate solution 
to the instructor shortage, however. Fewer technical 
graduates are entering the teaching profession because 
positions with industry are more attractive financially. 
Some companies are attempting to alleviate this situa- 
tion with direct grants to the universities and colleges, 
retaining faculty members as consultants, and making 
part-time employment available. Still, it appears that 
the present help for colleges by industry is insufficient. 
New ways should be found to insure that adequate 
competent instructors are available to properly train 
engineers for the industry 


Conclusions 


It is obvious that the great advance made by the oil 
industry could not have been possible without the won- 
derful cooperative spirit existing over the years between 
the oil producers, the service companies, the govern- 
mental agencies such as the Bureau of Mines and The 
Texas Railroad Commission, and the many industry 
technical groups including the AIME, API, and AAPG, 
plus the continued active support by the many excel- 
lent universities and colleges. 

This Formation Evaluation Symposium is a fine ex- 
ample of such cooperation, the system of free enterprise 
in action, with the entire Symposium being sponsored 
by the University of Houston and the AIME. With 
this type of cooperation, the industry cannot fail to 
continue to move forward—under forced feed of rapid 
technical advances. tok 
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FIRST PACIFIC COAST PERMANENT DRILLING ISLAND 
MODIFIED to ALLOW DRILLING of 70 WELLS 


E. E. PYLES 
ASSOCIATE MEMBER AIME 


lhe first permanent drilling island constructed on the 
Pacific Coast is now being modified to provide surface 
locations for the drilling of a total of 70 wells. The 
small circular Monterey Oil Co.-The Texas Co. filled 
island, 142 miles offshore from Seal Beach, was origin- 
ally designed for the drilling of only three wells. Now, 
through the construction of a raised floor to provide 
additional cellars, this small island is expected to pro- 
vide adequate surface drilling sites from which to de- 
velop the lease. 

When this modification has been completed, about the 
end of this year, the little island will provide probably 
the densest concentration of drilling sites anywhere in 
the world. 

The circular island is 75 ft in diameter and is located 
in 45 ft of water. A rectangular wharf 74 ft by 51 ft 
extends shoreward from the island. A group of four 
pipelines connects the island to a separation and ship- 
ping center onshore. Nine wells have been drilled from 
the island to date, and the tenth one is now being drilled. 

The outer rim of the island is formed of interlocking 
sheet steel piling driven into the ocean floor to a depth 
of 15 to 20 ft. The interior is filled with rock and sand. 
A series of parallel reinforced concrete beams across 
the island carries the working load. Each beam is sup- 
ported by a line of wood piling. The reinforced concrete 
floor rests on the beams and sand fill. The construction 
of the raised floor above the original concrete floor is 
now in progress. Through this addition the island, which 
was Once considered supplemental only, will become a 
very valuable surface property from which to develop 
the extensive oil reservoir. The operators have spent 
approximately $750,000 in connection with island con- 
struction. 


Early History of the Lease 


In 1944 the Marine Exploration Co., one of the pre- 
decessors in interest of Monterey Oil Co. and The 
Texas Co., obtained rights to drill! for oil offshore at 
Seal Beach by obtaining a state lease which was the 
predecessor of the present State Lease PRC 186.1. This 
lease begins at the Seal Beach shoreline and extends 3 
miles seaward. The drilling of the first two wells was 
done from a land location in the City of Long Beach, 
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MONTEREY OIL CO. 
LOS ANGELES, CALIF. 


County of Los Angeles, 2,500 ft inland from the wate: 
edge. 

The first well drilled from the land site was unsucces 
ful. The second well was also non-productive, but was 
redrilled to a location 9,271 ft seaward from the surface 
location and bottomed in the state lease, which is in 
Orange County. The deviation of this well reached 8 
degrees from vertical and the drilled depth was 12,187 
ft. The perforated interval was from 11,447 ft to 11,89! 
ft, but because of the very high deviation from vertic 
the pump was run to only 7,142 ft. The well was com 
pleted in Nov., 1948, pumping 41 BOPD. 

In a number of respects this well was unique. It pri 
duced oil from a greater distance seaward than an 
California well so far drilled. At present, after five ye 
production, it still yields 25 BOPD of 26° API gravity 
oil. 

While the second well was being drilled from the loca 
tion on shore, the California State Lands Commissio: 
by resolution instructed the company that all futur 


Fig. | — Members of the Southern California Petroleum 
Section of AIME view The Monterey Oil Co.-The Texas 
Co. filled island during a recent field trip. 





wells drilled into this state lease be drilled from a 
filled island offshore. 


Building of the Island 


The results of the second well justified further ex- 
ploration of the leased area and, therefore, the opera- 
tors decided to build a filled island on the lease in ac- 
cordance with the State Land Commission’s orders. 
Legal entanglements prevented the issuance of a per- 
mit by the United States Engineers until 1952. Finally, 
the permit was granted and the building of the island 
was started. 

The location of the island was chosen so that when 
enlarged later it would fall in line with the proposed 
extension of the Long Beach breakwater. Also, the 
island would present no obstruction to shipping and 
would be protected from heavy seas by the breakwater 
at the selected location. 

The island was designed by Capt. George F. Nichol- 
son of Long Beach to provide surface locations for the 
drilling of three wells. Sand was built up into a mound 
by a dredge, and then wood piles and the conductors 
for three wells were driven. At this point the work was 
stopped by injunction resulting from a claim by the 
City of Seal Beach to jurisdiction of offshore land to 
the 3-mile limit, and an ordinance preventing drilling in 
the City of Seal Beach. In a final decision, the Califor- 
nia appellate courts eventually ruled that the City of 
Seal Beach had no jurisdiction over the tide and sub- 
merged lands owned by the state. This litigation caused 
a delay of 14 months, during which time the partially 
completed island was maintained at a cost of nearl 
$100,000. 

The driving of the wood piles was completed and 
outside spuds were driven to support the circular beams, 
or waling, around the island. This waling was used as 
a guide for the driving of interlocking sheet steel piling, 
which was driven to refusal, and which formed the 
circular shell of the island. 

The space between the piling and three conductor 
pipes and the circular shell was then filled with 6,000 
tons of small Catalina rock and 6,500 tons of sand. 
Gunited wood piling was then driven for the 74 ft by 
51 ft rectangular wharf. Protection on the outside of 
the shell was provided by placing 14,600 tons of rip-rap 
around the shell to within 10 ft of water level. The rip- 
rap slopes outward to a distance of 45 ft from the island. 
Rip-rap was also placed under the wharf, but only to 
within 26 ft of the water level. 

Seven reinforced concrete beams were poured, sup- 
ported by the rows of wood piling. The concrete floor 
for the island and wharf was poured simultaneously. 
The floor thickness varies from 6 in. to 9 in. Space was 
provided for a cellar 6 ft deep by 8 ft wide and 11 ft 
long, within which were the conductor pipes for the 
three wells; space was also allowed for a 12 * 12 X 12 
ft mud pit and sump. The concrete floor was built with 
a slight slope so that liquids on the floor would drain 
toward the sump. 

A concrete wall extending 8 ft above the floor level! 
was built around the island except on the wharf side 
where it was reduced to a curbing. This wall was to 
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Fig. 2— Plan and cross-sectional drawings of The 
Monterey Oil Co.-The Texas Co. filled island. 


protect the crew from wind and spray during rough 
weather. 

The island was so designed that all the major loads 
to be carried would be borne by reinforced concrete 
beams, which were, in turn, supported by rows of wood 
piling. The type of design selected included such large 
factors of safety that there was no special need for 
concern about earthquake and storm loads on the struc- 
ture as a whole. A sandy ocean bottom and convenient 
source of rock at Catalina Island made the design eco- 
nomical as well as sturdy. The island is capable of last- 
ing 50 or more years, as may be required to produce 
the lease to an economic limit 


Equipment and Facilities on the Island 


The drilling rig was specially designed to operate in 
the restricted space available on the island. Before mov- 
ing any drilling equipment to the island, the drilling 
contractor laid out an area the exact dimensions of the 
island working area, and assembled and unitized the 
rig within this area. Reinforced mud tanks were set 
under the elevated regular rig substructure, thus provid- 
ing mud storage without using additional floor space. 
The additional height of the rig floor also provided room 
for mounting the mud pumps under the drilling engines. 
For further testing, the rig was moved to a Monterey 
Oil Co. lease and used to drill a well duplicating condi 
tions anticipated on the island 
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Fig. 3 — Laying of the submarine pipeline cluster. Four 
pipelines are bound together to form the cluster. 


The island has tanks for fresh water and fuel oil. 
Drilling mud was mixed on shore and taken to the 
island in barges; this was simpler and cheaper than 
separately transporting the dry materials and water. 
Waste mud was barged to shore for disposal. 

Cementing operations are conducted from cement 
trucks on barges; the slurry is pumped through leak- 
proof hoses to connections at the casing head. Well 
logging is also handled from barge-borne equipment. 
Extreme care is exercised to assure that no waste ma- 
terial of any kind, no matter how clean, is permitted to 
get into the ocean. 

Crew changes are made by boat from the Seal Beach 
pier on regular eight-hour shifts. 

The island has emergency housing for the members 
of the drilling crew in the event of water so rough as 
to make their departure from the island inadvisable. An 
office and emergency supply room is located under the 
elevated change room, and contains a desk, a telephone 
tied into the Long Beach exchange by a submarine cable, 
stove, sink, refrigerator, folding bunks for 12 persons 
and emergency food supplies for 12 persons for 48 
hours. 

For fire protection there are power pumps and chemi- 
cal equipment, hoses, fog nozzles, and portable extin- 
guishers. There is also a barge mounted pump and 
monitor. Life jackets are readily available and life rafts 
are placed so that they may be quickly released from 
either the top or bottom of the several escape ladders. 
The safety line from the derrick is attached to the outer 
corner of the wharf. Foghorns and lights are available 
and are used as required by Coast Guard regulations. 


The Submarine Pipeline 


During the drilling and completion of the first six 
wells there were no pipelines between the island and 
the shore, and the wells were cleaned up and tested into 
a barge. The crude oil produced during this comple- 
tion procedure was barged to Wilmington. The reserve 
blocked out by the first six wells was of sufficient mag- 
nitude to economically justify the construction of sub- 
marine pipelines from the island to a separation and 
shipping center on shore. 
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Pipeline construction was started in April, 1955. 1 
work consisted of the installation of one 8-in. and thi 
3-in. somastic coated steel pipes along with one 1|.275-i 
armored electrical communication cable. All were place: 
parallel in a common trench from the drilling islan: 
approximately 8,100 ft offshore to the tank batte: 
approximately 2,500 ft inland. The lines and cable o 
shore were laid parallel on the bottom of a trench. Fo 
the submarine section of the line, the lines were welde: 
on shore into sections approximately 750 ft long, th: 
joints X-rayed, and the assembled lines pressure teste: 
and coated. All three lines and the multiconductor ele« 
tric cable were banded together and then pulled out t 
sea. Another cluster 750 ft long was then added, ar 
this procedure was continued until the line reach 
the island 8,100 ft offshore. The cluster of lines w 
then embedded by means of a hydraulic jetting machine 
The pipe on land and oui io a water depth of 15 ft w 
covered a minimum of 3 ft. Cover from the 15 ft wate 
depth to the island was 5 ft minimum, in order to pri 
tect the line from dragging anchors of large ships. Th 
first oil was shipped through the line in July, 1955 
The commingled oil and gas from the wells is ship 
ped to shore through the 8-in. line, which has a capacit 
of 10,000 BOPD. The oil and gas are separated o 
shore. The individual wells are tested by producin 
through a gas-separator on the island, from which th 
oil is delivered to shore through one of the 3-in. lin« 
to a tank for gauging and the gas is metered out of th 
separator. The two remaining 3-in. lines are availab! 
for delivery of mud, fuel oil, fresh water, or other u 
as required. There is no crude oil storage on the islan: 


Present Modification Program 


The three originally-planned wells from the islan 
were successfully completed. Modifications to the cella 
permitted the drilling of four more wells. Then, by usin 
the mud pit as a cellar, three more locations were poss 
ble. The tenth well is now being drilled from the island 
The experimental island had proven a success, but ha: 
reached its capacity in cellar space to accommodate th 
required wellhead equipment and flow lines. 

The original plans, which provided for the buildin 
of a 600 ft square island, were again studied. Constru 





















Fig. 4 — Construction work is now proceeding on o 
raised floor which will provide additional cellar space. 
This modification of the island will permit the drilling 
of an additional 60 wells. 










tion of such an island would require about two years 
The original island was satisfactory except for the lack 
of the necessary cellar room. It was, therefore, decided 
to modify the original island. 


As the first step, 60 additional conductors were driven 
with approximately 5 ft between conductors. Also, 
added strength and permanence was given to the island 
by removing the sand from within the steel shell of the 
island, replacing it with gravel, and grouting this gravel 
with cement through pipe driven into the gravel to a 
depth of about 20 ft. This was done in a ring 4 ft wide 
adjacent to the outer steel sheath, thus forming a con- 
crete wall inside of the interlocking sheet piling. 

Work is now in progress on the construction of a 
new working floor which will be 9 ft above the forme: 
working surface of the land. The space between the old 
floor and the new floor will provide the cellar space 
required for the drilling of the additional 60 wells. 


Conclusion 


The experience that has been developed through the 
preliminary stage of California’s first real offshore drill- 
ing program is of benefit not only to Monterey Oil Co. 
and The Texas Co., but to the California petroleum in- 
dustry generally. The new California law makes offshore 
exploration more attractive and the other companies 
may be able to profit by the knowledge gained through 
this early development. Although the new California 
law no longer requires that offshore drilling be done 
from filled islands, building of this type of structure 
may still be the most practical under certain conditions. 

The Monterey-Texas interests are pleased to have 
been of service to the industry by helping to establish 
legal precedents concerning offshore drilling and by 
pioneering the building of the first Pacific Coast offshore 
drilling island. wok 
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UPHOLE FORMATIONS TESTED WITH ACCURACY 


Halliburton’s New 
OFF BOTTOM TESTING 
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Another Important Reason Why 


For Your Drill Stem Test! 

























Halliburton’s Best 





When you want to know what’s up hole, Halliburton 
can tell you accurately with its new technique for 
Off Bottom Testing .. . 


* Isolates Any Given Zone Regardless Of Location 
In The Well 


* Can Test Other Zones Without Removing From 
The Hole 


* Eliminates Drill Pipe Anchor From Bottom 
Packer To Well Depth 


* Saves Rigtime Required To Makeup And Break- 
out Anchor Pipe 





The procedure employs three unique Halliburton 
features ...Sidewall Anchor Assembly, Flush Joint 
Pressure Equalizer, and Expanding Shoe Wall 
Packer Assembly. These important tools, added to 
Halliburton’s Hydro-Spring Tester and Bourdon 
Tube, make it possible to test with ease and accu- 
racy the uphole formations that may have been 
previously overlooked. 

So, don’t pass up any zones—one trip in can test 
them all, regardless of location in the well. Talk it 
over now with your Halliburton Testing operator. 
Just phone the local or district office of Halliburton 
Oil Well Cementing Company. 


HALLIBURTON 


TESTING SERVICE 


250 SERVICE CENTERS—JUST MINUTES AWAY FROM ANY RIG 





OVER A BARREL ON © 
EMULSION BREAKING @ 





e@ If emulsion breaking chemicals have you over 
a barrel . . . as to treating cost, handling, or re- 
. check now with your Visco Man. Or, call 


MAdison 3-0433 ... or write Box 6826, Houston. 


sults . . 


We know that a fast, effective Visco formula can 


be fitted to your problems . . 


. and will make every- 
body happy (including us). 


VISCO PRODUCTS COMPANY 
INCORPORATED 
2600 Nottingham at Kirby « Houston 5, Texas 
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McCullough M-3 Guns 
Gain 760 B/D 


Previous Attempt Failed 


To Get Penetration 


Powerful McCullough M-3 Guns fir- 
ing improved Ogival Bullets again 
proved their superior penetrating 
power in a recent job in this Texas well. 

The well had first been perforated 
by another service company, but the 
operator was unable to break the 
formation down under 3500 pounds 
pump pressure. 

Depth of operation was 9800’. 5%” 
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ow 


Deep penetrating, accurately 

placed shots give you the results 

you want... better production. 
Always specify McCullough Perforating 
Service for MORE OIL! 








O.D., 17 lb. casing had been cemented 
through the pay zone. 

McCullough 4%” O_D., 30-shot M-3 
Perforators fired 84 improved 4” Ogival 
Bullets in the same interval previously 
perforated. The formation broke down 
at 1200 pounds pump pressure. After 
treating, the well made 760 barrels of 
oil per day—the best well in the field. 
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Literature Available on 
McCullough M-3 
Bullet Perforators 


improved Ogival Bullets. 


mented by actual job reports. 


Calif. 








1. Four-page illustrated bulletin, No. 
203, gives advantages, characteris- 
tics, etc. of M-3 Bullet Perforators and 


. Sixteen-page pocket brochure tells 
benefits of M-3 Perforators, supple- 


Write to McCullough Tool Company, 
5820 S. Alameda St., Los Angeles 58, 











OMPANY 


Cable Address 





M-3 Guns Shoot 400 Holes in 
73° Angle Well in Two Runs 


Unusual Perforating Job Could Only Be Accomplished 
By Simultaneous Firing, Deep Penetrating 
McCullough M-3 Bullet Perforators 


While unusual jobs such as this occur only infrequently, they serve to 
demonstrate the ability of McCullough Perforating Services to accomplish the 
tough jobs with the same speed, economy and satisfaction that characterize 


This well was directionally drilled from an onshore location out under the 
ocean bed. Due to the high degree angle of the well (73°) it was impossible 
to perforate by regular electric wire line methods. 


8%” O.D., 36 Ib. casing had been 
cemented in 12%” hole. The interval to 
be perforated was from 6696’ to 6796’. 

Ten 4%” O.D., 20-shot McCullough 
M-3 Guns (50°) were screwed together 
and run in on 4%” F. H. drill pipe. A 
swab cup was arranged on a 1%” O.D. 
probe, attached to the electric wire 
line and pumped down the drill pipe 
until it contacted a firing sub attached 
to the top of the guns. 

Two sets of M-3 Guns of the above 
length were run, firing 400 k” 
proved Ogival Bullets in the 100’ zone. 
The entire job took only 26 hours to 


im- 


complete which included two round 
trips with the drill pipe and two runs 
with the probe. No trouble, no delay 
was experienced. The job was satisfac- 
tory in every respect. 

Only McCullough M-3 Guns, firing 
all shots at one time, could successfully 
complete this unusual job. In addition, 
deep penetration through heavy casing 
cemented in large diameter hole was 
assured by the extra power of M-3 
Guns firing improved Ogival Bullets 
hardest shooting bullet perforators in 


the world 


LOS ANGELES 
HOUSTON 


MACTOOI EDMONTON 

















GUIBERSON 





_ WITH THE 


TYPE’ 


Finest Tubing 
Swab Cup Made 
For Lifting 


Heavy or 
Light Loads 





When you're swabbing 

from any depth or want to 
unload the well completely, 
you can do a better, more 
efficient job with the 
Guiberson Type L on the 
the regular Guiberson 

Swab. This unbeatable 
combination drops fast and 
seals perfectly for the lift — 
and has the flexibility and 
durability to do a top-notch 
job. Capable of lifting 

heavy loads, this versatile 
cup will also lift the last light 
loads from the bottom of 

the well. The L lifts full loads, 
not just part way but all 

the way to the surface, 
slipping by collars just as 
though they weren’t there. 
The next time you need a 

a general purpose swab cup, 
get the Type L. Look for 

it at any supply store, 
packaged in pairs in the 
brilliant blue box. 








when a blowout threatens... 


BE OUT OF THE WAY 


BUT IN CONTROL! 














See pages 2555-8 of 
Composite Catalog, 
21st Edition. 








Here’s the most effective combination 
of equipment ever designed for remote 
control of high-pressure wells. When 
a blowout threatens, you can depend 
on positive shutoff by simply moving 
one lever on the master control mani- 
fold located at a safe, remote distance 
from the danger zone, or by pushing a 
button at the driller’s control panel. 
The HyYpDRIL System of Blowout Pre- 
vention Control can save lives and 
equipment in an emergency. 

Contact your nearest HYDRIL repre- 


sentative today to get the full story. 





































































































HYDRIL COMPANY 


714 W. OLYMPIC BLVD., LOS ANGELES 15, CALIF 
FACTORIES AT 


LOS ANGELES; HOUSTON, TEXAS; YOUNGSTOWN, OHIO; ROCHESTER, PA. 


HYDRIL TAKES THE BLOW OUT OF BLOWOUTS 


HYDRIL SALES OFFICES 
California: Bakersfield, Los Angeles, V 
Louisiana: Harvey, New Iberia, New Orle 
Oklahoma: Tulsa; Pennsylvania: Rochester 
Texas: Corpus Christi, Dallas, He 

Odessa; Wyoming: Casper 
Canada: Calgary, Edmonto 








| 
| 

















JO& B. ALFORD 


Visit from Institute President 


A long standing custom in AIME is for the Institute 
president to visit as many local sections and student 
chapters as possible during his term of office. Following 
this custom the Petroleum Branch was favored by a visit 
from our current president, H. DeWitt Smith, to sev- 
eral sections during the week of Nov. 7-11. His talks 
before meeting of the sections visited and informal con- 
ferences with members were highly beneficial to the 
Petroleum Branch members in acquiring a better under- 
standing of AIME operation and organization. 

The tour included visits with the Mid-Continent Lo- 
cal Section, the Kansas Local Section, the Hugoton and 
Panhandle Sections in a joint meeting, the Texas Tech- 
nological College AIME Student Chapter, the South 
Plains Section, and an informal conference with some 
members of the Dallas Local Section. The writer accom- 
panied President Smith throughout the trip. Robert B. 
Gilmore, Petroleum Branch chairman, and William D. 
Owsley, member of the Branch Executive Committee 
and chairman of the Branch Advertising Committee, 
made a portion of the tour. 


Changing Membership Ratio 

The talks made by President Smith before the sections 
visited were concerned with the possible future organi- 
zation and operation of AIME. As shown in Table 1 
the ratio of AIME membership in the three Branches 
is now undergoing a considerable change and will prob- 
ably continue to do so in the next several years. (For 
example, the 1955 figures shown in Table 1 were on 
Jan. 1, and as of this writing the percentages are: 
Mining 43.34; Metals 24.54; Petroleum 31.71.) This 
changing ratio of membership indicates the need to 
study Institute organization for the benefit of all three 
Branches and their various professional divisions. 


TABLE 1 — CHANGING RATIO OF AIME MEMBERSHIP 
By Percentage in Branches 


1946 =——-1950 1955 1960 1963 
Mien! fs. +a 52 45 39 36 
SA ee ss, Sota 25 25 27 28 
i 23 30 34 36 


Accordingly, the Institute Board of Directors, in Feb., 
1955, authorized the appointment of a Long Range 
Planning Committee. Subsequently President Smith ap- 





? 


PANHANDLE SECTION OFFICERS—Back row, left to 
right: Donald E. Radtke, director and past chairman; 
William D. Owsley, member of Petroleum Branch Exec- 
utive Committee; Don W. Jackson, secretary-treasurer. 
Front row: Henry J. Rose, chairman; AIME President 
Smith. 
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REPORT and INTERPRETATION 





EXECUTIVE SECRETARY 
PETROLEUM BRANCH, AIME 


at 


> “— 
: 


HUGOTON SECTION OFFICERS — Back row, left to 
right: Clifford R. Horn, section delegate; Claud P. 
Pyeatte, director; AIME President Smith. Front row: 
R. C. Dickson, secretary-treasurer; William H. Flood, Jr., 
director. 


pointed Andrew Fletcher, past president for the Mining 
Branch, Leo F. Reinartz, past president for the Metals 
Branch, and Carl E. Reistle, Jr., president-elect for the 
Petroleum Branch, as members of this Committee. 

Smith reported that during the months since Feb- 
ruary the Committee has made good preliminary prog- 
ress on the study of Institute organization. Two basic 
concepts of organization have been explored: (1) To 
retain AIME as one central organization with three sub- 
sidiary branches or organizational units; (2) To re 
organize AIME as the American Institute of Mineral 
Engineers with three autonomous member institutes in 
the fields of Mining, Metals and Petroleum. 


Opinions on Reorganization 


Preliminary memoranda on these proposals have been 
circulated to various groups in the Institute and a good 
many individual members, soliciting their opinions. 
Smith reported that no enthusiasm whatever could be 
found in any group for the name “American Institute 
of Mineral Engineers.” These memoranda also de 
veloped strong opinion in all three Branches that the 
AIME should be basically maintained as an over-all pro- 
fessional society for, engineering in the mineral indus- 
tries in the United States. 


President Smith stated further that the Long Range 
Planning Committee will continue its study of Institute 
organization and will possibly employ an outside con- 
sultant to assist in this work. He then expressed his 
personal views that the AIME Board of Directors should 
be equally distributed between the three Branches and 
that the Institute presidency should be rotated between 
the Branches on an equal basis. Further, each Branch 
should be free to set up its own local sections along 
professional lines independent of the sections of the 
other Branches. Also, the Institute initiation fee 
should be reduced but dues should be at the same 
figure for all three Branches and not reduced in the 
immediate future. Each Branch should have the use of 
any surplus income after setting out a certain stipulated 
amount for the operation of the Institute and reasonable 
financial reserves. 

We hereby express appreciation to President Smith 
for his visit and the better understanding of Institute 
organization brought about by his discussion. week 












Eighteen Paper Meet 
Planned at Texas A&M 


Eighteen technical papers will be 
presented at a meeting to be held at 
A&M College of Texas on Jan. 9 
and 10. Sponsors of the meeting will 
be the AIME Student Chapters of 
the University of Texas and A&M 
College of Texas. The meeting will 
be held at the A&M Memorial Stu- 
dent Center, College Station, Tex. 


THE PROGRAM 


Evaluation of Fivid Input Wells by Pressure 
Build-up Analysis, A. R. Groeneman, Sinclair Oil 
and Gas Co Water and Gas Coning in Wells 
: Completed in Homogeneous Sands, K. E. Montague 
a % * and M. D. Noble, Sun Oj! Co.; Studies of Pres 
is a Comparison of Conventional Pumps soe, Sabian & Dimiad Seneiiie of dinate 

State, C. S. Matthews and H. C. Lefkovits, Shel! 


“end the Double-Displacement Rod Pump. Development Co.; and Comparison of Pressure Dis 


tribution During Depletion of Tilted and Horizontal! 


All pumps compared at twelve 20” S/M. Reservoirs, D. §. Howard, Jr., and H. H. Rachford 


jr., Humble Oil and Refining Co 





| 
A Reservoir Study of the Greater Clareton Area 


Wyoming, Core Laboratories, Inc.; Hydraulic Frac 
turing as Applied to the Chalk Rock, Caddo-Pine 
Island Field, Caddo Parish, Louisiana, B. F. Pat 
terson, The Texas Co Sampling and Analysis Pro 
cedures for O Reservoir Fluids, Frank O. Reude!l 


; T sure Distribution and Gas Mass in Pipe Lines, D 
and ©. Atkinson and Paul B. Crawford, TPRC, A&M 
| a ¥ t IMP—63.12 B/D College; and Secondary Recovery Potentialities of 
: ' eB 1 taco Buckrange Reservoir, Stevens Oil Field, Columbia 


: Nevado, and Ouachita Counties, Arkansas, Pav! 
. - TUBI re ; ‘ - Meadows ond Murphy Hawkins, U.S. Bureau of 
a! B/t 


; 


40 50 60 10 huber, Core Laboratories, Inc.; Generalized Pres 
- 


Mine 


Determination of Fluid Saturation in Porous 
Media with a Radiocobalt Tracer and Application 
to Petroleum Production Research, Leonard B. Lip 
son and L. S. Gowrnay, University of Houston 
Radicactive Isotopes in Petroleum Production Re 
search, J. W Watkins, H. WN. Dunning, F. E 
Armstrong, R. J. Heemnere, Lun Hsiao, U.S. Bureau 
of Mines; Scale Limitations in Potentiometric Mode! 
Consiruction, A. S. Odeh, H. B. Bradley, and J 
P. Heller, Magnolia Petroleum Co.; Bubble Forma 
tion in Petroleum Reservoirs, Denton R. Wieland 
ond Horvey T. Kennedy, A&M College of Texas; 
and Interfacia Tension of Water and Methane 
with Minor Components, E. F. Newman and E. W 
Hough, University of Texa 

Cloud Point Determination of Crude Oils under 
Dymanic Conditions, James E. Howell and Frank 
Jessen, University of Texas; The Significance of 
Mud Viscosity, W. C. Goins, Gulf Oil Corp.; A 
Lightweight, Low Woater-Loss, Oil Emulsion Cement 
for Use in O Wells, G. K. Dumbould, F. A 
Brooks, Jr B. E. Morgan, and G. W. Binkley 
Humble Oj! & Refining Co.; Dexet, a New O 
Well Cementing Material, Dowell, Inc 


CORRECTION 


Iwo errors were made in Tech- 
nical Note 316, “Fluid Mapper 
Model Studies of Mobility Ratio,” 


2 a! rUBII > UM —210.24 B/D published on page 49 in the Nov., 


1955, issue of JOURNAL OF PETRO- 
2%” Double-Displacement ROD PUMP — 286.56 B/D re sigh bmw he aug tn se 0 
. tion of the paper, the equations 
should have read 
The mobility ratio, M, 1s 
defined as 


HARBISON-FISCHER (displaced fluid ) 


M 
: (displacing fluid) 
Double -Displacemeant where k effective permeability and 


ROD PUMP eles 


“When the displacing fluid is mis- 

; P cible with the displaced fluid, the 
All of the production advantages of a full-volume tubing as tate ace maa ee 
pump PLUS the operating advantage of an insert pump bility ratio may be expressed as a 
ratio of viscosities 


ae "Best Pumps in the Oil Pate " ' . (displacing fluid) ,, 


OR mengigiaetan 


(displaced fluid) 


270° 280 290 300 
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WILLIAM F. SHELTON has joined 
Creole Petroleum Corp. He grad- 
uated with a degree in petroleum 
engineering from the Colorado 
School of Mines. Shelton is a senior 
petroleum engineer with Creole in 
the San Joaquin District of the East- 
ern Div. 

DUES BILLS IN MAIL 

Pursuant to Article II, Section 2, 
of the by-laws of the AIME, notice is 
hereby given that dues for the year 
1956 are payable Jan. 1, 1956, as fol- 
lows: Members and Associate Mem- 
bers, $20; Junior Members for the 
first six years of Junior Membership, 
$12, and thereafter, $17; Student 
Associates (including an annual sub- 
scription to a monthly journal), 
$4.50. 
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P E RSON A LS of a since sere 


THICKENED OIL 
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GaBor Dessau is now at the Tech- 
nion (Hebrew Institute of Technol- 
ogy) in Haifa, Israel. He will be 
there for one year as UNESCO ad- 
visor on the setting up of a geological 
engineering and mining department. 
After retiring in Sept., 1953, from 
Italian government service, he has 
continued as charge’ de cours of min- 
eral deposits at the University of 
Naples. From its inception, he has 
been in charge of the Uranium Pros- 
pecting Centre (Centro Ricerche 
Minerali Radioattivi) of the ENI- 
AGIP group. 

Dowell Incorporated has an- 
nounced changes involving three 
AIME members. Ear J. Davis has 
been named assistant to the general 
manager and transferred to the Tulsa 
general offices of the company. He 
formerly was area manager in 
Dowell’s Illinois District. JaMes R. 
PAUL, chairman of the AIME Kan- 
sas Local Section, has been trans- 


ferred to Carmi, Ill., from Wichita, 
Kans., to fill the vacancy created by 
the promotion of Davis. S. R. WEsT 
will replace Paul as district engineer 





NE RN I TT OS RY IT 


in Kansas for Dowell. He was for 
merly located in Great Bend. 


Ww 


RODNEY S. DuRkKEE, chairman of 
the board of Lane-Wells and a mem 
ber of the executive committee of 
Dresser Industries, Inc., has beer 
elected president of the California 
Manufacturers Association for the 
1955-56 fiscal year. Lane-Wells has 
also announced changes in positions 
of other AIME members. JaMes D 
HUGHES, formerly vice-president-sales 
will now direct the sales and servic« 
activities of the company. JoHN ( 
Stick, JRr., formerly the company’: 
chief geophysical engineer, has been 
appointed chief engineer in the ab 
sence of V. L. ForsyTH, who is tem 
porarily on special assignment as ad 
ministrative assistant to the president 
J. SHANNON BAKER has assumed the 
duties of the chief geophysical en 
gineer. J. M. HENDERSON, Gulf Coast 
Div. Manager, has been promoted 
to manager sales and service for the 
company. R. E. BusH has been ad 
vanced to manager of the Gulf Coas! 
Div. 


SERVICING 


OlLa GAS 
WELLS 





T-O (Thickened Oil) is The Chemical Process Co 
pany’s new low-cost fracturing fluid. Made 
location from crude oil, kerosene, or Diese! o 
T-O is ready for use within 30 minutes after mixi 
has begun. T-O has. 


HIGH VISCOSITYvvvyw 


for good sand carrying ability 


LOW FLUID LOSSvvy 


for deeper and better fractures 


INTERNAL BREAKERwW~y 


for faster, more complete clean-up 


LOW FRICTION LOSSw 


for lower injection pressure 


GOOD COMPATIBILITY 
with demulsifying agents or additive 
such as Adomite to lower fiuid loss 


Fracturing treatments with T-O are givi! 
excellent results. A gas well in Hansford Count 
Texas, increased from a potential of 9,000 
68,000 MCF after treatment with 20,000 gallo: 
of T-O (made from kerosene). T-O approaches t! 
ideal fracturing fluid. Your nearest Chemical Pro 
ess Company office will be glad to give you add 
tional information, or write to us at P. O. Box 83 
in Breckenridge, Texas. 


The Chemical Process Company 


There is no obligatior 
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A picture and biographical sketch 
of HENRY H. RACHFORD, JR., were 
published in the November, 1955, 
issue Of JOURNAL OF PETROLEUM 
TECHNOLOGY 


A. R. HENDRICKSON was gradu- 
ated from Beloit College, Wisconsin, 
in 1949 with a BS degree in organic 
LATCH-ON CENTRALIZERS chemistry. He did post graduate ; 4 
work at the University of New W. P. SCHULTZ received a BS de- 
with the NEW : al , ‘ 
Hampshire, receiving his MS degree gree _in chemical engineering from 
KON-KAVE BOW in 1950. He joined Dowell Incor- Rice Institute in 1947. Until 1955 
porated, Tulsa, as a chemist shortly he was employ ed by Humble Oil & 
thereafter. Refining Co. He is now employed 
by Core Laboratories, Inc., and is a 
member of the Engineering and Con- 
sulting Dept. in Dallas. 


R. A. SEBAN is an associate pro- 
fessor of mechanical engineering at 
the University of Colifornia, Berk- 
eley. He holds BS and PhD degrees 
received from the university in 1938 

and 1948, respectively. His research 
Son thee ee ee me aud teaching experience primarily 
cross - section, highest cerns heat transfer. 


ay ME ag G. L. Foster holds a BS from 


a Ottawa (Kansas) University and a 
— in’ forraing and MS degree from Kansas State Col- 
KAVE BOWS. . lege. He was employed by Dowell 
Incorporated as a chemist in 1952, 

FIRST and was promoted to his present po- 
sition of laboratory group leader at 


IN THE Tulsa last year. 
FIELD! 


@ Strongest — greatest FRANK G. MILLER is professor and 
; ide thrust 2 head of petroleum engineering at 
Serer ee - % Stanford University. He holds MS 
and PhD degrees awarded by the 
University of California, specializing 
in oil-reservoir engineering. 


@ Most resistance to 
deformation ... 


@ Easiest to install and run. 


R. B. Rosene is laboratory direc- 
tor for Dowell Incorporated at Tulsa. 
He received his BS degree from 


RB and W in C Whittier College in 1944 and his doc- 
° torate from Purdue University in 


Well Completion Specialists 1951. First serving as research chem- 
WEST COAST ist for the Plantation Div., U. S. 
arene, Rubber Co. in Sumatra, he joined 
s Angeles 54, Co Dowell as a research chemist in L. C. VoGEL graduated 
Phone DA-4-1106 a : p 
1952. UCLA in 1949 with a degree in 
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AUTOMATIC, REGULATED 
FILL OF THE CASING 
AS IT IS BEING RUN 





































ECONOMICAL 


Baker Flexiflow Fill-Up Equipment not 
only permits the casing to be run more 
rapidly with greatly reduced pressure 
increases or surges on the formations 
below the shoe, but in addition, by 
eliminating the surface-filling operation, 
saves an additional amount of rig time. 
When these advantages are considered, 
Flexiflow Equipment is much more 
economical to run than conventional 
floating equipment. 







XpROTECTION AGAINST 
DESTRUCTIVE 
HIGH-PRESSURE SURGES 









BAKER 


LEXIFLOW 


FILL-UP SHOES AND COLLARS 


Baker Flexiflow Fill-Up Equipment, 
a new companion line to Baker Differ- 
ential Fill-Up Equipment is designed to 
be run in wells of moderate depth. 


FEATURES 


SELF-ADJUSTING FLOW 
CONTROL DIAPHRAGM 


Baker Flexiflow Fill-Up Equipment 
contains a variable valve (Self-Adjusting 
Flow Control Diaphragm) that auto- 
matically opens wider or closes down 
to provide MINIMUM BOTTOM-HOLE 
PRESSURE INCREASES and regulated fill 
of the casing. 


POSITIVE BALL-TYPE BACK- 
PRESSURE VALVE FOR CEMENTING 


A Ball-Type, Back-Pressure Valve for 
floating or cementing can be brought 
into action at any time by applying a 
pressure of approximately 700 psi to 
the casing. See Figs. 2 and 3 below. 


COMPLETELY DRILLABLE 


The internal parts of Flexiflow Fill-Up 
Shoes and Collars consist of readily 
drillable concrete, a small amount of 
Bakelite and synthetic rubber. 


FIG. 2 

Ball 
pumped 
through 
diaphragm 





el lle surges over 


5,060 psi can be created 
each time a joint of 
casing equipped with 
conventional floating 
equipment is lowered 
under certain slim hole 
clearances. 


High-Pressure surges can 
fracture and weaken a 
formation sufficiently to 
cause lost circulation or 
the failure of a primary 
cement job. 


These destructive, high- 
pressure surges can be 
reduced... 


by increasing the lower- 
ing time per joint 


by increasing the casing 
clearance 

or 
by running Baker Filexi- 


flow (or Differential) Fill- 
Up Shoes and Collars 


NOTE: 


Baker Fill-Up Equipment 
contains a variable valve 
that automatically opens 
up and adjusts to main- 
tain minimum boftom- 
hole pressure as the 
casing is being run. 














BAKERS 


OIL Ma INC. 


ANGELES-- NEW 


YORK © 


POLICIES FOR 


Pursuant to Article X of the by- 
laws of the AIME, the following in- 
formation is hereby given as to the 
“conditions, prices, and terms under 
which the various classes of mem- 
bers, and Student Associates, sever- 
ally, shall be privileged to obtain 
publications of the Institute during 
the ensuing year.” 

Publications authorized for issue 
in 1956 include the following: Min- 
ing Engineering, published monthly, 
containing material, including tech- 
nical papers, of interest to those en- 
gaged in exploration, mining geol- 
ogy and geophysics; metal, nonme- 
tallic, and coal mining and beneficia- 
tion; and fuel technology. The Jour- 
nal of Metals, published monthly, 
containing material, including tech- 
nical papers, of interest to those en- 


n makes it 


gaged in nonferrous smelting and 
refining, iron and steel production, 
and physical metallurgy. The Jour- 
NAL OF PETROLEUM TECHNOLOGY, 
published monthly in Dallas, con- 
taining material, including technical 
papers, of interest to those engaged 
in petroleum and natural gas drill- 
ing and production. 

Annual subscriptions to any one of 
the above journals will be provided 
all members in good standing with- 
out further charge, a subscription 
credit of $6 for members and $4 for 
Students Associates being included in 
the dues paid. (A member ceases to 
be in good standing if current dues 
are not paid by April 1). If more 
than one of the monthly journals is 
requested, $4 extra will be charged 
for an annual subscription, or 75c 


Competition is the final answer in any field of 
activity. Competition is the dynamo that creates 
fresh new ideas, products and services. 
Competition makes oil well cementing service 
better, too. On your next job, call in the company 
that put quality competition in the business 

with its new engineered approach to cementing. 


Engineered answers to | 


oil well cementing 


BJ SERVICE, INC. 


For further details see page 32 


JOURNAL 


1956 AIME PUBLICATIONS ESTABLISHED 


for single copies of regular issues 
and $1.50 for special issues. The non- 
member subscription price for each 
journal is $8 in the Americas and 
U.S. possessions; foreign, $10; for 
single issues 75c in the Americas; 
$1.00 foreign for regular issues; $1.50 
for special issues. Student Associates 
will be entitled to the same privileges 
for all publications as members. 
AIME members subscribing to more 
than one of each of the three month- 
ly journals will be billed at the non- 
member rate of $8 per year, domes- 
tic; $10 foreign, for the extra sub- 
scription(s) 

Three volumes of “Transactions’’ 
are authorized for 1956 publication, 
as follows: No. 202, Mining Branch; 
No. 203, Metals Branch; and No. 204, 
Petroleum Branch. Volumes 202 and 
204 will be available to members at 
$3.50 each for a first copy if paid for 
in advance with dues; otherwise at 
the nonmember rate of $7 less 30 
per cent. Nonmembers $7 in the 
United States; Foreign $7.50. Vol- 
ume 203 will be available to mem- 
bers at $4.50 each for a first copy if 
paid for in advance with dues; other- 
wise at the nonmember rate of $9 
less 30 per cent. Nonmembers $9 in 
the United States, foreign $9.75. 

Special volumes now planned for 
publication in 1956 include the fol- 
lowing: Open Hearth Proceedings, 
Vol. 39, price to AIME members $7; 
nonmembers, $10. Blast Furnace, 
“oke Oven, and Raw Materials Pro- 
ceedings, Vol. 15, AIME members 
$7; nonmembers $10. Electric Fur- 
nace Steel Proceedings, Vol. 13, 
AIME members $7; nonmembers 
$10. Statistics of Oil and Gas Devel- 
opment and Production, Vol. 9, cov- 
ering data for the year 1954, mem- 
bers $5, nonmembers $10. Proceed- 
ings of Nigara Falls Conference on 
Reactive Metals tentative; price to 
be determined 

If dues are paid subsequent to 
January 31, back issues of Institute 
publications will be supplied only 
if adequate stocks are on hand. A 
member is not entitled to receive a 
or a spe- 
cial volume, in lieu of a monthly 
journal, free of charge on member- 
ship. Members in arrears for dues 
are not entitled to special members’ 
prices for publications. 

Rocky Mountain Members may 
have their choice of an annual sub- 
scription to one of the monthly jour- 
ke 


volume of “7ransactions,” 


nals on request 
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EMPLOYMENT 





NOTICES 





* 





The JOURNAL OF PETROLEUM 
TECHNOLOGY will post notices of 
men and jobs available. Companies 
and AIME members are invited to 
use this space, for which there is no 
charge. Except as noted, address re- 
plies to: Code (appropriate number), 
JOURNAL OF PETROLEUM TECHNOL- 
ocy, 800 Fidelity Union Bldg., Dal- 
las 1. These replies will be forwarded 
unopened and no fees are involved. 
Replies for notice listed by code 
W2555(a) are to be sent to: Engi- 
neering Societies Personnel Service, 
8 West 40th St., Room 1207, New 
York 18, N. Y. The ESPS collects 
a fee from applicants actually placed. 


POSITIONS OPEN 


yy Research mathematician or theo- 
retical physicist. Major oil company 
research laboratory has opening for 
PhD or MS in the field of reservoir 
engineering research. Previous expe- 
rience in applied mathematical re- 
search desirable but not necessary. 
Excellent working environment with 
opportunities for advancement. Re- 
plies confidential. Send detailed edu- 
cation and experience record, age, 
and state salary expected. Code 590. 


yy Combination Petroleum-Geolog- 
ical Engineer or Petroleum Engineer 
competent in geology, with one to 
two years’ experience; Gulf Coast ex- 
perience preferred. Reserves and 
evaluation work for expanding gas 
transmission company located in 
Houston. Reply stating experience, 
desired salary and approximate avail- 
ability date. Code 594. 


WANTED 


Petroleum Engineer be- 
tween 28 and 35 years of 
age with field experience 
and interest in reserves and 
valuation work. New York 
financial institution offers 
opportunity for experience 
and advancement. Please 
send resume of personal 
history, experience, and re- 
cent photograph. Replies 
will be considered strictly 
confidential. Write P. O. 
Box 226, Church Street Sta- 
tion, New York 8, N. Y., 
Dept. 12. 
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yy Petroleum Engineer, four or five 
years’ experience with emphasis on 
valuation and reservoir engineering. 
Must be free to travel with Dallas as 
headquarters. Initial letter application 
should set out past experience in de- 
tail, education and salary require- 
ments. Code 595. 


yy Professor of Petroleum Engineer- 
ing to be responsible for undergradu- 


ate curriculum and graduate courses 
MS or PhD. Interest in research de 
sired and full cooperation of univer 
sity assured. Rank and salary depend 
ent upon education and experience 
Code 596. 


yy Assistant or Associate Professor o! 

Mining Engineering to teach courses 

in mining and petroleum engineering 
(Continued on Page 37) 





This revolutionary process — 
thoroughly field-proven — 
combines heat and pressure 
to weld together 
thick-walled products. 











drill collars and kellys renewed 


the economical way by 
BAASH-ROSS PRESSURE 


* Produces a perfect weld—one as homogene« 
strong as the base metal itself. 


* No metal added—nor are the physica! p: 
of the welded product changed in any way 


* Use it to “stub out”’ short drill collars for ad 
service ...to renew kellys by reversing thei 
faces... and for many other operations 


Literature available. Or ask your Baash-Ross representative 


BAASH-ROSS TOOL COMPANY 


DIVISION OF JOY MANUFACTURING COMPANY 
Los Angeles, Calif. « Houston, Texas « Odessa, Texas * OBlahoma | 
Casper, Wyo. « Olney, Ill. « Edmonton, Alb., Canada « New York City 
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with a 1%" LTANEOU 
RADIATION LOGGING INSTRUMENT 
FOR 2” TUBING COMPLETIONS 





Kemembler.. 


PGAC—AND ONLY PGAC—GIVES 
YOU A SIMULTANEOUS GAMMA 
RAY and NEUTRON/NEUTRON 
LOG..:ON A SINGLE TRIP IN THE HOLE! 


WHY USE 2 SEPARATE TOOLS 
REQUIRING 2 TRIPS IN THE HOLE? 


® All of the superior logging standards of calibration, sta- shneeh OGRE salle tie maaatien % loner 
bility and absence-of-drift — for which operators have hailed wathdes Bm gp rae ng 
PGAC’s. 352” logging instrument for years —have been fully 

retained in this new 1's Simultaneous Radiation Logging 

Instrument now available for tubing completions. 


ak ne 
—— 


® PGAC’s Simultaneous Radiation Log is still the only Gamma 
Ray and Neutron/Nevtron Log available. It utilizes a pure 
Neutron/Neutron to avoid the faulty phenomena of detrimental 
gamma rays inherent in all other “neutron” logs. 


® Moreover, PGAC originated and was the first company — 
years before any other —to deliver immediately in the field a 
simultaneously recorded radiation log. 


® Furthermore, PGAC has successfully logged and perforated 
the world’s deepest (17,570 ft.) and hottest (349°F.) wells. Such 
know-how leads to PGAC’s constantly improved instrumenis — 
tools without equal in accuracy, safety, operating speed and ; 

4 Y Y P itil PGAC Simultaneous Radiation Log being recorded 


efficiency. on instrument in field unit. Complete log is 
immediately delivered in field to operator 


© .. PERFORATING GUNS ATLAS CORPORATION 
Oo Generel Offices ond Mein Plot: 7790 Seott' Street — Selec Offtes: Anan 


CALL FOR PROMPT SERVICE — ALWAYS READY TO SERVE YOU 
Cy TEXAS: Abilene — Alice — Beaumont — Colorado City — Corpus Christi— Dallas — Fort Worth — Gainesville — Graham 


OIL 
FIELD 
SERVICES 


Houston — Longview — Odessa — Pampa — Victoria — Wichita Falls 
LOUISIANA: Houma — Lafayette — Lake Charles — Shreveport. KANSAS: Great Bend — Liberal 
OKLAHOMA: Ardmore — Healdton — Oklahoma City — Pauls Valley — Perry — Seminole — Tulsa. NEW MEXICO: Hobbs 


AFFILIATE COMPANIES: CANADA — Perforating Guns of Canada, Ltd; Edmonton, Alberta 
GERMANY — Atlas Deutsch- Amerikanische Olfelddienst G. m. b. H.; Kiel VENEZUELA — Servicios Tecnicos Atlas, C.A., Caracas 















































































President-Elect Reistle 
Speaks at Arizona Meet 


AIME President-Elect Carl E 
Reistle, Jr., was the featured speaker 
before some 700 dinner guests at- : , 
tending the annual Arizona Section Nik —firs} /n Pressure Lontro/ 
meeting in Tucson, Ariz., on Dec. 5 
The Humble Oil and Refining Co 
vice-president spoke on “Offshore 
Developments in the Oil Industry 
Reistle also spoke on the morning 
of Dec. 5 to the AIME Student 
Chapter at the University of Arizona, 
telling the students what manage 
ment looks for when hiring engi- 
neers. During both the morning and 
evening talks, President-Elect Reistle 
briefly outlined recent developments 
in AIME affairs. \ 
Other guests representing the Pe 
troleum Branch at the annual An 
zona Section meeting were Branch 
Chairman Robert B. Gilmore, De 
Golyer and MacNaughton; Joe B 
Alford, executive secretary; aid R 
W. Taylor, editor of JOURNAL of! 
PETROLEUM TECHNOLOGY. 
USC Chapter Activities 
The University of Southern Cali- 
fornia Student Chapter of AIME re- 
cently took a field trip to the Union 
Oil Co. water flood at Dominguez 
Oil field. Several other field trips are 
planned for this chapter with 24 ac 
tive members. 
The chapter holds monthly meet 
ings, featuring speakers from indus- 
try. The highlight of the year will be 
the annual Spring Meeting on the 
USC campus, for the students and 
men from industry. 

| Officers for the 1955-56 schoo! 
year are: Leland Cecil, president 
Allen Nelson, vice-president; Bruce 
Melkonian, corresponding secretary 
George Heckler, recording secretary 
and R. L. Perrine, California Re 
search Corp., counselor. 





Student Chapter Reorganized 
The Missouri School of Mines 
AIME Student Chapter has been re 
organized and now consists of three 
active sections, Mining, Petroleum, 
d and Metallurgy. Each section has its 
own slate of officers and own pro 
| grams. Joint meetings of the sections 
featuring discussions of subjects ot 
mutual interest are also held. Until! 
this reorganization the petroleum 
students had their own petroleun 
club, with mining students dominat 
ing the AIME Student Chapter 
Petroleum Section officers are: | 
W. Ware, chairman; B. Malson, vice 
chairman; J. Lutz, secretary; and 
Henry Lamb, treasurer. Prof. J. Pau! 
Roston is faculty section advisor 





DECEMBER, 1955 


Bod -built 


to give you a 


Betier Jo 


111-T High Pressure Cementing Unit 


' OTHER BJ-ENGINEERED 
OIL FIELD SERVICES: 


CALIPER SURVEYS 


» DOUBLE INTEGRATOR 
VOL‘  E COMPUTATION 


TEV ~~ RATURE SURVEYS 


SIMULTANEOUS 
RADIATION LOGGING 


PERFORATING... JET AND BULLET 
ACIDIZING . 
FRACTURING 


; 


ONE OF THE BIG reasons why BJ Service, Inc., 
has become the dynamic new source for oil well 
cementing is BJ-engineered equipment. These 
BJ designed and built units have set new stand- 
ards in high-powered, high-efficiency cement- 
ing. The BJ Service PL-7—built in 1946—was 
the first power unit capable of cementing wells 
at pressures to 10,000 lbs. Today’s BJ equip- 
ment continues to be the pacemaker in the 
industry. 

The result—when you order a BJ job-engi- 
neered cementing service, you know that the 
equipment assigned will be the best in the 
oilpatch. 

Next time, call for the service that gives 
you engineered answers to specified needs... 


call for BJ Service, Inc 


PL-7 High Pressure Cen 


enting Unit 


rd 


BJ Service.Ine. 


GENERAL OFFICES: 
6505 Paramount Bivd., Long Beach 5, Calif. 
Los Angeles phone: NEvada 6-2493 
Long Beach phone: 20-1411 


SERVICE STATIONS THROUGHOUT THE 
PACIFIC COAST, ROCKY MOUNTAIN AND 
PERMIAN BASIN AREAS 


Service 


SUBSIDIARY OF BORG-WARNER CORP 
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Feb. 19-23 in New York City 





PETROLEUM PROGRAM 


Six papers have been definitely 
scheduled for presentation at Petro- 
leum technical sessions at the AIME 
Annual Meeting in New York. Ad- 
ditional papers are to be added to 
the program, which will be held at 
the Statler Hotel during the course 
of the Annual Meeting, Feb. 19-23 


Petroleum Papers 


Names of papers and authors 
scheduled are: 
Problems Related to Negotiation of 
Foreign Oil Concessions, by Burney 
Brady, Continental Oil Co. 
Gilsonite and Its Relation to Oil, by 
Park L. Morse or R. E. Nelson, 
American Gilsonite Co. 
The Importance of Foreign Oil to 
Our Military Services, by Admiral 





HAVE YOU CHANGED 
YOUR ADDRESS? 


In order that publications and cor- 
respondence may reach you promptly, 
the Petroleum Branch, AIME, should 
be advised as soon as possible of any 
change in your address, preferably a 
month before the change becomes ef- 
fective. For the AIME directory and 
for the Personals column of the 
JourNAL oF PETROLEUM TECHNOLOGY, 
additional information is desired. The 
form below is provided for your con- 
venience, and should be sent to Petro- 
leum Branch, AIME, 800 Fidelity 
Union Bldg., Dallas 1, Tex. 


Name 
Membership No. 
Old Address 


New Address for Publications 


Address for Directory Listing and Job 


Title. 


List below your former title or com- 
pany position, nature of your new po- 
sition, or other information of interest 
to your associates for publication in the 
JOURNAL OF PETROLEUM TECHNOLOGY. 


One month normally required for 
change of address. 
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E. E. Elber (USN, retired), Director 
of Engineering for Bucknell Univer- 
sity. 

Foreign Service and Equipment, by 
A. W. McKinney, president of Na- 
tional Supply Co. 

The Future Growth and Financial 
Requirements of the World Petro- 
leum Industry, by Joseph E. Pogue, 
consultant for Chase Manhattan 
Bank. 

Long Range Supply and Demand, by 
Charles Hedlund, manager of Co- 
ordination and Economics Dept., 
Standard Oil Co. (N. J.). 

The program will probably con- 
sist of 10 papers when completed. 
Selected speakers are being requested 
to present papers on the subjects of: 
application of electric computing to 
the oil industry, the habitat of oil, 
the economics of well spacing, and 
sonic drilling. 


Host Section 


Development of the technical pro- 
gram for this meeting and arrange- 
ments for the speaker at the Petro- 
leum Branch Dinner are being han- 
dled by the New York Petroleum 
Section. Zeb Mayhew, Standard Oil 
Co. (N.J.) is program chairman, re- 
placing Henry E. McAuliffe, Socony- 
Mobil, who was recently transferred 
to Colombia. Dale Nix, Arabian 
American, is chairman of the New 
York Petroleum Section. 

This is the first year in which the 
New York Petroleum Section has 
had charge of the Petroleum activi- 
ties for the Annual Meeting. The 
formation of the New York Petro- 
leum Section in 1954 made it pos- 
sible and desirable for this section 
to develop the program. Prior to this 
time the Petroleum Branch Tech- 
nology Committee has handled all 
Annual Meeting programs. 

The Petroleum Branch Executive 
Committee is pleased with the pro- 
gram which is being developed. 
Papers selected to be presented on 
this program are of particular inter- 
est to people of the New York area; 
therefore, a larger attendance than 
in past years is expected. When held 
in New York on alternate years, this 
meeting will serve as somewhat of a 
regional meeting for the New York 
area, and as such it can be of great 
value to the New York Petroleum 
Section. 








BEING DEVELOPED for AIME ANNUAL MEETING 


Tentative Petroleum Program 
Saturday, Feb. 18 
10:00 a.m.—Meeting of Council « 
Section Delegates 
Sunday, Feb. 19 
2:00 pm.—Open Meeting o 
AIME Board of Dire 
tors 
Monday, Feb. 20 
9:00 a.m.—Open Meeting of P« 
troleum Branch Execu 
tive Committee 
Noon—A IME Welcomin 
Luncheon 
Afternoon—Petroleum Technic 
Session 
Tuesday, Feb. 21 
Morning—Petroleum Technica 
Session 
Afternoon—Annual Production Re 
view Session 
Evening—Petroleum Branc 
Dinner 
Wednesday, Feb. 22 
Morning—Petroleum Technica! 
Session 
Evening—AIME Annual Banquet 
and Dance 





Bob Cannon (left), vice-president 
of the University of Oklahoma of Pi 
Epsilon Tau, presents a bronze key 
plaque to William D. Owsley, vic« 
president of Halliburton Oil Well 
Cementing Co. and a member of th« 
Petroleum Branch Executive Com 
mittee. Through presentation of this 
plaque, Owsley was made an honor 
ary member of the national honorary 
petroleum engineer’s fraternity. A! 
this meeting he delivered a talk en 
titled “The Petroleum Engineer 
His Way of Life.” 


New AIME Director 
Named 


At the recent AIME Board meet 
ing, Charles R. Kuzell was named to 
succeed E. 1. Oliver, deceased, fo: 
the latter’s unexpired term as AIM! 
Director, until February, 195° 
Kuzell’s regular te:m as director fo 
three years begins at that time. 









In New Member Applications 





BRANCH SHOWS 43 PER CENT RISE 


VIRGINIA B, DAGGETT 


November was another record- 
breaking month for membership, 
with 193 applications received. The 
previous high for a single month set 
in October was 186. As we entered 
the last month in the 1955 member- 
ship contest, 1,111 applications had 
been received, up 43.4 per cent from 
775 at the same point last year. 

The outstanding increase of the 
Gulf Coast Local Section continued, 
boosted by membership committee 
activity at its two-day Formation 
Evaluation Symposium in late Octo- 
ber. Ninety-nine applications were 
credited to Gulf Coast in Novem- 
ber, giving the section a 34.7 per 


Proposed for Membership, 


California 

Bell—Marquis, Duane Elvin (M). 
Huntington Park — Subt, Frederick George 
(A). 

La Habra—Hildebrand, Alexander (M). 
Los Angelee—Marshall, Ira Allison, Jr. (M); 
Stanley, Charles Calvin (A); Travers, Wil- 
liam J. (R, M). 

Taft—Murphy, Denis Michael (J). 
Ventura—Larson, Robert Reynold (C/S-J-M). 
Colorado 

Denver—Gordon, Archie (M); Sneed, Robert 
W. (C/S-J-M). 

Illinois 

Mt. Vernon—Dost, Raymond Matthew (M). 
Kansas 

El Dorado—Halstead, Louis Norman (M). 
Wichita—Gary, John Thomas (C/S-J-M): 
Steves, Herbert Burdett (M). 

Louisiana 

Baton Rouge—Graves, Ernest Doyle, Jr. (J). 
Franklin—Baribault, Jules Donat (M). 
Lafayette—Baker, Joe B., Jr. (A); Haube, 
Edmund Dale (M); Thomas, Earl F. (M). 
New Iberia—Young, Ola Olympia (A). 
New Orleans—Favier, Pierre Pascal (M); 
Keating, Charles Crews (J); Norman, Ben- 
jamin Franklin, Jr. (M). 
Shreveport—Bynum, Winfred Wilson (M); 
Haynes, James Edward (M); Martin, Monty 
G. (M);: Thompson, Robert Earl (J); Wood. 
Jesse Winston (M). 
Venice—Dawsey, Albert 
Suggs, William Pat (J). 
Michigan 

Gra Rapids 
New Mexico 
Hobbe—Dennis, Jack Milton (J); Frick. 
Kenneth Dare (M); Luechi, Harold Joseph 
(J); Palmer, Aubrey King, Jr. (R, 
C/S8-S-J). 

New York 

Bedford Village—Davis, Dwight Ed (R, M) 
New York—Black, Thomas Oscar, Jr. (J); 
Volkman, William Edward (M); Wright, Al- 
bert (M). 
Pound Ridge 
Oklahoma 
Ardmore—Hill, John James (M). 
Drumright—Frazier, Ernest Richard (A). 
Lindsay—Fisher, Charles Albert (J). 
Oklahoma City—Lightner, Earl Miller (R, 
C/S-S-M). 

Purcell—Michie, Troy Willard (J). 
Tulsaa—Bridwell, Harold Cloyd (J); Gardner, 
Dan Lewis (R, J); Hunt, Elton B., Jr. 
(C/S-J-M); Moody, Wayne Charles (M); 
Newcombe, Jack (M): Snow, Robert Dale 
(M); Stoore, Philip, Jr. (A); Tillerson, Roy 
Bernard (M). 

Pennsylvania 

Bradford—Wood, Quentin Eugene (M). 
Texas 

Abilene—Brogdon, Bill R. (M); 
Owen (M). 

Amarillo—Tomlinson, Guy H., Jr. (M). 

Bay City—Kilpatrick, Robert Emmett (M). 
Baytown—Capshaw, Thomas David (J); 


louis, Jr. (J); 


Lazar, Albert (M). 


Coenen, Edward Joseph (M). 


Townley. 
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ASSISTANT SECRETARY 
PETROLEUM BRANCH, AIME 


cent increase for the year to date. It 
is interesting to note that the larg- 
est Petroleum Branch section now 
has the second highest per cent in- 
crease, exceeded only by two per- 
centage points by the Panhandle Sec- 
tion. Contest standings for all sec- 
tions are shown elsewhere on this 
page. 

All applications postmarked by 
midnight, December 31, will be eli- 
gible to count in this year’s member- 
ship contest. The winning section in 
each of the four groups will be 
awarded a slide projector, an AIME 
banner, or other suitable prize. 


Petroleum Branch, AIME 


Spriggs, Robert Paul (J). 
Beaumont—Bateman, Sammie J. (M); 
Eckles, Wesley Webster (J); Hayes, Edwin 
Garwood (M); Maddux, Henry Sherman, Jr. 
(J); Noone, Edward James (M); Pruett, 
Gayle Moreland (M). 
Beeville—Giese, Lester 
Withers, Bruce, Jr. (J). 
Bellaire—Dawson, Arch (R, M); Dieckman, 
John June (J); Fuqua, Marjorie (M); Gore, 
George Houston (J); Hottman, Clarence 
Edward (J); McCurtain, Edmund Elum 
(J); Simmons, Franklin Page (J); Smais- 
trla, Robert N. (J); Stelzer, Roland Brunn 
(J): Townley, John Paul (J). 
Boling—Paden, James Root (J). 
Dallas—Clote, Paul J. (J); Melsheimer, Al- 
bert Henry (R, M); Reeves, John Lloyd 
(J); Shaw, Roland Clark (A); Sinex, James 
Mitchell (M); Slusser, Marion Liles (M); 
Tompson, Gilbert Carr (C/S-J-M). 
Fort Worth—Britton, Stanley Grant (M); 
Powers, Richard Ear! (M); Robertson, John 
James (A); White, James Robert (M):; 
Yager, Charles Ervin (M). 
Galena Park—Ramey, Darol Keith (J); 
Shannon, William L. (M). 
Genoa—Davies, William James (R, C/S-S-J). 
Houston—Akkerman, William Herman (M); 
Ballard, Junius Chambers, Jr. (J); Banks, 
Overton Barnes (M); Barrett, Joe Bill (J): 
Beard, James Ijams (J); Blumberg, Randolph 
(J); Boyd, Walter Andrew, Jr. (J); Brixey, 
Austin Day, Jr. (M); Broussard, Douglas 
Eugene (M); Chisholm, James P. (J); 
Cobb, Lewis Allen (J); Colby, John Byron 
(A); Connelly, Francis Becker (J); Craig, 
J., Jr. (J); Cralle, George Leroy 
; Crain, Marvin G. (J); Derrick, A. M., 
(C/S-J-M); Dewan, John Tierney (M); 
Diener, Milton Lawrence (M); Doyle, James 
Joe (M); Drummond, William Temple (J); 
Ekey, John Walter (J); Evans, John Billy 
(R, M); French, William Lee (J); Gimblet, 
Ernest Carlyle (M); Glasgow, Herschel Mil- 
ler (J); Grimes, Erwin Elwood (J); Gris- 
sett, John Neal (J); Haddox, R. Putnam 
(C/S-J-M); Heiny, Rems Chandler (J); 
Henderson, Paul William (J); Hinerman, 
David Arthur (J); Hodgson, Glen Reid 
(M); Humble, James M. (J); Kellough, 
John Stevens (M); Koch, William Martin 
(M); La Croix, Sam Henry (J); Lawrence, 
William Jackson (M); Malek-Aslani, Morad 
(M); Manson, Donald Humphrey (M);: Mar- 
shall, Ray W. (A); Martin, Richard Wendell 
(J); MeKeever, Ira Everett, Jr. (J); Mur- 
phy, Arthur Gilman (J); Nelson, Paul Mor- 
wood (J); Nelson, William Paul (M); 
Oliver, Raymond Alexander (J); Phillips, 
William Lawrence (A); Pianta, James 
Eugene (J); Raabe, William Louis (J); 
Rehm, William Arthur (J); Reiter, William 
Arthur, Jr. (M); Rhodes, Roy Marvin (J): 
Riep, Otto (M); Robb, D. Douglas, Jr. (M); 
Robertson, Wayne W. (J); 
James (M); Stovall, George Alston (J); 


Edgar, Jr. (J): 
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PETROLEUM BRANCH MEMBERSHIP 
CONTEST 
Eleven-Months Report 
Covering Applications Received 
Jan. 1-Nov. 30, 1955 


Jan. 1 Applico 
Member tions 
Sectior ship Received 


Per Cent 
Increase 


Group A 
] PANHANDLE OU 36 
2. Hobbs 8& 34 
3. Mississipr 5 20 
4. Billings 
Petroleun 6 19 
Hugoton 2 18 
6 North Texa 15 17 
7. South Pia 2 
Group B 
] KANSAS 41 y 20 
2. Denver 
Petroleum ! g 19 
East Texas 1 19 
3 Lou-Ark 16 
4. West Central 
Texas 135 2 8 
5. Wyoming 8 
5. Illinois Ba 
Petroleum } 5.§ 
Group C 
1 DELTA 30 18 
Venezuela 
Petroleur i. 16 
3. Fort Worth 181 14 
4. Oklahoma City 14 
5. Permian Basin 11. 
6 New York 
Petroleum 212 ! oy 
7. San Joaquin 
Valley 169 7 
Group D 
1. GULF COAST 929 
2. Mid-Continent 553 
3. Dallas 433 
a Southwest Texas 461 
5. Southern Califor 
nia Petroleum 597 


Total Applications 
from Sections 
Miscellaneous 44 
Foreign 47 
Total Applications Re- 
ceived Jan. 1-Nov.30 1111 











Starr, George 


Thomas, Bill Eugene (J); Thomas, Sherwood 
Lindsay (J); Trimble, Richard L. (M): 
Trowbridge, Charles Eugene (J): Vercellino, 
Wilton Canova (M); Vestermark, Lyle Ar- 
mand (M); Wherry, R. Britton (J). 
Katy—Burley, James D. (J); Corbell, Billy 
Joe (J); Gorden, Cecil Dennis (M); Hilton, 
Byron Lee (J); Thomas A. D. (M). 
Liberty—Daniels, Whitfield Floyd (M). 
Midland Davis, Donald William (J); Fitz- 
gerald, Thomas Allen (J); Fitzgerald, Wil- 
liam Patton J) Holley, Joseph Wilson 
(M) Jones, Samuel Clayton (M); Rock- 
wood, Spencer Henry (R, M); Vickery, Jack 
Hill (M Whitfield, Marion Anderson (A); 
Woodward, Maurice E. (M); Yarbrough 
Richard William (J). 
Monroe City—Baugh, Edgar 
McCune, John Stanley (M); 
Paul (J). 

Odessa—Langford, Billy J. (J). 
Pampa—Michae!, Dan Gauntt (J); 
Irwin Thomas (J) 

San Antonio—Mitchell, John K. (A). 
Thompsons—Grinstead, Fred Edward (J) 
Tomball—Peterson, Gerald Lee (J). 
Tyler—Williams, James Douglas (M). 
Wharton—Eastman, Frank Allan (M); Mce- 
Culloch, Robert Andrew (M). 

White Oak—Haynie, Moss (A). 

Wichita Falls—Graf, Arthur F. (C/S-J-M) 
Martin, Joseph Frederick (A); O'Neal, J. A 
(A). 


Gratton (J); 
Pipes, Kenneth 


Winter 


Winnsboro—Wofford, Chancey Jackson (J) 
British West Indies 
Trinidad—Howe, Sta 
Canada 
Calgary—Buckles, Re 
Regina Foste Gex 
M) 

England 

Birminghan Coop Alan Charles (J) 


bert Sheldon (J). 
rge Michael Hendersor 


Germany 
Hamburg 
Persian Gulf 
Qatar—O'Connell, Patrick D. A. (M) 
Venezuela 
inaco—Bates, William Clifford (J) 
Rhoads, Daniel Wayne (M) 
Cobb, Dow C. (A) 


Hehlert, Otto Ernst (M) 


Barcelona 
Caracas 
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MEET THE AUTHORS 
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chemistry. Following graduation he 
worked as research chemist on well 
completion for Oilwell Rescarch, 
Inc. Joining Union Oil Co. in 1951, 
he is engaged in well completion 
work with the production research 
group of the Research Dept. 





R. F. KRUEGER holds a BS degree 
in physics from Ripon College, 
Wisconsin, and a MA from the Uni- 
versity of Illinois. Since 1948 he 
has worked for Union Oil Co. on 
secondary recovery and well com- 
pletion problems and is_ section 
leader of the Well Completion and 
Field Service Group. 





R. L. NEWACHECK is a senior 
mechanical engineer with Tracerlab, 
Inc., Western Div., joining the firm 
in 1950. He holds a BS degree in 
mechanical engineering from Kansas 
State College. After graduation he 
joined the Pacific Gas & Electric 
Co., and is a former instructor in 
blueprint reading. 





E. E. ANDERSON received a BS 
degree in chemistry from San Diego 
State College in 1950. He joined 
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(Continued from Page 26) 


Tracerlab, Inc., Western Div., as a 
radiochemist following work with 
Consolidated Aircraft Corp. He is 
presently carrying out development 
of a number of commercial radio- 
isotope preparations and industrial 
applications. 





et 


JEROME KOHL was graduated 
from California Institute of Tech- 
nology in 1940 with a major in ap- 
plied chemistry. Prior to joining 
Tracerlab, Inc., in 1948, he was em- 
ployed by Tide Water Associated 
Oil Co. He is now chief engineer for 
Tracelab’s Western Div. 





MOHAMED Mortapa received a BS 
degree in chemical engineering from 
the University of Cairo in 1946 and 
MS and PhD degrees in petroleum 
engineering from the University of 
California, Berkeley. Joining Mag- 
nolia Petroleum Co. in 1954, he is 
a senior research engineer for Mag- 
nolia’s Production Div. 





J. W. McGuezE is employed in the 
secondary recovery section of the 
Cities Service Research and De- 
velopment Co. He joined Cities 
Service after graduating from Okla- 


homa A&M College with a BS « 
gree in chemistry in 1951. 





Jack NEWCOMBE joined the Cities 
Service Research and Development 
Co. in 1947 after graduation fron 
Oklahoma A&M College with a: 
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INVESTIGATION of VARIOUS REFINED OILS for FORMATION FRACTURING 


ABSTRACT 


The properties of some of the 
common fracturing oils, such as 
sand-falling rate, fluid loss and 
viscosity, are a function of the thre« 
major constituents of the oil. These 
are: asphaltic material, paraffin 
waxes and a base oil. The base oil 
can be either high or low viscosity. 

Asphaltic material increases vis- 
cosity by thickening the oil, decreas- 
ing the rate of fall of sand through 
the oil, at temperatures up to 
200° F. 

When large amounts of paraffin 
wax are present, much of the wax 
is not dissolved, but is present in 
the form of microscopic crystals. 
This, of course is limited to tem- 
peratures below the melting points 
of waxes. As a result, paraffin waxes 
do not function primarily as thick- 
eners, but rather as bulking ma- 
terials. 

The function of wax, 
asphaltic material, in a low viscosity 
base oil is uniquely unlike that of 
wax with asphalt coexistent in the 
oil. The resulting difference in wax 
crystal size governs the sand-falling 
rate to a large extent. 


without 


An ideal tyve of fracturing oil is 
one containing neither asphalt nor 
waxes. In such cases, the viscosity 
and resulting low sand-falling rate 
are apparently due to soluble resin- 


Revised manuscript received on Oct. 24, 
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ous hydrocarbons, which are com- 
patible with formation oils. 


INTRODUCTION 


Since the introduction of hydrau- 
lic fracturing as a method of stimu- 
lating production from oil and gas 
wells, a variety of fluids have been 
employed to open the fractures and 
carry the propping agent. Such fluids 
include gels, emulsions, and oils. In 
the latter category, the use of refined 
oils has increased tremendously, and 
it is probably safe to assume that 
such oils represent more than one- 


half of the total fluid used in fractur-. 


ing. 

The general acceptance of refined 
oils for fracturing can be accounted 
for on the basis of their relatively 
low cost and general applicability. 
To minimize transportation costs, 
such refined oils are generally ob- 
tained within the local areas of oper- 
ation. As a result there are several 
kinds of oil used in different areas, 
with a wide variation in the basic 
composition. These types may vary 
from a fair grade of lubricating stock 
to a highly viscous cut-back asphalt, 
or road oil. 

Originally, little attention was paid 
to the composition of the oils used 
for fracturing, and any oil which 
would carry sand was considered 
satisfactory. Recently, however, serv- 
ice companies and well operators 
have become more concerned with 
the type of oils used for fracturing. 


DOWELL INCORPORATED 


With this in mind, the present work 
was undertaken to evaluate and com 
pare many widely used fracturing oils 
to determine which constituents 
were most desirable for fracturing 
purposes. 

Some of the important properties 
of fracturing fluids are sand-falling 
rate, viscosity, and fluid loss. Other 
properties, that are actually of no less 
importance, include flash point, pour 
point, emulsifying tendencies with 
formation brines, and compatability 
with the formation oils. This discus 
sion has been limited, however, to the 
first mentioned properties in order to 
give some idea of the general com 
position of the refined oils, how they 
function, and why. 


PROPERTIES OF FRACTURING 
OILS 


The properties of fracturing oils 
are a function of the three majo: 
constituents of the oil. These may be 
designated as: (1) asphaltic material, 
(2) paraffin waxes and (3) a base oi! 
which may be of either high or low 
viscosity. The asphaltic material clas 


TABLE 1 — ANALYSIS OF TYPICAL FRACTURING 
OIL WITH REGARD TO ASPHALT AND WAX 


CONTENT 
% Asphalitene” 
% Total (in asphalt %, Base 
Oil Asphalt! fraction) _% Wax? Oil 
a 33 5 3 64 
B 30 2.2 21 49 
Cc 5 trace trace 95 
Dd 17 3.1 10 73 
E 20 1.1 33 47 


1Distillation Residue ot 360°C (10 mm pressure 
Precipitation with 10:1 Pentane (85°F) 
*'Dewaxed with Acetone dilution 7:1 (32°F) 









sification, for the purpuse of this 
paper, includes any high molecular 
weight naphthenic constituents in the 
oil. The variation in the composition 
of different fracturing oils is indi- 
cated by the data shown in Table 1. 

Low as well as high percentages 
of both asphalts and waxes are found 
in various oils. In Table 1, the sec- 
ond column indicates the total as- 
phaltic material present in the oil, 
while the third column shows the 
percentage of asphaltene materials 
present in this fraction; the fourth 
column indicates the high molecular 
weight paraffinic material present; 
and the percentage of base oil pres- 
ent is shown in the fifth column. 
(This base oil consists of low molec- 
ular weight naphthenic materials and 
paraffins. ) 

The effect of the various constit- 
uents on the viscosity of the oil is 
illustrated in Table 2. In order to 
keep the phenomenon of gel strengih 
from interfering with measuring ac- 
tual viscosities, the latter were deter- 
mined by an Ultra-Viscoson (Model 
XV, Rich-Roth Laboratories). Vis- 
cosities in the column headed Dis- 
tillate show the effect on viscosity 
of removing the asphaltic material. 
Figures in the Base Oil column rep- 
resent the viscosities of the base oils, 
after both asphalts and wax have 
been removed. It can be seen from 
Tables 1 and 2 that Oil A achieves 
its viscosity almost entirely from as- 
phaltic material. OU C, by contrast, 
achieves its viscosity mainly froma 
the use of a highly viscous base oil. 

The viscosity of most fracturing 
oils with various fractions omitted 
demonstrates that an oil of low vis- 
cosity is usually used as a base, when 
appreciable amounts of either asphalt 
or wax are used in compounding a 
fracturing oil. With asphalts and 
waxes present in only small amounts, 
a low sand-falling rate is achieved 
by using a high viscosity base oil, 
similar to the heavy weight lubricat- 
ing oils. Sand-falling rate is usually 
measured by observing the rate at 
which graded sand, 20-40 mesh, falls 
through the oil. 


This is generally considered to be 
related to the viscosity of the fluid. 
A low rate is of importance for a 
fracturing oil since the sand must 
remain suspended in the oil during 


TABLE 2 — VISCOSITY* (CP @ 80°F) OF FRAC- 
TURING OILS 
Base Oil 

rr" Original Disti’'‘ate (asphalt and wax 
oil _ Oil (asphalt removed) removed) 

A 215 13 

8 184 20 

Cc 215 185 

ie) 200 35 

E 124 60 17 
*Viscosity measured on Ultra-Viscoson 


40 


the treatment. A rate of 5 ft/minute 
or less has been found satisfactory 
for current pump capacities. The 
sand-falling rate in water, for com- 
parative purposes, is approximately 
27 ft/minute. As pointed out later, 
it has been observed that there is not 
necessarily a direct relationship be- 
tween viscosity and sand-falling rate. 


ASPHALTIC MATERIAL 

Asphalts or asphaltic material are 
usually present in the form of micel- 
les, or aggregates of macro-molecules 
within the oil base. The asphaltic 
material reduces the sand-falling rate 
by two methods. In addition to in- 
creasing the viscosity of the base oil, 
the micellular aggregates of asphalt 
further decrease the sand-falling rate 
by a bulking action. It appears prob- 
able that this bulking action results 
from the physical dimensions of the 
micelles. 

Asphaltenes are oxidized high mo- 
lecular weight polycyclic aromatic 
hydrocarbons. They are not actually 
soluble in oil. According to Sachanen’, 
the insoluble asphaltenes readily ad- 
sorb resins and aromatic hydrocar- 
bons and are dispersed in the oil. 
Thus, the resins and aromatic hydro- 
carbons colloidally protect the as- 
phaltenes, as all three materials exist 
as miscelles in equilibrium with the 
hydrocarbon constituents of the oil. 
The resistance offered by these 
micelles serves to slow the fall of the 
sand particles. 

The asphaltic miscelle effect is rep- 
resented graphically in Fig. 1, show- 
ing the viscosity-sand-falling rate re- 
lationship. The curve of an oil of 
Newtonian properties, containing no 
asphalt, follows a nearly straight line 
function. The effect of asphalt is to 
impart lower sand-falling rates to a 
base oil of comparable measured vis- 
cosity. For a specific viscosity, the 
sand-falling rate of the blended oil 
is lower than that of the base oil (a 
true solution). An oil with a low 
sand-falling rate and low viscosity 
might be of advantage since friction 
pressures encountered during a frac- 
ture treatment using such an oil 
would be lower than when using a 
base oil with equal sand-falling rate. 

The effects of asphaltic material 
are not all favorable, however. Re- 
fineries and pipeline companies are 
becoming more concerned with the 
dilution or contamination of crude 
oils with “foreign fluids” such as 
fracturing oils. This concern for com- 
patability, in many cases, stems from 
the use of asphaltic fracturing oils. 
The colloidally dispersed asphaltenes 


! Reference given at end of paper. 
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MICELLE EFFECT OF 
ASPHALT IN OIL. 


may be upon dilution 
with This is especially 
true for high gravity crudes. When 
the fracturing oil is diluted, the resins 
and aromatic hydrocarbons tend to 
dissolve, the micelle equilibrium is 
destroyed, and the asphaltenes coag- 
ulate and precipitate. The precipita- 
tion of a crude asphalt can be ex- 
tremely troublesome, so far as plug- 
ging and fouling tubing and lines is 
concerned. 


precipitated 


crude oils. 


Some refinery processing includes 
the use of propane dilution for as- 
phalt removal. This same precipita- 
result of dilution, may 
occur during a fracture treatment. 
A highly asphaltic fracturing oil 
could, therefore, plug the flow chan- 
nels of a formation and offset any 
advantage gained as a result of the 


tion as a 


fracturing treatment. 


PARAFFIN WAXES 


Whereas asphaltic materials are 
usually present in the form of micel- 
lular aggregates, paraffin waxes, pres- 
ent in the oil, exist either in solution 
or in the form of wax crystals. Such 
wax crystals in the refined oil used 
for fracturing treatments, display a 
bulking effect which, in certain re- 
spects is similar to asphalt. The crys- 
tals of wax delay the fall of the sand 
particles without increasing the ac- 
tual viscosity of the base oil. Fig. 2 
demenstrates and compares this ef- 
fect. The curve representing fractur- 
ing oil with paraffin wax is more 
divergent than that for the asphaltic 
oil, since the waxes have a definite 
melting or crystallization range. 
Above the melting point, the wax 
effect does not function; therefore, 
the curves coincide. To further illus- 
trate this point Fig. 3 shows the ef- 
fect on sand-falling rate, at various 
temperatures, of the various com- 
ponents of a typical fracturing oil 
(Oil E) 
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Fic. 2 —- CRYSTAL EFFECT OF WAX 
IN OIL. 


The size of the wax crystals has a 
unique role in the characteristic 
property of sand-falling rate control. 
The presence of asphalt with the wax 
in the oil results in a smaller particle 
size than that obtained without the 
asphalt present. Apparently the wax 
particles become coated with the res- 
inous asphalt, and are thus restricted 
from material growth. Fig. 4 shows 
three photo micrographs: one of an 
asphaltic, wax-base, fracturing oil: 
and the second of the base oil frac- 
tion, with the asphaltic material re- 
moved. By omitting the asphalt, the 
wax crystals were 2 to 3 times as 
large. This increase in crystal size 
actually reduced the sand-falling rate, 
although the viscosity was decreased 
by the omissions of the asphalt por- 
tion. As the third photo micrograph 
shows, addition of the asphalt, with 
heating and recooling, caused a re- 
duction in crystal size of the wax 
particles, increasing both the viscos- 
ity and the sand-falling rate. 


A refined oil that is high in par- 
affin waxes could be a source of 
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troubl: if used as a fracturing fluid 
for low temperature wells. In wells 
with temperatures below the melting 
points of the waxes, it is possible 
that wax deposits might plug the flow 
channels. Wax deposits also might 
plug tubing and lines, or accumulate 
in tank batteries. 

Wax particles, as well as as- 
phaltenes, are hydrophobic in nature. 
Because of this, these materials have 
been observed to stabilize water-in- 
oil type emulsions. This is not of too 
much importance in most fracturing 
treatments; however, it should be 
considered since, in isolated instances, 
emulsion troubles do plague recovery 
operations. The absence of both wax 
and asphaltic material in the fractur- 
ing oil should prevent aggravating 
any potential emulsion problem. 


FLUID LOSS 

Another important property in a 
fracturing oil is fluid loss. If the loss 
of the fracturing oil into the forma- 
tion is too great, sufficient pressure 
build-up to open a fracture is dif- 
ficult to obtain. If most of the oil 
was lost into the formation while 
the fracture is created, the sand 
would have no carrying media. 

Low fluid loss in a fracturing oil 
may be accomplished in two ways: 
(1) by incorporating a plugging ma- 
terial into the oil, and (2) by in- 
creasing the viscosity of the oil. As- 
phaltic micelles, which are actually 
insoluble, as well as wax crystals, 
show a tendency to plug formation 
pores, reducing fluid loss. As a rule, 
fluid loss reduction by this type of 
plugging is not desirable, as the 
plugging materials may be difficult 
to remove from a formation. This, 
of course, could be detrimental, in a 
fracturing treatment. The reduction 
of fluid loss by increasing the vis- 
cosity of the fracturing oil, on the 








other hand, avoids plugging of fo: 
mation pores. 

Fluid loss properties of a fractu: 
ing oil are, in part, a function of 
relatively insoluble asphalts as shown 
by the data in Table 3. It would 
appear that the asphalts are not com 
pletely dissolved and therefore tend 
to plug, thereby reducing the loss o! 
fluid. This plugging effect, that oc 
curs at the face of a formation dur 
ing a fracturing treatment, would 
tend to minimize contamination with 
in the formation itself. The amount 
of insolubles governs, to a large ex 
tent, the penetration of the fluid into 
the formation flow channels. The as 
phalts could then precipitate, result 
ing in a decrease in potential pro 
ductivity. The extent of precipita 
tion will depend on the nature of 
the native oil in the formation 
Formation crude containing large 
amounts of light ends will precipitate 
more asphaltic material from a frac- 
turing oil than a crude with smal! 
amounts of light fractions. 


CONCLUSIONS 


In summation, the three majo: 
constituents of refined oils used as 
fracturing fluids, may function in 
several ways. The asphalts, waxes, 
and base oil individually and collec 
tively affect the viscosity, sand-fall- 
ing rate, fluid loss, and compatibility 
of the oil with produced crude oils. 
Asphaitenes, in the asphaltic portion, 
and wax crystals impart low sand- 
falling rates at lower viscosities than 
would be required for a base oil, 
alone. A significant disadvantage of 
the asphalt and wax is their tendency 
to plug minute flow channels, and to 
contaminate produced crude by pre 
cipitation of insolubles, when diluted 
in the formation. To induce low 
fluid loss by the use of plugging 
materials is not necessarily harmful, 
if the material used to accomplish 
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TABLE 3—EFFECT OF CONSTITUENTS OF A FRAC- 
TURING OIL ON FLUID LOSS 


Composition: 15% Asphalt, 20% Wax, 65% Base 

Oil (20 cp at 80° F measured with Ultra-Viscoson) 
Fluid Loss (API Code 29) 
@ 80° F 

75 m1 /30 min 


Test Fraction of Oil 
@ 200° F 
1 Original oil 100 m1/30 min 
2 with asphalt 
removed 
3 with wax 
removed 
4 wifwax & as- 
pholt re- 
moved 400 m1/4 min 
5* High Viscosity 
base oil 400 m1/13 min 
(no asphoit, 
no wax) 


*Viscosity of 140 cp at 80° F (Ultra-Viscoson) 


70 m1/30 min 400 m1/2 min 


75 m1/30 min 90 m!/30 min 


400 m1/2 min 


400 m1/5 min 


this is easily removed by dissolution, 
decomposition, or other assured 
means. Removal is usually difficult 
to accomplish in the case of asphalts 
and waxes. 

The most satisfactory refined oils 
for use in fracturing should be ones 
without appreciable amounts of 
waxes, and especially asphaltic ma- 
terial. A more desirable fracturing 
oil would be one whose viscosity, 
sand-falling rate, and possibly fluid 
loss are a function of the base oil 
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itself. Using this type of oil should 
result in more complete compatibil- 
ity with the crude oil. In addition, 
a minimum of foreign deposits would 
be introduced to restrict flow through 
the newly fractured formation, and 
maximum benefit could be realized 
from the treatment. 
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SOME USEFUL TABLES for APPROXIMATING SMOOTH CURVES 


by FIFTH-AND-LOWER DEGREE POLYNOMIALS 
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W. P. SCHULTZ* 
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INTRODUCTION 


The use of computing machines to solve physical 
problems has made it imperative to represent physical 
data in a form computing machines can use. Although 
curve-fitting is an old and well-practiced art, there still 
seems to be some interest in questions regarding what 
methods can be used and why. Even though there are 
many good references on this subject, it may be of some 
value to discuss a particular method which has been 
found to be useful, and to make available the associated 
tables for its use in the hope that there can be saved 
a few of the tedious hours necessary to reduce many of 
the suggested methods to practice. Although a very brief 
discussion of the mathematical background of the tables 
is included immediately below, those interested only in 
their use can skip directly to the section on Use of the 
Tables. 


DIRECT FIT OF A SET OF POINTS 


A frequently encountered problem is: Given data 
which can be represented by a smooth curve (x), to 
find an easy-to-evaluate function F(x), with only a few 
coefficients, which is a good approximation of $(x) on 
the interval 0 < x <. x,. Consider F(x) as a linear com- 
bination of m + 1 arbitrary, single-valued, continuous, 
linearly independent functions, f,(x), 

F(x) = aff(x) + afd{x) + ... + aifi(x) 
Ae a Bs ae js amet oe 
If the wnbes of (x) is known or can ie read from the 
curve at the points x,;: (x.), (x,), o(x;) . 
o(x,), and if nm = m, then all the coefficients of F(x) 
can be determined, and F(x;) — $(x;) = 0 for j = 0, 1 
. .n. This follows from the solution of the set of linear 
simultaneous equations. 


$(x.) = afo(x.) + af,(x.) Y Gelate) 
o(x,) = af. (x,) D3 af, (x,) a ee EH 
$(x,) = Fitawes (x,) + er ) 4 a+ Anbm(X) 

(2) 


for the coefiicients a. 


Although the set a, from the solutions of the set of 
Equations 2 yields an F(x) which gives an exact value 
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of }(x;) for the n + 1 points, there is no restraint on 
F(x) in the interval x, << x < x,+,. Indeed, if the points 
(x,) are not sufficiently smooth, very serious differ- 
ences may arise between F(x) in the interval and the 
smooth curve $(x). 


LEAST-SQUARES FIT 


It has been suggested’ that if n is made to exceed m, 
and the set a, is found to minimize the function 


j=n 


a,,) = } [F(x,) d(x), 
j=0 


the function F(x) is much less sensitive to small errors 
in #{x,), and usually provides a much better fit to the 
smooth function ¢(x). This is known as the best fit in 
the sense of the least-squares error. 

To minimize y, set 


WB < s 


ov 
me = 0,i= 0,1, Sa eee. OD 
ad, 
The m + 1 equations from Equation 3 are 
k=n k=n 
a, ; f(x) fix) + a, > 
k=0 k=0 
k=n k=n 
h(x) Fi) eee ee fm (Xx) fi(x%.) = 
k=0 k=0 
hl Xy): Sak Re esa : i Goto.” Lan) 
one equation being obtained for each i=0,1,...m 


Inasmuch as the arbitrary functions f,;(x) may be evalu- 
ated for x,, kK = 0, 1, . n, and the ¢(x,) are known 
values of the curve to be fitted, the simultaneous Equa- 
tions 4 may be solved for the set a,. 


PARTICULAR FUNCTIONS 


Given the functions /,;(x) and the set #(x,), the solu- 
tions for a; require evaluating a number of sums of 
products a as well 2 as 5 the solution of anim + 1} order set 


'Reference given at end of paper. 










of simultaneous equations. Much of the work can be 
done once and for all if there is chosen a standard set 
of f(x) for the approximating functions and a standard 
set of (x,) where ¢ is to be approximated. Then the 
Equation 4 can be solved for the set a, to give the 
solutions 


j=n 


a,= ; b,; o(x;), i = 0,1,...m, 


j=0 
requiring at most m + 1 sums of n + 1 products to 
evaluate all the coefficients in F(x) for a given set 


$(x;). 


USE OF THE TABLES 


This description is inciuded to assist readers who have 
no taste for the foregoing discussion but who, nonethe- 
less, should like to fit graphical data with polynomials of 
degree five, or lowe 

Given a smooth curve ) 
interval on the must be selected for which the 
curve y(x) is to be fitted. By translation of the origin 
and change of scale, the interval along the x-axis must 
be made to be 0 < x < |. This interval must be divided 
into 10 sub-intervals by 11 equally spaced points, x = 0, 
0.1, 0.2 1.0, and the 11 corresponding values of y 
must be read off the curve y(x). Now the coefficients of 


as some function of x, an 


¥-aXIS 


a polynomial to approximate ¢ in the interval can be 
obtained using the appropriate table. 

Consider finding the coefficients of the third degree 
polynomial 

F(x) Gt Gn 
by use of Table 3 and the 11 points from a curve y(x). 
The value of a, is found as the sum of the products of 
the 11 values on the top (x’) line of Table 3 with the 
11 values of y(x) taken from the curve at the points of 
x corresponding to the values in the heading of the 
tables. The value of a, is found similarly by using the 
second (x') line of Table 3. The polynomial with the 
four coefficients evaluated by using Table 3 in this way 
is the unique third degree polynomial which minimizes 
the sum of the squares of the 11 differences between the 
value of the polynomial at the points x hd, Ose, 
1.0 and the y-values used in forming the sums of 
products. 

It might be well in conclusion to point out that many 
smooth curves y(x) are not approximated well by poly- 
nomials. It is always advantageous to evaluate the poly- 
nomial at the 11 input points to check the closeness of 
fit. In the event y(x) is not approximated well! by a poly- 
nomial, log y, y", or some other easily evaluated func- 
tion of y may yield a better fit. Judgment and experience 
will help in selecting the best function of y to use. 


The validity of Equation 5 may be seen most easily 
by the use of matrix notation. Equation 4 may be writ- 
ten as 


[D]a, [G]d(x;) (6) 


aX A,X 


where the element of [D] is d,,, = f(x) f(x), the 


k=o 
waa 


element of [G] is g, , f,(x,), @ is the column vector 


of m + 1 coefficients, and ¢(x,) is the column vector 
of m + 1 values from the curve (x). Premultiplying 
Equation 6 by [D]” gives 


a, = [D]"* [G]¢(4,) = [B]¢(s;), (7) 
which is equivalent to Equation 5 where the coefficients 
b,, are the elements 5, , of the matrix [B] = [D]' [G] 
of Equation 7. 


NUMERICAL VALUES FOR [B] 


The functions f,;(«) were chosen as f(x) x' and the 
interval 0 < x < 1 was divided into 10 equal parts by 
the set x, j/10, j7 = 0, 1, ... 10. Tables 1 through 5 
present the numerical values of the elements of [B] for 
m = 1, 2, 3, 4, and 5, respectively. An estimate of the 
precision of the individual elements has not been made, 
but tests with simple functions, e.g., y = x, y = c, yield 
approximations which are better than one part in 10,000. 
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TABLES 1-5 — NUMERICAL VALUES OF ELEMENTS OF [B] FOR m 1, 
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-0. 13636364 
45454545 


-0.09090909 
36363636 0. 


-13636364 -04545455 
-09090909 0.09090909 -1618181 -272 0 


- 18181818 
- 18181618 


- 31818182 
~ 45454545 


«27272727 
- 36363636 


+22727273 
+27272727 


-01398602 
- 33566434 
- 69930070 


-12587413 
- 29370630 
- 74825175 


-1048951¢ 
-88111688 
.-69930070 


-11188811 
- 13980014 
-04895105 


-08391608 
- 16550117 
- 16550117 


-02097902 
+ 95804196 
-04895105 


-07692508 
+51748252 
+ 69930070 


- 58041958 0. 37762238 
-20279720 -1.06293706 
+ 74825175 0.69930070 


+ 20979021 
-15617715 
-11655012 


Ill 


055944 
-002331 
447552 
- 165501 


-083916 
-039627 
-99.3007 
-827506 


-055944 
-674437 
-002331 
-467754 


-013986 
- 415695 
-030303 
.719503 


055944 
+288267 
- 526806 
-273504 


-083916 
-073058 
400932 
.407754 


- 118881 
.513598 - 165501 
- 128205 - 165501 
-719503 0 


-083916 


Iv 


-069930 
- 368.307 
-405598 
- 799532 
7.482518 


-041956 
- 330994 
-860128 
137522 
- 482518 


+209791 
-974339 
902078 
130538 
482518 


- 069930 
658913 
. 379969 
~238541 
482516 


-041956 
-205116 
20. 687627 
965031 
7.482518 


-069930 
~3.329441 
599063 
-029529 
-655012 


+104896 
-801495 
-043142 
«295260 
913752 


-034965 
400146 
-440552 
-590523 
655012 


v 


015730 
191151 


045446 

-5.675974 2. 
- 109307 -22.397372 
12417 75.824619 
. 88802 ~102.70301 
148218 48.074109 


0.973771 
~14, 426829 
69. 000090 


-006987 
«700070 
077492 
~ 91956 
-59201 
-098812 


094415 

-018550 
416930 
05414 
85294 
148218 


-041952 
- 3.205087 
20.687646 
-53948 - 965055 
-90190 7.482476 
096812 0 


-089159 
-166224 
. 799092 
74.225850 
544256 
-024703 


027978 
-228024 
+635413 
65370 
33315 
-098612 


041952 
-429619 
-117073 
~145.75444 
137.66779 
~48 .074109 
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THE CONDUCTION of HEAT INCIDENT to the FLOW 
of VAPORIZING FLUIDS in POROUS MEDIA 


FRANK G. MILLER 
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RALPH A. SEBAN 


INTRODUCTION AND PURPOSE 


Problems relating to thermal methods of oil recov- 
ery have been given increasing attention during the past 
year. The nature of the physical and chemical processes 
underlying thermal recovery are not yet well under- 
stood. The need for pertinent basic information is recog- 
nized generally. One of the processes involved is the 
transfer of heat by conduction in oil reservoir rocks 
containing moving reservoir fluids. This technical note 
deals with heat conduction in a laboratory sand column 
filled with a moving hydrocarbon and therefore should 
be contributory to the presently inadequate fund of 
knowledge on the over-all subject. 

In an analysis of the flow of vaporizing propane 
through a horizontal column of sand, one of the pres- 
ent authors’ demonstrated correspondence between the 
measured values of pressure along the length of the 
column and those predicted by a theory which post- 
ulated that the enthalpy of the propane was invariable 
along the length of the column. This postulate was 
implied as a deduction from the principle of conserva- 
tion of energy but the details of the deduction were 
not given. The specific purpose of this technical note 
is the presentation of these details both as an adjunct 
to the earlier analysis and to supply some additional 
details on the processes involved in this type of flow. 
Results disclose that rates of heat transfer by conduc- 
tion, in the direction of fluid motion, generally would 
be smail for systems of the type studied. 


FUNDAMENTAL ENERGY EQUATION 


The principle of conservation of energy may be ap- 
plied to a section of the porous material on the basis 


“References given at end of paper. 
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that heat transfer, like fluid flow, is possible in the « 
direction only. Transfer of heat in the y and z direc- 
tions is prevented by suitable insulation or by taking 
the medium as infinite in directions perpendicular to 
the x axis. For steady state heat conduction and condi- 
tions of steady fluid flow the principle of conservation 
of energy specifies that any reduction in the enthalpy 
rate increases the flow of heat in the same proportion. 
Stated mathematically, 
dH KA dT 


—_- = — : re eo es ee ee (1) 
dx W dx’ 


KA dT. : : a 
where — W dx’ the flow of heat in the positive x 


direction. 


TEMPERATURE DISTRIBUTION IN DIRECTION 
OF FLOW 


The column of sand used in the previously described 
research is divisible into three different flow regimes, 
each of which may be analyzed independently. In the 
entrance region the liquid approaches the saturation 
state. A transition region is postulated following the 
transverse section at which the saturation state of the 
liquid is attained. In this transition region gas and liquid 
phases are both envisioned as present, with no vaporiza- 
tion occurring and with fluid movement restricted to 
the liquid. The upstream end of the exit region is at 
the transverse section of the column where the frac- 
tion of the pore space occupied by gas first becomes 
great enough for gas flow and where vaporization of 
the liquid begins. Fig. 1 indicates schematically the 
thermodynamic processes undergone by the fluid in the 
three flow regimes. 


A nomenclature list appears at the end of paper. 












ENTRANCE REGION 


The relation between the enthalpy of the flowing 
fluid and the pressure and temperature of this fluid 
is given by 

me ere + eee. Cf ee (2) 

Equation 2 is applicable regardless of whether the 
flowing fluid exists as a gas, a liquid, or as a gas-liquid 
mixture. Liquid is the only fluid phase present in the 
entrance region. 

In this analysis the effect of the pressure on the en- 
thalpy of this liquid is neglected, so that dH = c,dT is 
assumed to apply in the entrance region. If this value is 
substituted for dH in Equation 1 and the result inte- 
grated with the boundary conditions specified as T = T 


d 1T 
at x = O and~ =(5 at x = x,, the temperature 
dx dx J: 


distribution for the entrance region will be found to be 


r=1.+(=) (: = ‘). ees (3) 
dxJ, ae**1 


TRANSITION REGION 


In the original paper’ the boiling point at x, was 
considered to be a mathematical sheet where vaporiza- 
tion began and the pressure gradient became discontinu- 
ous, the ratio of this gradient at the upstream side of 
the boiling front to that at the downstream side being 
estimated as equal to (k,/k)., or about 0.8. This con- 
cept of the boiling front was based on the conclusions 
of Wyckoff and Botset’ regarding their critical “equilib- 
rium gas saturation.” Accordingly, the gas content at x, 
would be taken as this equilibrium value, but an anom- 
aly results. 

Vaporization occurs if both phases are flowing, and 
then c, > o . If only the liquid is flowing c, is 
finite, however, and refers to liquid approaching the 
saturation state. As gas and liquid are both present in 
the postulated transition region, even though only one is 


dP an: 
— = 7 ——, in which 
dx "dx 

» is the slope of the vapor pressure curve. Thus, any 
discontinuity in the pressure gradient would be reflected 
in the temperature gradient. But a discontinuity in the 
temperature gradient is not plausible, because the rate 
dT 
dx 


flowing, Equation 2 may be written 


of heat transfer by conduction is proportional to 


and the energy-balance basis of Equation 1 reveals no 
heat sources that could bring about abrupt changes in 


dT 
dx” 


equilibrium gas content over zero length of sand column 
is not tenable physicaliy, the gas content is assumed to 
increase gradually in the direction of fluid motion from 
zero at x, to the critical equilibrium value at x,, a finite 
distance farther downstream. 


For this reason and as attainment of the critical 


As no gas moves within the transition region, Fig. | 
depicts the thermodynamic process within this region 
and the saturated liquid line as coincident. Thus, the 
mass fraction of gas in the flowing stream is represented 
correctly as zero even though the in-place gas content 
is finite, this gas having been formed during the tran- 
sient flow period preceding the establishment of steady 
flow. 
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Fic. 1 — THERMODYNAMIC PROCESSES UNDERGONE BY 

A SINGLE-COMPONENT FLUID THAT VAPORIZES WHILE 

MOVING THROUGH A COLUMN OF THERMALLY-ISOLATED 
Porous MEDIUM. 


Darcy’s law, the saturation relation of pressure and 
temperature, and the transition region concept indicate 


that at x, 
lfd7 dP W wr 4 
oo hg 9 lle > oe Ce 


The first and last members of Equation 4 may be used 


dT : ce Ra , 
to evaluate ( ] ) and the result substituted in Equation 3 
GX fa 


to give the temperature difference between the inlet 
end of the column and x,: 


] W wr ) ee 1 (5 
n Pen >) 
( kA ( ae** 


The pressure difference over the same length of column 
is given by integration of Darcy’s law: 


P P 7 (6) 


The length x, of the entrance region in which subcooled 
liquid moves is that which satisfies Equations 5 and 6 
in conjunction with the saturation relation P, P, (7). 
In the transition region all states are saturation states 
and the presence of stationary gas results in the per- 
meability of the sand to liquid being reduced from k at 
x, to about 0.8 k at x. Because of the saturation con- 
dition and because dH c, dT for liquid, approaching 
the saturation state Equation | may be written 
dP KA dP 
dx We, dx 
and differentiating Darcy’s law for the liquid with re- 
spect to x gives 


(7) 


d'P W uv dk | dP dk 
dx Ak? dx k, dx dx 


(8) 


If the righthand side of Equation 8 is substituted for 


&P , “ 
ie in Equation 7 and the result integrated the length 
ax 


of the transition region is found to be 
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KA 1 (ki); 
We, r (ky). 
where /n [(k,),/(kz) 2] = 0.22. 


ExiT REGION 

The vaporization which begins at x, initiates the exit 
region in which the fluid flows in two phases. In this 
region the enthalpy of the flowing mixture is given by 


Qn y= (9) 





WH=WH,+Wih . ... =. (10) 
Integration of Equation 1 with the boundary conditions 
1T 
specified as H = H, when Ad = {<‘ yields 
dx dx J: 
WH = KA\2 — = + WH, .. €14) 
dx dx }» 


Equating the righthand sides of Equations 10 and 11 to 
each other and setting W, = fW and W, = (1 — f) W, 
results in the energy equation for the exit region: 


_KA|dT (dT 
dx 2 


fH, + (l— f) A. ==, — 
ES S| SS ere 


Darcy’s law for the exit region, in terms of the quantity 
F used in the original presentation, may be stated 


dP Boe eo d within thi a ee dP 
z -— AF’ an 1thin is region de n a; « 
Hence, Equation 12 may be put intw the form 
= ( ] l ) 
H, — H, K\F F, 
SED - neiteendiaaaene 1 ne ase ws 13 
f H, ae Hy, H, H, ) 


When the second term in brackets approaches zero, 
Equation 13 approaches that which would apply for 
isenthalpic flow. 


APPLICATION OF RESULTS 


Data pertaining to “Experiment V” of the earlier 
paper are: W = 2.9 X 10° lb sec’; k (k;), l 
10” ft; A = 2.0 X 10° ft’; (k.). = 0.8 (k,),; K 
2.4 X 10° Btu sec” ft* °F’ (Herbert et al’). The prop- 
erties of the sub-cooled propane approaching the satura- 
tion state at x, were yn = 3 X 10° °F Ib” ft* (Sage 
et al‘); (uv), = 6.2 X 10° ft sec; c, 0.68 Btu Ib” 
°F” (Sage et al‘). As per Fig. 3 of the earlier paper, 
7» = 3 X 10° °F Ib” ft” for the regions of the sys- 
tem in which the propane existed in two phases, and 
according to Fig. 6, F and F, were of the order of 5 
10° ft* sec* and 13 X 10° ft’ sec”, respectively; H. 
H,, was of the order of 10 Btu Ib”. 

Inspection of Equation 5 with the appropriate numer- 
ical data indicates that the quantity a was so large (4.1 
< 10° ft’) that the reduction in temperature in the 
entrance region was negligible, and the postulate of 
isothermal flow in that region sustained. 

Equation 9 indicates that the length of the transition 
region was 5.4 X 10° ft or approximately equal to the 
size of the openings in the No. 80 sieve of the Standard 
Tyler Sieve Scale. About half the sand comprising the 
column in the propane flow experiments was coarser 
than Tyler No. 80, so that the computed length of the 
transition region probably did not exceed at most a 
few pore diameters. Neglect of this region in the earlier 
paper therefore was justified, although the comments 
made there regarding discontinuities in the pressure 
gradient were not justifiable in the limit. 

Finally, the second term in brackets in Equation 13 
may be computed with the foregoing numerical data 
to be of the order of 0.001, which is negligible. 
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These numerical calculations disclose that the heat 
conductive effect in the propane flow experiments was 
of no appreciable consequence. 

The energy-balance principles applied mathematically 
in this note would apply also to multicomponent fluid 
flow systems. Research on such systems, as an exten- 
sion of the present work, appears desirable in the inter- 
est of gaining a better understanding of heat transfer 
in oil reservoirs. 


NOMENCLATURE 


distance along column of porous medium 

rate of heat transfer by conduction 
weight rate of flow of gas plus liquid 

thermal conductivity of fluid-filled porous 

medium 
bulk cross-sectional area of column of porous 
medium 
I = absolute temperature 
= enthalpy per unit weight of fluid 
= Joule-Thomson coefficient (also slope of vapor 
pressure curve) 

c, = heat capacity at constant pressure 

P = fluid pressure 

j. = absolute viscosity 

v = specific volume 
k., = effective permeability to gas 
k,, = effective permeability to liquid 

k = homogeneous-filled permeability 

f = weight of gas per unit weight of gas-liquid 
mixture passing a transverse section of sand 
column in unit time 


II 


ATOx 
i 


I 


aN 


Subscripts: 
G = gaseous phase 
L = liquid phase 
0 = inlet end of column 
| = upstream end of transition region 
2 = downstream end of transition region 
e = critical equilibrium gaseous-phase content o! 
porous medium 
Combination terms: 


_ GW 
wipe | 
F = F(P) = (=A + =) 
PLY PeVa 
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‘Way out in West Texas where the big core grows, Core 
Lab has put into service a fleet of Freeze Trucks. Com- 
pleted only recently, these sparkling new dry ice units 
are already in great demand. 


In 1949, before the advent of widespread diamond coring 
with routine recovery of up to fifty feet of large diameter 
core, Core Lab inaugurated its Frozen Core Service with 
the now-familiar Freeze Box. A supply of dry ice soon 
became standard “equipment” in each local lab. 


Today’s Core Lab Freeze Truck provides special five-foot 
tempered aluminum trays for receiving core on the rig 
floor. The trays are then placed in proper sequence in the 
truck-mounted, insulated chest (containing 250 pounds of 


dry ice) and rapidly transported to the nearest lab where 
analysts begin processing samples directly from the con- 
venient containers. The overall result— more efficient 
preservation of critical fluid contents and less handling 


of the core. 


From thirty-one installations covering all active explora- 
tion areas, Core Lab keeps an eye out for ideas that keep 


services “up to snuff”. 


Whether it’s core analysis, well logging, reservoir fluid 
analysis, or engineering and consulting, Core Lab is con- 
tinually proving out techniques and procedures to do your 


next job even better. 
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An ANALOG COMPUTER for STUDYING HEAT TRANSFER 
during a THERMAL RECOVERY PROCESS 


ABSTRACT 


A design is presented for an elec- 
trical analog computer which can 
solve non-steady state heat transfe) 
problems in an extensive radial for- 
mation containing a moving cylin- 
drical source. The computer is used 
to simulate a simplified thermal oil 
recovery process in which heat trans- 
fer from a moving, constant tem- 
perature source takes place radially) 
by conduction only. Temperaturé 
distribution curves are shown for 
several different assumed modes of 
travel for the heat source. The data 
are used to estimate the residual 
fuel requirements necessary to main- 
tain a_ self-propagating isothermal 
front for the particular system being 
studied. 

Although the computer is design- 
ed to represent an unique system, 
conversion factors can be adjusted 
to show the effects of changing the 
assumed values of system constants, 
such as thermal conductivity of the 
formation and temperature of the 
source isotherm. Examples are given 
to illustrate the relative influence the 
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assumed values of these variables 
will have uvon the quantity of fuel 
required to maintain the source. 

The maximum effect of heat 
transfer by the flowing air stream 
is estimated mathematically, and 
shown to be markedly reduced. 

Inasmuch as the examples discuss- 
ed refer to a simplified process, the 
results are not directly applicable to 
a practical field operation. It is be- 
lieved, however, that the data illus- 
trate some general trends which are 
important in thermal recovery proc- 
esses. 

Data laboratory and _ field 
experimentation can be used to 
modify the computer to take into 
account the influence of heat trans- 
fer by various mechanisms, for ex- 
ample, by injected air and by fluids 
ahead of the hot zone. With such 
refinements it would be nossible to 
estimate more accurately the limit- 
ing values for the rate of travel 
of the high temnerature front and 
for the required air injection rates. 


from 


Use of this computer should aid 
in evaluating the economic feasibil- 
ity of oil recovery by various 
thermal processes. It is honed this 
paper will stimulate further work 
by others to helo accomplish this 
objective. 





UNION OIL OF CALIFORNIA 


INTRODUCTION 


As the costs of exploration fo 
and development of new oil field 
rise, increasing the recovery of o! 
from established fields becomes more 
and more important. In fields con 
taining low gravity crude oil, ult 
mate recoveries by primary mean 
are sometimes less than 10 per cent 
and often total only 10 to 20 pe 
cent of the oil in place. Conventiona 
secondary recovery methods are not 
always effective in increasing th 
total. Methods which have beer 
proposed for increasing the recove! 
of viscous oils include :hermal prox 
esses for increasing the mobility 
the oil. One method is to provide 
moving heat source generated with 
the formation to heat the oil and 
flowing gaseous medium (for e 
ample, air) to sweep the oil to th 
producing wells. The possibility © 
such a process has intrigued en: 
neers for many years,” ** but litth 
or no information has been availab! 
in the literature relative to the tech 
nical or economic feasibility of thi 
method. Recently, however, interes! 
has been revived by publication o! 
laboratory and field test results ind 
cating the feasibility of maintainin 


References given at end of paper. 









a self propagating heat source in an 
oil formation. Very high ultimate 
recoveries were indicated to be pos- 
sible***. 

To date, numerous laboratory 
tests have been devoted to studying 
the many variables involved, such 
as oil and oil-sand characteristics, 
pressure, and air flux. Of necessity, 
such tests are usually made in linear 
cells. The application of laboratory 
results to a radial system of suffi- 
cient scale to evaluate the economic 
feasibility of a thermal oil recovery 
process is extremely expensive and 
subject to considerable trial and 
error. One intermediate step is to use 
mathematical and analog methods to 
study means of utilizing the experi- 
mental information in the most ef- 
fective manner. 

To a large extent, the evaluation 
of economic feasibility depends upon 
a knowledge of such factors as sweep- 
out efficiency, well spacing, and air 
injection rates and pressures. A study 
of areal sweep of a combustion pat- 
tern has been reported in a recent 
paper by H. J. Ramey, Jr., and G. 
W. Nabor.’ 

It is proposed in the present paper 
to describe a design for an electrical 
analog computer which can calculate 
the temperature profile resulting from 
non-steady state heat transfer from 
a moving source in an extensive ra- 
dial system. A method is presented 
for using such profiles to estimate 
the fuel required to maintain a mov- 
ing heat source. The effect of changes 
in the assumed values of different 
variables is demonstrated through 
changes in scale factors in the com- 
puter. 

In the present study the computer 
was used to simulate a simplified ther- 
mal process in which the single fac- 
tor of heat transfer by conduction 
alone was considered. Inasmuch as 
the nature of the thermal oil recov- 
ery process may predetermine the 
quantity of fuel available for generat- 
ing the heat source, and the heat 
transfer characteristics are dependent 
upon the physical operating condi- 
tions, it may be concluded that the 
optimum mode of travel of the source 
will be affected by additional factors 
not included in the present design. 
Through the utilization of data ob- 
tained in laboratory and field experi- 
ments, the computer described herein 
may be modified to evaluate different 
hypothetical mechanisms for the re- 
covery of oil by thermal means. 


APPARATUS AND PROCEDURE 


Mathematical or experimental so- 
lutions of unsteady state flow prob- 
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lems are often difficult and time 
consuming. A common method of 
obtaining a rapid solution to such 
problems is often found by use of 
electrical analogs.*’”” One such net- 
work” has been adapted in our labo- 
ratory to study heat conduction from 
a moving, cylindrical heat source in 
the earth. 

The electrical model makes use of 
the similarities which exist between 
the flow of heat in a rigid body and 
that of charge in a non-inductive elec- 
trical circuit. 

The method is based on the iden- 
tity in form between the following 
electrical and heat flow equations: 


Heat Flow Electrical Flow 
7 
= A R, — AV 
q I 
C. = 2 C.= _Q. 
AV ave 


R, = 


AT ave 

It follows from this identity in 
form of the fundamental defining 
equations that all mathematical re- 
sults based on them will be similar. 

In the analog computer discussed 
in this paper, the simulated forma- 
tion was divided into 100 concentric, 
cylindrical tubes of unit thickness, 
each varying from the next one by 
an increment of radius. Because of 
the symmetry of this two-dimensional 
problem, it was possible to simplify 
the electrical network by adjusting 
the values of the components to rep- 
resent the thermal capacitance and 
resistance associated with the radial 
location of each thermal element. To 
approximate an infinitely large me- 
dium, the last element in the circuit 
was designed to correspond to a for- 
mation volume which was large com- 
pared to the entire radial system 
being studied. 

In order that the computer could 
be operated at a convenient voltage 
and transient time interval, scale fac- 
tors were chosen for the conversion 
from thermal to electrical units, as 
follows: 

1= 0.18 v/°F; that is, 180 v is 
equivalent to 1000°P. 
0.006 sec/hour; that is, one 
minute in the computer is 
equivalent to 10,000 hours 
actual time. 
3.04 * 10” farads/Btu/°F: 
that is, one microfarad is 
equivalent to 3290 Biu/°F. 

The following thermal constants 
were assumed to be representative of 
a fluid saturated, porous rock: 


thermal conductivity, k = 1.6 
Btu/sq ft hour °F/ft 

specific heat, c = 0.2 Btu/Ib °F 

density, p = 140 Ib/cu ft 

diffusivity, « = 0.057 sq ft/hour 


The scale factors may easily be 
changed without invalidating the re- 
sults as long as the relationships be- 
tween the two systems (electrical and 
thermal) are maintained. It is possi- 
ble, within certain practical limits, to 
study the effect of changing the as- 
sumed values of the variables in the 
thermal systems by making suitable 
changes in the scale factors. This is 
important in considering the useful- 
ness of the computer. 

For a given position (x) in the 
thermal system, the thermal resist- 
ance is: 

l 
k A,/ Ar, 
0.625 4" F hour 
A. Btu 


fo convert to electrical units: 


\R, 


fl 
R AR 
im 
_ 
ohms 
| 


x 


1.233 X 10 


where A 
area. 
For the same element, the thermal 
capacitance is: 
AC a(? rh i 6<¢ 
r) Btu/°F 
To convert to electrical units: 
¢.. mAC 
2.66 10° (r 


is the logarithmic mean 


87.6 (r 


r) faradays 

A schematic diagram of the com- 
puter circuit is shown in Fig. 1. The 
body of the computer, a conventional 
RC network made of non-inductive 
resistors and oil-filled condensers, is 
shown with each of the elements 
connected to the bars of a large 
commutator. The voltage source is a 
180 v “B” battery, chosen to give full 
scale deflection of the oscilloscope 
beam when connected directly to the 
deflection plates 

To simulate the movement of a 
heat source, the voltage source is 
connected to a rotating brush which 
contacts the commutator bars. A sec- 
ond brush, rotating at a high speed 
relative to the primary voltage source, 
is used to transfer the condenser volt- 
age at each element to the oscillo- 
scope plates. By means of an adjust- 
able selector switch synchronized 
with the source rotor, a triggering 
signal is sent to a commutator bar 
adjacent to the last RC element to 
initiate the oscilloscope sweep. This 
results in a visual representation of 
the voltage profile in the computer 
network for the selected location of 
the simulated heat source. A photo- 
graph of the assembled computer and 
associated equipment is shown in 
Fig. 2. 
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Fic. 1 — SCHEMATIC DIAGRAM Of 


BY CONDUCTION IN A 


To operate the computer, the rate 
of travel of the source rotor is se- 
lected by changing the gear ratio on 
the drive motor. After a suitable 
warm-up period, the positions of the 
trigger selector switch are adjusted 
and synchronization is checked with 
the oscilloscope sweep circuits. The 
camera is attached to the oscillo- 
scope and the network voltage pro- 
file photographed at the preselected 
positions. 

If the computer is used to simu- 
late a constant temperature source 
moving through a radial formation, 
a temperature profile similar to that 
shown in Fig. 3 is observed on the 
oscilloscope screen. The vertical dis- 
tance between grid lines represents a 
temperature difference of 100°F. The 
curvature of the grid lines calibrates 
the beam deflection for the entire 
screen. The voltage of each segment 
of the computer network is shown as 
a spot on the screen at a vertical 
position corresponding to its numeri 
cal value. 

This analog computer can be 
adapted to the solution of heat trans- 
fer problems involving steady state 
or non-steady state conditions and 
moving or stationary heat sources. In 
order to check the network design, 
analog solutions were determined for 





Fic. 2—-ASSEMBLED COMPUTER AND 
ASSOCIATED EQUIPMENT. 
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problems involving a stationary heat 
source for which the theoretical 
curves could be obtained by conven- 
tional mathematical methods. In gen- 
eral, it was concluded that the net- 
work satisfactorily represented a cy- 
lindrical system and that the solution 
of non-steady state problems in a 
simulated infinite cylindrical medium 
was feasible for short times. For 
problems involving long computing 
times at a given position of the source, 
the error increases with time.” For 
such problems the accuracy could be 
improved by using a DC amplifier at 
each lump in the analog.””” How- 
ever, for problems involving a mov- 
ing source, the large additional ex- 
pense involved in such a refinement 
was not considered necessary. 

As described herein, the computer 
has been designed to study heat con- 
duction problems involving a mov- 
ing, constant temperature source in 
a radial formaticn. It should be 
pointed out, however, that the gen- 
eral design is applicable to a variety 
of heat transfer problems involving 
a moving source. With » me modifi- 
cation it would be possible to study 
non-steady state heat transfer by 
gases, liquids, and solids, either in- 
dividually or as an interacting sys- 
tem, in media of different geomet- 
rical shapes. In addition, the moving 
source rotor may be designed to rep- 
resent a heat source with any of a 
number of unique properties. 


ILLUSTRATIVE PROBLEMS 


To illustrate the type of problem 
which can be studied with this com- 
puter, several assumed modes of 
travel of the heai source will be dis- 
cussed. It should oe emphasized that 
these problems are intended to be 
illustrative only and are not neces- 
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sarily representative of a practical 
field operation. 

For the following illustrations, 
is assumed that a constant temper: 
ture heat source is produced in th« 
formation. Although a number of! 
possible mechanisms for the produc 
tion of such a heat source may be 
visualized, for convenience a mecha 
nism similar to that postulated by 
Kuhn and Koch* has been assumed 
According to their hypothesis, a com 
bustion zone is produced by burning 
residual carbonaceous material d: 
posited out of the native crude oi! 
ahead of the high temperature zone 
In the idealized case, the tempera 
ture of this combustion zone could 
be maintained at a constant optimu 
value by adjustment of the rate « 
air injection into the formation 
Although reference will be mad 
throughout the ensuing discussio! 
to the “combustion front” and th 
associated fuel requirements, any 
means of obtaining a constant tem 
perature front would result in ident) 
cal heat transfer characteristics fo 
the particular conditions postulated 

In one set of tests, a constan/ 
temperature source was moved radi 
ally through the formation at variou 
constant speeds which ranged fron 
0.874 ft/day to 0.097 ft/day. Tem 
perature profiles were photographed 
at convenient positions. In Figs. 4 
5, 6, and 7, temperature data fron 
the above tests have been plotted o1 
a semi-log scale to show the com 
puted temperature distributions fo 
different positions of the heat source 
It is immediately apparent that (1 
at a constant rate of travel, preheat 
ing of the formation ahead of th 
source isotherm increases with di 
tance of travel and (2) at a given 1 
dius, preheating increases for slowe 
rates of travel. 


These temperature distributio 
curves may be used to calculate th 
theoretically required amount of res 
dual fuel to maintain a self-propagat 

















Fic. 3—TYPICAL VOLTAGE PROFIL! 
IN ANALOG COMPUTER NETWORK 
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Fic. 4— TEMPERATURE DisTRIBU- 
TION AHEAD OF COMBUSTION FRONT. 
Conduction only in radial for- 
mation. Rate of movement of 


front: 0.874 ft/day. 


ing combustion front at a given rate 
of advance. In order to make such a 
calculation, it is necessary to know 
the instantaneous heat loss from the 
combustion front. If this heat loss 
can be supplied by burning the avail- 
able fuel, the temperature of the 
combustion front will remain con- 
stant and the source will be self-per- 
petuating as long as sufficient air is 
supplied. If there is insufficient fuel, 
the temperature will decrease and, 
under certain conditions, the process 
will cease. 

The instantaneous heat loss can 
be determined from the equation, 
q kA 0T/ér. The term, @7/@r, is 
simply the instantaneous temperature 
gradient at the heat source and can 
be readily determined directly from 
the figures as 

In T, 


In T. 7 


T /ar 


Then the fuel required to maintain 
the combustion front will be 
k aT/ér 


s dr/dt 


WwW lb/cu ft of formation 
where s is the heat of combustion of 
the residual fuel (assumed as 17,000 
Btu/Ilb) and dr/dt is the rate of 
movement of the front in feet per 
hour. 

Through use of the above method, 
Fig. 8 has been plotted to show the 
residual fuel required to maintain a 
moving heat source (temperature 
1,000°F above ambient) in a forma- 
tion in which heat is transferred ra- 
dially by conduction alone. 

A suitable change in conversion 
factors (/ 0.18, n 0.006, m 
6.08 XK 10”, r rx/V2), can be 
made to show the effect of changing 
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Fic. 5——- TEMPERATURE DIsTRIBU- 
TION AHEAD OF COMBUSTION FRONT. 
Conduction only in radial for- 


mation. Rate of movement of 
front: 0.353 ft/day. 


the assumed value of the conduc- 
tivity. In Fig. 9 the fuel requirements 
to maintain a moving combustion 
front are indicated for k = 0.8, A? 

1,000. 

Similarly, through the choice of a 
different set of conversion factors 
(/ 0.30, n 0.006, m = 6.08 
i. r./\/2), the effect of a 
change in the temperature of the 
combustion front can be found. Fig. 
10 shows the fuel requirements if 
k 0.8, AT 600. 

In general, Figs. 8, 9, and 10 show 
that the fuel requirements are higher 
for the slower rates of travel of the 
source. Inasmuch as the physical na- 
ture of the assumed process makes 
available only a given amount of 
fuel, which may be affected by for- 
mation and oil characteristics, there 
should be a minimum rate of ad- 
vance possible for a self-propagating 
combustion front. The curves indi- 
cate that this critical rate will de- 
crease with distance from the well- 
bore. This characteristic is a result 
of the greater preheating by conduc- 
tion as the combustion front travels 
away from the wellbore and should 
be true when the additional effect of 
iniected air is considered. 

From the data presented above. it 
is apparent that the differences be- 
tween the actual value of k and that 
used in the computer are not critica! 
in determining the reauired fuel con- 
tent. provided & is of the right order 
of magnitude. Of greater importance 
in using the data is the value as- 
sumed for the temperature at the 
combustion front, inasmuch as the 
required fuel content was shown to 
be directly proporiional to the tem- 
perature. This indicates one area in 
which experimental data from labo- 
ratory tests would be desirable. 


In order to maintain a constant 
advance of the combustion 
front as shown in these figures, it 
would be necessary to provide for an 
increasing volume rate of injection 
air as the front travels away from 
the wellbore. To determine the re- 
quired air rates, it would be neces- 
sary to know the position of the 
front at any time. This is not likely 
to be practicable for field operation. 
A simpler operational method would 
be to inject air at a constant rate. If 
both the fuel supply per unit volume 
of formation and the air injection 
rate are assumed to be constant, it 
would appear that because of geo- 
metrical factors the minimum rate of 
movement of the combustion front 
should be inversely proportional to 
the radial position. 

It is possible to study the feasibil- 
ity of this mode of operation by 
either of two methods. In one case, a 
family of curves, each of which is 
radius, can be ob- 


rate of 


proportional to | 
tained with a constant rate drive. In 
the second case, the computer can be 
set up so that the source rotor can be 
moved at a rate proportional to 
1/radius 

To illustrate the first method, the 
family of curves in Fig. 8 has been 
used to examine a mode of travel in 
which the rate of movement of the 
combustion front is equal to 8/ra- 
dius in feet per day. Then for the 
constant rates of 0.874, 0.353, 0.211 
and 0.097 ft/day, the front would 
be located at approximately 9, 23, 38 
and 83 ft, respectively, in order to 
he traveling at the designated rate. 
If these points are plotted on Fig. 8 
they will all fall at a constant fuel 
requirement of approximately 2.3 
lb/cu ft of formation, indicating that 
the front would be self-propagating 
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under the above conditions. 

If the computer is set up so that 
temperature profiles are obtained at 
radii, r,, at the corresponding rate, 
8/r,, the required residual fuel can 
te calculated from the computed 
temperature gradients in the manner 
described previously. In this case a 
single curve would be obtained as 
shown in Fig. 11. Here again the re- 
quired residual fuel is found to be 
2.3 Ib/cu ft of formation. 

In general, it may be concluded 
that for this mode of travel the nu- 
merical value of both the tempera- 
ture gradient and the radial velocity 
of the combustion front must be 
changing at the same rate. The results 
indicate also that in the radial sys- 
tem described a constant rate of fuel 
consumption will result in a constant 
temperature heat source only if the 
rate of movement is inversely pro- 
portional to the radial position. 


In the above discussion, heat trans- 
fer by the injected air has been neg- 
lected. It is possible, however, to use 
the available data to estimate the 
minimum fuel requirements in a sys- 
tem in which no heat is retained be- 
hind the combustion zone. This can 
be done by assuming that all of the 
heat stored in the rock as the mov- 
ing source reaches and passes a given 
radius will be transferred instantane- 
ously by the flowing air stream back 
to the combustion front. Thus the 
minimum fuel requirement in such a 
system can be approximated simply 
by subtracting the equivalent fuel 
content of the corresponding volume 
increment of heated rock from the 
residual fuel content caiculated from 
the data shown in Figs. 4, 5, 6, and 
7. This has been done for a radial 
thermal system in which k = 1.6 and 
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Fic. 8 — FUEL REQUIRED TO MAIN- 
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Thermal conductivity: 1.6 Btu/ 
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1,000°F above Ambient Tem- 
perature. 


curves have been plotted in Fig. 12. 
The values for the residual fuel con- 
tent shown in Figs. 8 and 12 repre- 
sent two extremes of operation. In 
the one case (Fig. 8), it is assumed 
that all of the heat stored in the rock 
behind the combustion front remains 
unavailable for raising the tempera- 
ture of the formation. In the second 
case (Fig. 12), all of this heat is uti- 
lized to the fullest extent. 

As a final illustration of the use- 
fulness of the analog, there is shown 
in Fig. 13 the analog computation 
of the temperature dec = a radial 
formation after a hea ource has 
stopped moving. In ths problem, the 
combustion front ws. raoved through 
the formation at a constant rate of 
0.353 ft/day to a radius of 40 ft and 
then stopped. This might correspond 
to a system (with heat flow by con- 
duction only) in which air injection 
has stopped because of equipment 
failure. From the slow decay of tem- 
perature shown, it might be con- 
cluded that once a hot zone has been 
well established, re-ignition of the oil 
sand following a prolonged shutdown 
should not be a serious problem. 

This problem is similar in nature 
also to that investigated mathemati- 
cally by Jenkins and Aronofsky.” The 
temperature increase after 840 days 
of 70° F at a radius of 100 ft com- 
pares reasonably well with their theo- 
retical value of 75° F after two years. 

From the above illustrations, it 
would appear that the computed tem- 
perature distribution curves should be 
useful for studying various factors, 
such as the best start-up technique 
and the optimum mode of travel of 
the heat front. It is obvious, however, 
that the simplified mechanism set up 
for the foregoing problems does not 
simulate closely the heat transfer 
process existing in a practical field 
operation. In order to obtain the 
maximum utility from this machine, 
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it would be desirable to modify it to 
include additional factors such as 
heat transfer by the injected air and 
by the movement of hot fluids ahead 
of the combustion front. It seems ap 
parent that modification of the com 
puter to study heat transfer by vari 
ous assumed mechanisms must bh: 
coordinated closely with results ob 
tained from laboratory and field ex 
perimentation. 


CONCLUSIONS 


A design has been presented fo: 
an electrical analog computer whic! 
can solve non-steady state heat con 
duction problems pertaining to an ex 
tensive radial system containing 
moving, cylindrical heat source. 

It has been the intent of this pape: 
to show the type of problem the ana 
log computer can be adapted to 
study, and to indicate how its us 
can be integrated with the result 
obtained from laboratory data to in 
crease Our understanding of possib| 
field applications of oil recovery b’ 
thermal methods. A method has been 
shown for using the temperature dis 
tribution curves obtained with the 
computer to calculate the fuel | 
quired to maintain a self-propagatin 
heat source. 

Inasmuch as the computer desig: 
presented herein takes into account 
only those heat losses resulting fron 
radial conduction, the data 
not be construed as applying direct! 
to field operations. Modification o!| 
the computer to include other heat 
transfer mechanisms depends pri 
marily upon an understanding o! 
the particular thermal! oil recovery 
method being studied. Laborato: 
scale studies of the thermal proces 
itself aid in determining the mos! 
satisfactory mechanism to use. 


Because of the complex nature o! 
thermal oil recovery processes. i! 


should 
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would appear that continued efforts 
of research workers in many differ- 
ent groups will be required to eval- 
uate the numerous possible ramifica- 
tions in a reasonable time. It is hoped 
that this paper will stimulate the pub- 
lication of other work on thermal 
recovery processes. 
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NOMENCLATURE 


Thermal units 
R,=thermal resistance, 
F hour/Btu 
C, = thermal capacity, Btu/°F 
QO. =heat stored in C,, Btu 
q = hourly heat flow, Btu/hour 
AT = rise in temperature of C, 
due to Q,, °F 


COMBUSTION FRONT MOVES 0.353 FT/DAY 
FOR 40 FT, THEN STOPS 
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Fic. 11—-FUEL REQUIRED TO MAIN- 
TAIN A MOVING COMBUSTION FRONT. 

Thermal conductivity: 1.6 Btu 


ft hour °F/ft. Combusion tem- 
perature 1,000°F above Ambient. 


= distance in the thermal cir- 
cuit (radial), ft 
average cross sectional area 
of each increment Ar,, 
sq ft 
- number of the element 
(Ar) under considera- 
tion 
time in the thermal circuit, 
hour 
k ag 
thermal diffusivity, 
sq ft/hour 
heat capacity per unit 
mass, Btu/Ib °F 
= thermal conductivity, Btu 
sq ft hour °F/ft 
p = density of formation, Ib 
cu ft 
Electrical units 
R, = electrical resistance, ohms 
= electrical capacity, farads 
Q. = charge stored in C,, cou- 
lombs 
electrical current, amperes 
- rise in voltage of C., due 
to Q., volts 
time in the electrical cir- 
cuit, seconds 
= distance in the electrical 
circuit (radial), ft 
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The authors are to be compli- 
mented for their ingenious solution 
to a problem having an important 
bearing on future development of 
thermal oil recovery processes. The 
interpretation of results from this 
heat transfer study relative to an in- 
situ combustion oil recovery process 
illustrates the utility of this type of 
theoretical work. Undoubtedly, their 
electric analog computer will have 
many other uses in reservoir tech- 
nology. 

With regard to the results of the 
study, it should be emphasized that 
for the more conventional well pat- 
terns, radial flow conditions exist for 
only 20 to 30 per cent of the distance 
from an injection well toward the 
nearest production well. The results 
from this study are, therefore, most 
applicable and important in consid- 
eration of the ignition phase of this 
type of thermal oil recovery process. 
As the authors implied, the fuel re- 
quirements presented are based on 
assumed thermal properties and an 
idealized mechanism. Actual fuel con- 
sumption can differ greatly from the 
theoretical values obtained from this 
study. Nevertheless, the trends illus- 
trated by the authors permit impor- 
tant observations. 

For the authors’ case of a constant 
combustion front velocity, the mini- 
mum fuel requirement appears espe- 
cially critical when the burning front 
is near the injection well and/or 
when the front is moving at a slow 
rate. In contrast, the results for the 
case where the velocity varies in- 
versely with radius of the burning 
front indicate a fuel requirement that 
does not vary with position of the 
burning front. While this result tends 
to contradict the conclusions drawn 
from the case of constant velocity, 
it is believed that most actual opera- 
tions will fall between these condi- 
tions, and the fuel requirement wil! 
be substantially higher near the in- 
jection well. 
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During investigation of this ther- 
mal recovery process, we computed 
fuel requirements by an analytical 
solution of substantially the same 
heat conduction problem. Two cases 
were considered: (1) heat conduction 
from an infinite cylindrical source 
moving radially in an infinite medium 
of constant thermal properties, and 
(2) the same problem with the in- 
clusion of the boundary conditions 
resulting from consideration of the 
ignition wellbore. 

The solutions to Case (1) for a 
burning front velocity inversely pro- 
portional to the radial position of the 
front showed the fuel requirement to 
be entirely independent of position 
of the front. And for the same ther- 
mal conditions assumed by Vogel and 
Krueger, computed results were in 
good agreement with the values from 
the electrical analogy. A fuel require- 
ment of 2.4 Ibs/cu ft was obtained 
as compared with 2.3 reported by 
Vogel and Krueger. This establishes 
a direct verification of the electrical 
analog for a moving source. 

A solution for Case (2) at a radial 
position of 0.5 ft indicated a slightly 
higher fuel requirement of 2.7 Ibs/cu 
ft. But this is attributed to the effect 
of the wellbore. 

The analytical solutions mentioned 
above are as foliows. 

1. Case (1) — Constant velocity 


of source (7 ry. 
T-T 
. ? x (t t’) 
2k V(t ?’) ‘ 
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Where 7, is the initial temperature, 
T is the temperature at radius r at 
time ¢, q is the rate of heat genera- 
tion per square foot of frontal area, 
U is the constant velocity of the 


burning front, k is the thermal con 
ductivity, ¢ is time, and ¢ is the time 
the burning front is located at radius 
r’, « is the thermal diffusivity, and 
I, is Bessel’s function of the first kind 
of an imaginary argument. If Equa 
tion | is integrated for r r’ att 

the fuel requirements may be found 
from the expression: 
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where s is the heat of combustion of 
the fuel. 

2. Case (1) — Velocity of source 
is inversely proportional to the radius 
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The fuel requirements may be found 
by integrating Equation 3 for r } 
at ¢ and using the expression: 
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In Equations 6 and 7, ¢ is the spe- 
cific heat of the saturated sand, p is 
bulk density, r.. is the radius of the 
wellbore, J, and Y, are Bessel func- 
tions of the first and second kind, 
respectively, of order zero, and x is 
a variable of integration. 

A further refinement of the prob- 


AUTHORS’ REPLY to JENKINS and 


We should like to thank Rodman 
Jenkins and H. J. Ramey, Jr., for 
their pertinent comments. It is grat- 
ifying to see the good agreement re- 
ported between the analyiical and 
the analog solutions to the prob- 
lem of determining the fuel require- 
ment in a radial system in which 
the burning front velocity is in- 
versely proportional to the radial 
position. This same agreement has 
been noted by R. G. Hawthorne of 
our laboratory. 

The presentation of the equations 
used in studying the simplified prob- 
lems discussed in the paper serves 
to illustrate the somewhat compli- 
cated nature of analytical solutions 
to heat transfer during thermal oil 
recovery processes. Equations 3, 6, 
and 7 illustrate the complexity added 
to the simple case by considering 
the wellbore effects. In the com- 
puter, wellbore effects are au- 
tomatically taken care of in the 
design of the electrical lumps near 
the wellbore. It would be interesting 
to see the analytical solutions to the 
more complex heat transfer problem 
involving the additional factors of 
vertical heat losses and heat trans- 
fer by the driving gases and con- 
densing fluids. 


lem treated by authors Vogel and 
Krueger is the consideration of verti- 
cal heat losses. This would require 
the expensive modification of rebuild 
ing the RC network in two dimen 
sions, and considering a burning front 
of finite vertical dimension. Vertical 
heat loss can be important to the 
maintenance of self-sustained com- 
bustion when the sand section is thin 
and the well spacing is large. If ver- 
tical heat losses were considered, it 
would also be desirable to account 
for the lower thermal diffusivity of 


In general, analytical solutions of 
this nature are tedious and time con- 
suming unless attacked with modern 
digital computer techniques. In any 
case, considerable expense is in- 
volved in arriving at numerical an- 
swers. It would be of considerable 
value to the industry to have further 
work published on solutions to the 
more complex problems. On _ the 
basis of our present thinking, we 
feel that most such problems can be 
treated somewhat more simply and 
inexpensively by the modification of 
the analog computer. 

Jenkins and Ramey have sug 
gested that vertical heat losses be 
considered as a refinement to the 
problems treated in the paper. As 
they have noted, vertical heat losses 
can be important when considering 
the problem of maintaining a com- 
bustion front in thin oil zones. We 
believe, however, that the applica- 
tions in which vertical heat losses 
are the controlling factor are likely 
to be the exception rather than the 
rule. For thicker zones there is some 
evidence from analytical studies that 
the effect of vertical heat transfer 
on fuel requirements is considerably 
less than the effect of heat trans- 
fer by injected air (or inert gases) 


the “burned” sand section of the 
reservoll 

We agree with the authors’ state 
ment that the developed by 
paper should not be 
construed as applying directly to field 
Operations, inasmuch as the front is 
Known to progress over a range of 
temperatures, velocities, and fuel con- 
sumptions in accordance with incom- 
pletely defined rate phenomena oc- 
curring in a reservoir. Many chal- 
lenging problems along these lines 
remain to be solved. 
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and possibly by the flowing liquids. 
both of these latter 
factors can be studied rather easily) 
by inexpensive modifications of the 


The effects of 


analog compute! 

As Jenkins and Ramey have em- 
phasized, radial heat flow conditions 
exist for only 20 to 30 per cent of 
the distance from the injection well 
toward the producing well. At the 
believe this region to 
require the most critical considera- 
tion when attempting to evaluate the 
feasibility of starting a successful 
thermal oil recovery process. How- 
feel that a study of heat 
later stages of the 
achieved through 

modifications of 
appears that ana- 
involving interfer- 
wells might become 


moment, we 


ever, we 
transfer during 
process can be 
relatively simple 
the computer. It 
lytical solutions 
ence between 
quite complex 
It is obvious from a consideration 
of both the comments of Jenkins 
and Ramey and the authors’ reply 
that there is much to be studied and 
learned in this field. Both analytical 
and analog studies can be helpful in 
solving the many problems involved. 
However, as mentioned in the paper, 
their maximum utilization will be 
achieved through correlation with 
laboratory and field studies. wk 
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LOCATING CASING SHOE LEAKS with RADIOACTIVE ARGON 
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ABSTRACT 


Radioactive materials were used underground before 
the advent of the nuclear reactor. Radium containing 
ores have been mixed with cements and other materials 
for later location in oil wells by gamma ray logging 
since 1940.' Perforating bullets have been marked with 
radium salts, and radioactive materials have been added 
to drilling muds to find the cause of “lost circulation.” 
More recently colloidal suspensions of radioactive 
cobalt-60 and silver-110’ and solutions of I-131* have 
been used in waterflooding operations to determine 
injection profiles, flow paths and water velocities. 

In this paper a recent and interesting underground 
tracer application is discussed. Development work and 
field engineering was carried out under the direction of 
Tracerlab, Inc., with assistance from personnel of the 
client. A gaseous tracer, argon-41, was used in an at- 
tempt to determine the cause of leakage from an under- 
ground gas storage reservoir. 


STATEMENT OF PROBLEM 


This gaseous tracer application was carried out for 
the Natural Gas Storage Co. of Illinois. This company 
has develope? an underground gas storage field 1,750 
ft below th. surface of the earth near the village of 
Herscher, Ill., about 60 miies south of Chicago. The 
total volume of gas stored in the field was estimated to 
be about 13 billion cu ft in the spring of 1954. The 
field has a maximum estimated storage capacity of 90 
billion cu ft. The gas is contained in the Galesville sand- 
stone, previously a water bearing formation of about 25 
per cent porosity. The water as it is pushed out by the 





1References given at end of paper. 
Original manuscript received in Petroleum Branch office on July 
8, 1955. Revised manuscript received Nov. 6, 1955. 
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gas forms a barrier for the gas around the edge of the 
reservoir. Over the top and beneath the bottom of the 
reservoir are located several layers of impervious rock, 
which barring cracks or unplugged and unknown wells 
should act as seals. 

Twenty-one injection and withdrawal wells were 
drilled into the Herscher reservoir and gas injection was 
begun during the summer of 1953. After a few weeks 
of gas injection, evidence of gas leakage was detected at 
the surface; water wells began bubbling and losing 
prime, and dry spots in the fields indicated gas leakage 
from uncharted oil wells. In every case, the leakage 
could not be traced back any further than the water 
bearing formations a few hundred feet below the sur 
face. 

An extensive search was made to locate uncharted oi! 
wells which were believed to be possible channels fo! 
the gas from the storage formation to the water contain 
ing formation. Long-time residents in the area assisted 
in the location of several abandoned oil wells. The wells 
when located were cleaned out and either used as vent 
wells to relieve the gas pressure in the water formations, 
or were plugged with cement. None of the wells even 
approached depths which would endanger the soundness 
of the storage reservoir. 

All of the injection wells were tested for leaks in the 
casings by conducting neutron surveys and by othe! 
techniques. 

To detect any leakage up through the cement-filled 
annulus surrounding the casing, a section of the casing 
wall in one well was milled out. Then, plugs were in- 
stalled above and below the cut-out section and pressure 
applied to the annulus, which proved conclusively that 
the cement was sound in that well. The well was then 
plugged with cement and abandoned. In another well, 
ammonia was injected into the well and a temperature 
log was made. It was hoped that a leak would be re- 












vealed by a temperature rise resulting from heat being 
generated by the ammonia passing through the leak and 
dissolving in formation water. 


At the request of the Natural Gas Storage Co. a de- 
tailed laboratory study was made by Tracerlab to deter- 
mine the feasibility of using radioactive gas tracers to 
locate casing shoe leaks. In the course of the study avail- 
able references were consulted” *. Several gaseous tracers 
were considered, and methyl bromide tagged with bro- 
mine-82* was selected as the best available gamma- 
emitting gaseous tracer. The quantity of tagged methy! 
bromide tracer required was determined. Techniques for 
preparation of methyl bromide were worked out* 
Then, methods of monitoring water supplies and a 
method of injection were solved. A schematic diagram 
of the setup devised to locate casing shoe leaks is illus- 
trated in Fig. 1. 


After thoroughly satisfying themselves as to the ab- 
sence of any hazards to residents of Herscher, or to 
their own employees, the Natural Gas Storage Co. of 
Illinois decided to proceed with the tests. Arrangements 
were made with the U.S. Atomic Energy Commission 
to obtain the required Br-82. Equipment was shipped 
to Herscher by Tracerlab for synthesis of the methy] 
bromide and for monitoring, and a test program was 
set up. 
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Fic. 1—EQUIPMENT PLACEMENT FOR LOCATING CASING 
SHOE LEAKs BY RADIOACTIVE TRACE! ER TECHNIQUE. 


*Bromine-82 [Br-81 (n, y) Br-82] .55 to 1.31 MEV y, .465 MEV 8, 
35.5 hr half-life. 

**The methyl bromide was prepared by reacting KBr, H2SO, and 
ROH; distilling overhead the methyl bromide, and condensing the 
methyl bromide in a cold trap maintained at —7°C. 
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PREPARATION OF EQUIPMENT 


Because of its short half-life, argon-41+ had not 
been considered as a tracer for this work. However, 
L. D. Marinelli of the Radiological Physics Div. of 
the Argonne National Laboratory, who had been called 
in as a health physics consultant, suggested that Argonne 
might be able to prepare this isotope for the tests. 
Argonne was able to prepare the A-41, approval for 
this preparation having been obtained from the Isotopes 
Div., U.S. Atomic Energy Commission. Because of its 
chemical inertness, low health hazard, and ease of use, 
A-41 was selected for the tests rather than Br-82. The 
two-hour half-life of argon-41 dictated that the com- 
plete operation, from obtaining the argon-4i from the 
reactor to its injection into each well, must be car- 
ried out on an exacting time schedule 

Irradiation capsules with a 90 cc volume were con- 
structed of aluminum. A chamber for filling the irradia- 
tion capsule with stable argon and for removing the 
radioactive argon was constructed as shown in Fig. 2. 
The method of filling the capsule was as follows: 

The capsule, with the screw cap loose, was placed 
in the chamber. By manipulation of valves, the cham- 
ber and capsule were successively evacuated and then 
filled with 15 psig of argon four or five times. After 
the final filling of the capsule with argon the pressure 
was dropped to atmospheric and the capsule cap tight- 
ened. The capsule was checked for leakage by placing 
it in boiling water and observing if bubbles formed. 

The capsule containing the argon was placed in the 
CP-£ reactor at the Argonne National Laboratories. It 
was subjected t~ a flux of 2.3 x 10” thermal neutrons 
per cm’/sec for . irradiation time from 12 to 16 hours 
(overnight), which produced approximately 1,600 mc 
of argon-41. The reactor was shut down at 8:00 a.m. 
and the ix:.<diated argon was removed after a 15-minute 
cooling off period which permitted decay of the 2.4 
minute A1-28. 

The capsuie in a 350-lb lead shield was transported 
by automobile to the Herscher compressor station. The 
radioactive argon was transferred to the injector (Fig. 
3) in a hood by placing the irradiation capsule (after 
removing the screw cap with a long handled tool) into 
the chamber which was used for filling the capsule. The 
injection tank, having previously been evacuated, was 
connected to the chamber through a manifold as shown 
in Fig. 2. Four in. of lead shielding were placed around 
both the chamber and injection tank. By opening the 
injection tank valve, the radioactive argon was drawn 
into the tank until the pressure equalized. The tank 
valve was then closed and the chamber pressured to 
15 psig with stable argon gas. These operations were 
repeated until a survey meter located near the injection 
tank ceased to indicate an increase in radiation. The 


TABLE 1— TYPICAL SCHEDULE WITH RADIOACTIVITY VALUES FOR EACH 
OPERATION 
Quantity of A-41 


Operation Present 


Remove argon-41 from reactor 1470 mc 


Injector and Detector in Positio 356 me 


Injection timer opens valve to 286 mc 
release tracer 

Injection timer closes valve (186 me injected) 

43 mc total 

(27 mec injected 
16 mc in tank) 


Test conclusion — 
Withdraw injector 


0147 mc total 
0093 mc injected 
0054 mc in tank 


24 hours after injection 


+Argon-41 [A-40 (n, ¥) -41] 1 MEV vy, 109 minute half-life. 
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Fic. 2 —- LABORATORY SETUP FOR TRANSFER OF RADIO- 
ACTIVE A-41 To GAs INJECTOR. 


injection tank was then pressurized with stable argon 
to 2,000 psig. At the time of transfer, the activity of 
the empty capsule was estimated to be over 1 curie 
due to the impurities in the aluminum (principally cop- 
per-64, 12.88 hour half-life), even though the purest 
available commercial aluminum was used for capsule 
fabrication. 


TEST PROCEDURI 


The injection apparatus shown in Fig. 3 comprised 
a 12.3 cu ft, 2,000 psi tank, a pressure regulator which 
determined the injection flow rate, and a timing mech- 
anism which opened and closed the release valve at 
predetermined times. 

The base of the well had been enlarged to present 
more formation area to the injection gas, as shown in 
Fig. 1. This space was packed with gravel to conserve 
the radioactive gas and support the injector. The injec- 
tor, containing the radioactive gas, was lowered to the 
base of the casing shoe in the injection well by use of 
a shooters hook. At the bottom of the well, the shoot- 
ers hook was released from the injection apparatus and 
withdrawn from the well. A gamma ray detector (3-in. 
diameter by 3-ft long high-pressure ionization cham- 
ber) was lowered inte the well and set at the elevation 
where the indicated radiation level was about one-third 
to one-half scale due to direct radiation from the injec- 
tion tank (approximately 12 ft above the injector). The 
arrangement of the detector and injection apparatus 
within the well is shown in Fig. |. The gas flow into 
the well was regulated to restrain the radioactive gas, 
when injected, from flowing upwards in the well. 

At the time of injection, the radiation level at the 
detector decreased, indicating that the radioactive gas 
was moved from the injection tank down to the base 





VOL. 204, 1955 


of the well. Since the well pressure was approximately 
700 psig and the tank pressure was 2,000 psig, only 
two-thirds of the radioactive gas could be injected. 


A typical schedule with radioactivity values for each 
operation is given in Table 1. 


The radiation level approximately 12 ft above the 
casing shoe was recorded for a minimum of four hours 
after injection. At 30 minutes and at 60 minutes after 
injection, the radioactive gas was brought up to the 
detector by venting gas from the well through a positive 
displacement meter. As soon as the radiation level at 
the detector indicated that the radioactive gas was up 
to the detector, a sufficient amount of gas was metered 
into the well to move the radioactive gas back to the 
bottom of the casing shoe. 


It was found that very small gas leaks in the well 
casing or out the top of the wellhead could be quan- 
titatively measured by timing the uphole movement of 
the radioactive gas. In several instances (with all known 
gas inlet valves closed) the radioactive gas was found 
to move down in the well. This indicated a down hole 
flow in the well, caused by leakage through the valve 
used to isolate the injection well from the gas injection 
line. 

Three hours after injection of the tracer, the radio- 
active gas was drawn 150 ft up the casing by venting 
gas from the well. A volume of gas twice the quantity 
of the vented gas was then injected into the well. This 
procedure was followed to assure that the volume of 
radioactive gas was sufficiently large to fill the cavity 
at the base of the well and to increase the probability 
of the tracer moving up any shoe leak. At the time of 
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Fic. 3 — RADIOACTIVE Gas INJECTOR. 









venting the gas, about 30 mc of released A-41 remained 
in the well, This amount was sufficient to give a read- 
ing more than seven times full scale when the gas 
flowed past the detector. 

At the conclusion of each test, a complete gamma 
ray log was made of the well. This log was compared 
with a background log made before the radioactivity 
testing program in order to check for new radiation 
abnormalities. 

Even though the radioactive gas was mixed with 150 
linear ft of casing gas, some question still existed as 
to whether the volume of radioactive gas was sufficient 
in size to assure that some of the tracer would reach 
any leak. To dispel this doubt several wells were tested 
by mixing the radioactive argon with 500 linear ft of 
casing gas, and in other wells the radioactive gas was 
mixed with a full length of casing gas, 1,750 ft or 
approximately 600 cu ft at 700 psi for the 8-in casing. 

To assure a relatively even mixture of gas throughout 
these large volumes, a controlled mixing system was 
devised. The radioactive argon was injected at the top 
of the well by using a purge bottle in place of the 
injector. A metering arrangement as shown in Fig. 4 
permitted regulating the ratio of radioactive gas to 
natural gas. An added advantage of this method of 
introduction was that a larger percentage of the radio- 
active argon could be injected into the well. 

Prior to injecting the radioactive argon into the well 
the gamma ray detector was lowered to a level 5 ft 
above the casing shoe. The radiation level at the detec- 
tor was above full scale throughout the passage of the 
full casing length of gas. 

In the tests where the down-hole injector was not 
used, it was possible to log completely to the bottom 
of the well with the gamma detector. In two instances, 
a high radiation level was observed a few feet above the 
bottom of the casing shoe. In both cases where the 
radioactivity was indicated behind the casing it was 
noted by a study of core samples and elevation data 
that the casing shoe extended into the porous formation 
a few feet rather than terminating at the bottom edge 
of the cap rock as it should. The radioactive gas was 
apparently trapped behind the casing just below the 
cap rock. 


CONCLUSIONS 


All 20 operating injection wells were tested once; 
one well was tested twice; one well tested three times, 
and the well most suspected of leaking was tested five 
times. No indication of leakage was observed in any 
of the wells. Sufficient information was obtained by 
the tests to conclude that there was no leakage up 
through the casing annulus in the vicinity of the casing 
shoe. The loss of gas from the reservoir was apparently 
occurring elsewhere than at the casing shoes. 

The 27 tests were carried out over a period of 34 
months. Most of them were conducted at a rate of three 
tests per week. A total of over 40 curies of argon-41, 
and nearly as much activity in the aluminum capsules 
was handled under field conditions. The radiation ex- 
posure experienced by the field engineers handling the 
radioactivity was under 10 mr per test. The presently 
accepted tolerance level is 300 mr per week or 50 mr 
per day. 
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ABSTRACT 


A practical analytical method is developed in this 
paper which provides the practicing reservoir engineer 
with a handy method for analyzing oilfield interference 
problems. The procedure employed entails the principle 
of superposition of individual effects. A set of charts 
has been developed and can be used to evaluate the 
individual effects. The charts represent non-steady state 
pressure distribution in an extensive aquifer due to a 
constant rate of production from an oilfield which 
can be approximated by circular geometry; they are 
constructed from solutions to the diffusivity equation 
considering a circular sink as distinguished from a point 
sink. The solutions are obtained partially from math- 
ematical considerations and partially from analog and 
digital computations. 

By -virtue of the method of images, the charts can 
be used also to treat the case of production near a 
fault. In addition, provisions are made for treating oil 
fields, which are produced by partial water-drive. 

Treating an actual field case is deferred to a sub- 
sequent writing; however, illustrative example is pre- 
sented to demonstrate the method of application. 

The pressure distribution in linear and limited radial 
aquifers is currently under study. 


INTRODUC™’ON 


With the increasing difficulty of discovering new oil 
reserves, attention is necessarily focused on the efficient 
development and production of existing reservoirs. Effi- 
cient production depends on the availability of reliable 
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means for predicting the pressure-production behavior 
of oil fields. The prediction is usually performed by 
analytical methods or computers. Analytical methods 
are often limited to systems which can be approximated 
by simple geometry. However, analytical methods have 
the distinct advantage of being readily performed with- 
out special equipment. Analog and digital computers 
lend themselves to a variety of problems, but are not 
always locally available and consequently their use is 
limited. 

One particular situation in which analytical tech- 
niques are needed involves studies of multiple, water- 
drive oil fields sharing a common aquifer. Such oil fields 
are in hydrodynamic communication. Therefore, pro- 
duction from any of the oil fields is accompanied by 
a pressure decline which is transmitted through the 
aquifer to the other oil fields and manifested as pres- 
sure interference. The rate of propagation of the pres- 
sure decline is such that the pressure may be signif- 
icantly reduced many miles away from a producing 
pool. An example of pressure interference is found in 
the Woodbine formation in East Texas. There, a dozen 
or more pools, including several of major size, draw 
on the Woodbine formation for water drive. The pres- 
sure in the water-bearing portion of the formation has 
been significantly reduced over an area of several thou- 
sand square miles. This reduction of pressure has been 
such that when the Hawkins pool was discovered in 
1940, the pressure in it had already been reduced 280 
psi below the initial value in the aquifer. This resulted 
from prior production from other Woodbine sand pools, 
including, in particular, the East Texas field some 20 
miles to the east of Hawkins.’ 


The above example indicates that the interference 
effect can be significant and should be taken into con- 


‘References given at end of paper. 


sideration in predicting the pressure-production behavior 
of oil fields in a common aquifer. This preblem is be- 
coming more important because of the increasing fre- 
quency of its occurrence. In addition, it usually con- 
cerns more than one operator, in which case unitization 
may be deemed desirable. 

Analog computers have been used successfully to 
treat this type of problem’*’, but such computers are 
not always available to practicing reservoir engineers. 
In view of this, a practical analytical method has been 
developed and is presented in this paper. It provides 
the practicing reservoir engineer with a handy tool for 
analyzing oilfield interference problems. 


PRELIMINARY CONSIDERATIONS 


SUPERPOSITION THEOREM 


The procedure employed to treat oil fields in inter- 
ference entails the principle of superposition. Math- 
ematically, the superposition theorem‘”’ states that the 
linear combination of particular solutions to a linear 
and homogeneous differential equation is a solution to 
the differential equation. The superposition theorem is 
a useful tool for treating systems upon which involved 
boundary conditions are imposed. The general solution 
is the summation of the particular solutions obtained 
by treating one boundary condition at a time. 

The transient flow of compressible liquid in a uni- 
form porous matrix is described by the diffusivity equa- 
tion, which is linear and homogeneous and therefore is 
subject to superposition. It follows, that the performance 
of multiple oil fields in a common aquifer can be 
evaluated from the separate solutions obtained by deal- 
ing with one field at a time. In essence, if the pressure 
change associated with producing each and every field 
were computed individually (i.e. ignoring the presence 
of all oil fields but one) for a time instant, ¢, and at 
some arbitrary point in the reservoir system, then the 
total pressure change at this point and time instant is 
given by the sum of all individual changes. For that 
matter, the arbitrary point may well be the effective 
center or any other appropriate point in an oil field, 
and the time instant may represent some assigned future 
date at which it is desired to predict the behavior of 
the system. 

For illustration, assume two oil fields, A and B, 
located in a common aquifer. Producing oil field A is 
accompanied by a pressure change in A and a relatively 
smaller pressure change in B. Similarly, producing oil 
field B is accompanied by simultaneous pressure changes 
in B and A. 

Let P,, = pressure change in A due to A’s production, 
P,, = pressure change in A due to B’s production, 
and define similarly P;,, and P,,. Then the total pres- 
sure Change in A, 
Py, = Py + Pra 
and the total pressure change in B, 
Py = Poy + Pru. 
Px, Py., and similar terms will be labeled hereafter 
“interference pressure change,” as distinguished from 
P,,, Ps», and similar terms which represent the pressure 
change in an oil field due to its own production. 

The pressure change terms depend on the rates o! 
production and the physical characteristics of the res- 
ervoir system. The case of an extensive aquifer is 
treated by solving the diffusivity equation. The solu- 
tion represents the pressure distribution in the reservoi! 
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system and permits the evaluation of the individual 
changes. 


METHOD OF IMAGES 

The aquifer mentioned earlier was assumed exten- 
sive. As a modification, it is proposed to consider the 
case of an oil field located in an extensive aquifer, one 
side of which is limited by a straight extensive fault. 
In this case the benefit derived by the oil field from the 
aquifer is curtailed by the presence of the fault, and 
the resulting effect can be expressed as pressure inter- 
ference. This situation is usually referred to as pro- 
duction near a fault* and is treated as a special case of 
oilfield interference by applying the method of images.’ 

A fault is defined here as a boundary perpendicular 
to which no flow will take place. This fact can be rep- 
resented analytically by placing a mirror image of the 
oil field in offset position and assuming the image to 
be produced in a manner identical to the oil field. The 
mirror image represents a fictitious oil field which is 
placed symmetrically with respect to the fault. This 
procedure reduces the component of flow perpendicular 
to the line of symmetry to zero and therefore fulfills 
the boundary condition required by the fault. The prob- 
lem is thus transformed to a simple case of interference 
between two identical oil fields located in an extensive 
aquifer. The pressure change observed at the oil field 
will be the sum of the effects of the oil field and its 
image. Evidently, the effect of the image will be small 
in the early life of the oil field and will grow larger 
with time, depending on the distance between the oil 
field and the fault 


ASSUMPTIONS AND LIMITATIONS 


Before proceeding to the solutions of the diffusivity 
equation, it would be advisable to review some charac- 
teristics of reservoir systems to be treated in the fol- 
lowing sections of this paper. The reservoir system is 
assumed to consist of an extensive aquifer in which a 
number of oil fields are located. It is subject to the 
following definitions, assumptions, and limitations. 


AQUIFER CHARACTERISTICS 


1. An extensive aquifer is defined as an aquifer, the 
exterior boundary of which is sufficiently removed from 
the oil fields that it does not influence their perform- 
ance during the period of analysis. For all practical 
purposes, this assumption is valid for an aquifer having 
a size about 600 times that of the largest oil field* in 
the reservoir system. 

2. In applying the method of images to treat the 
problem of production near a fault, the aquifer is as- 
sumed extensive and bounded by a straight fault. This 
aquifer is represented geometrically by an infinite half 
plane and the assumption is justified on conditivn that 
the size of the aquifer is about 300 times that of the 
oil field and the length of the fault is about 25 times 
the field’s effective diameter 

3. The compressibility coefficient, which is defined 
as the sum of the water and the rock effective com- 
pressibilities, is assumed constant and independent of 
the pressure. Note that while in the following analysis 
only the expansion due to the water compressibility is 
taken into account explicitly, superposed effects due to 
the compressibility of the rock, which may be appre- 
ciable, readily fall within the scope of the analysis. For 
it may be shown’ that if the porous medium, as well as 
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the entrained fluid, be considered as compressible, the 
diffusivity equation remains unchanged, except that the 
compressibility coefficient is to be considered as the 
sum outlined above. 


4. The porous matrix is isotropic and uniform with 
respect to permeability, porosity, and thickness. 


5. All flow through the aquifer is macroscopically 
laminar, two-dimentional, and in a horizontal plane. 
Thus, the flow is governed by Darcy’s law and the 
effects of gravity are neglected. 


OILFIELD CHARACTERISTICS 


1. Each oil field is considered as a “ring” of produc- 
tion with influx equal to the water influx into the field 
and with pressure equal to the pressure at the original 
oil-water interface. A “ring” of production is used to 
imply a circular sink, as distinguished from a point 
sink, and represents an oil field the areal configuration 
of which can be approximated by a circle. Note that 
radial flow obtains when the effect of each field is con- 
sidered individually. 


2. The rate of water influx is equal to the production 
rate, provided that the expansion of the oilfield fluids 
and rock resulting from pressure changes is negligible. 
If this assumption were not justified, reliable predictions 
could be achieved by employing an iterative scheme as 
shown in a following section. 

3. The interference pressure drop computed at the 
effective center of an oil field represents the areal aver- 
age pressure drop, if the oilfield radius were small, 
and/or the distance between the oil fields were large. 
If neither condition were satisfied, the interference pres- 
sure drop should be calculated for more than one loca- 
tion in the oil field. The number of locations depends 
on the desired degree of accuracy. 


The above assumptions should not imply that actual 
oil fields are strictly circular and that aquifers possess 
isotropy and uniformity. However, these assumptions 
provide a convenient approximation for many actual 
systems. The ensuing results, as evidenced by past 
experience, will give good estimates even when applied 
to systems that deviate considerably. 


SOLUTIONS OF DIFFUSIVITY EQUATION 


The individual components required to evaluate the 
total pressure change are obtained from solutions of 
the diffusivity equation considering one ring of produc- 
tion. The diffusivity equation describes non-steady state 
flow of a compressible liquid in uniform porous media 
and is derived” by combining Darcy’s law with the 
continuity equation. By dealing with one ring of pro- 
duction at a time, radial flow obtains and the diffusivity 
equation takes the form 

ap + ee fuCw ep (1) 

: ar | or Or 2. 
where f, k, » and C,. represent porosity, permeability, 
water viscosity and compressibility, respectively. The 
term Loan represents the diffusivity of the aquifer and 
is advantageously eliminated by introducing dimension- 
less time, ft, pressure change, P(r, ¢) and radius, r, such 





*See Letter Symbols for Reservoir Engineering, J. Pet. Tech., 
Jan., 1955, 38. Symbols not included in above mentioned list are 
explained in Nomenclature. 
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that, 
ieee ae ee eee. (2) 
P(r, t) = 2aRAPO, O/mg. sw ee. 3), 
and r=1/te. Britten. sae : wee, (4), 


where r, is the radius of the ring Ps production and 
represents the effective radius of the oil field. Equation | 
then becomes 

oP 1 oP 

or r or ot ; 

The pressure change in an oil field due to its own 
production (P,,, Pr, ... etc.) is obtained by solving 
Equation 5 for r = 1 and has been treated by Van 
Everdingen and Hurst." The interference pressure change 
(Px, Peo, ... etc.) in an oil field due to producing an 
adjacent oil field is obtained by solving Equation 5 fo 
r equal to the dimensionless distance between the cen 
ters of oil fields. 


ed 


~ (5 





The solution of Equation 5, for all practical values 
of r and ¢, is developed in this paper for the constant 
rate case and the following boundary conditions 
(1) P(r, 0) = O indicating uniform initial pressure 
over the entire aquifer; (2) P(r, t) > Oasr— ox, 


ere ; ; oP 
indicating an extensive aquifer; and (3)/ 3 ) 
r 


indicating a constant volume rate of water influx of 
27k/» per unit thickness. 


This solution represents the pressure distribution 
P(r, t), in an extensive aquifer due to a ring of con 
stant production. Once P(r, t) is obtained for a con 
stant rate of production, it can be easily computed for 
any rate which may vary with time. Again, this is 
accomplished by the superposition theorem. It will be 
demonstrated in the following section of this paper 

The above-mentioned solution is treated explicitly in 
the Appendix. The solution is obtained partially from 
mathematical considerations employing the Laplace 
transform and the concept of a continuous point sink 
and partially from analog and digital computations. Th« 
former provide accurate expressions of P(r, ¢) for 10 

t > 500. However, in treating oil fields (in which 
case r. can have a large numerical value) the value of 
t can be easily less than 500 even for large values of 
absolute time, t. Therefore, analog and digital com 
puters were used to bridge the gap between f 01 
and 500. 


Analog computation was performed by the use of a 
resistance-capacitance type analyzer. Digital computa- 
tion was accomplished by translating the diffusivity 
equation into a set of difference equations which wer« 
solved simultaneously with the aid of a digital comput 
ing machine called “Datatron.” The digital computation 
provided values of P(r, t) which are accurate to thre« 
decimal places. This is adequate accuracy considering 
the degree of uncertainty associated with most prop 
erties of reservoir systems. 

The solution of the diffusivity equation for the con 
stant pressure case was treated by Andrew Gemani 
for the following boundary conditions: (1) P(r,0) = 0 
indicating uniform initial pressure over the entir 
aquifer; (2) P(r, t) > 0 asr > o, i mee an eX 
tensive aquifer; and (3) PCI, , indicating con 
stant unit pressure drop at r = 1. 


The rate of water influx asociated with this set of 





boundary conditions was discussed in detail by Carslaw 
and Jaeger’, and Van Everdingen and Hurst.” 


PRACTICAL CONSIDERATIONS 


Values of the dimensionless pressure change, P(r, ¢), 
for unit rate of water influx are presented in graphical 
form. Figs. | to 3 cover the ranges of ¢t from 
10° to 10° and r from 1 to 64. These include essen- 
tially all ranges of r and ¢ encountered in practice. Val- 
ues of P(r, t) for 10° > t > 10° and r > 6% can be 
accurately calculated from Equation 9 of the Appendix. 
The charts are readily used to evaluate the pressure drop 
at any point in an extensive aquifer due to some pre- 
scribed rate of vil production.* Presenting the charts 
in dimensionless form permits the evaluation of the 
pressure drop for any reservoir system independent of 
its physical characteristics, but within the limitations 
outlined earlier. 

If at zero-time the rate of production of an oil field 
A, of radius r.,, were some value q, then the pressure 
drop at its boundary over a period of time f is given 
by the relation, 


*For water-drive reservoirs the rate of oil production is equal to 
the rate of water influx plus the rate of expansion of hydrocarbon 
content. The rates of oil production and water influx will be used 
interchangeably hereafter, since one is easily converted to the other 


1x to 


P,,(t) = P(re., t) = 887.6 —— (q,) P(1, 10), 

~ l 
and the interference pressure drop in an adjacent oil 
field B is given by, 


P,.(t) = P(r. re. t) = 887.6 (gs) Pra; t). 
1 


2k} 
Where ¢ is the dimensionless equivalent of 1, 
kt 
fuC,a’ 
and r, is the distance between the centers of A and B 
in multiples of r.,. The variables in these and the fol- 
lowing equations are expressed in units common to oil- 
field usage. Thus, 
P = pressure change in psi, 
ft = viscosity in cp, 
C,.** = water compressibility vol/vol/psi, 
q = production rate in reservoir B/D, 
¢ = permeability in md, 
= thickness in ft, 
time in days, 
porosity, fraction, and 
oilfield area in acres 
Furthermore, if during the period of interest, f,, the 
production rate of A varies, then the pressure drop in 


t = (4.56) (10°) 


‘*In case the rock compressil 
sent the sum of the water and 


ippreciable, Cw should repre- 
ve compressibilities. 


ix 10 1x10 
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psi at any point in the aquifer is given approximately 
by the sum 

P(r, t.) = 

a= 7% 


8 : LL | a (b4+3) —, q(t) P(r, ) aw t,) . (6) 
— a - 
1 = 


where ¢, represents time instants (/,, f,, 2... ft) and 
qa(t,) represents average production rate during the 
time interval f,-, to f,. 


Inserting the proper values of r and the corresponding 
values of P(r, t) in Equation 6, gives the pressure drop 
in A and the interference pressure drop in B due to 
production from A. The same approach is employed to 
yield the effects of production from B, thus obtaining 
the total pressure drop in A and that in B. This proce- 
dure is used to predict the pressure behavior of two 
oil fields in a common aquifer for preassigned produc- 
tion schedules. 


Let it be required to predict the pressure-production 
behavior of oil fields A and B if A’s production and 
B’s pressure were preassigned and denoted by q, (1) 
and P, (t), respectively. The procedure in this case 
consists of evaluating the pressure drop in A and that 
in B due to q, (t) by using Equation 6 as pointed out 





duction, is determined next and equals the preassigned 
pressure drop, P;,(t), less the pressure drop at B due to 
q,(t). Once Py,(t) is determined as a function of time, 
the cumulative production from B is predicted from 
the following equation: 


Qs(t.) ” 


ng -- 5 

15516 fC,Ah S [Patt 4 Pati) O(t, t,) 
i=0 
iz 


° ° . ° ( } 
The variables in this equation are expressed in units 
common to oilfield usage and Q is dimensionless wate! 
influx term for unit pressure drop at r = |. The appli 
cation of Equation 7 was illustrated by A. T. Chatas 
with the aid of Van Everdingen and Hurst’s” constant 
pressure solution. The rate of production from B, g,(1), 
it then used to evaluate P,,(t) by employing Equation 
6. The sum of P,,(t) and P,,(t) gives the total pressure 
drop in A. 

Another combination arises if the pressure drop in 
both A and B were preassigned and denoted by P,(1) 
and P,(t) respectively. To predict the rate of produc 
tion of each field, it is necessary to resolve the total 
pressure drop in each into the components, P,,(1), 
Px,(t), Pr» (t) and P;,(t), which are related by the fol 
lowing four algebraic equations: 

Pelt) se Faghte), 4 Besbta)iencvs- - par) 








Ai. 
7] 


















































earlier. The pressure drop at B, P,,, (;), due to its pro- Pate) = Pate) + Pats). . - (9) 
' 1x10 1xto® 1x10 
1.2 1.2 
= 
it 
} 
z / / 
1.0 7 
| 
NOTE / / : 
- SEE FIGURE | FOR / 
a= DEFINITION OF TERMS 
Li 8 8 
a.| 
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P(t.) = 


Pati) — Pelt.) [PCr !, — t\) . 





Pos(ty) = 
—- j 


n 
; [Patt ~ Pats) PCr — th). . ° ° (11) 


i=0 

P (1, t) is dimensionless pressure drop in the aquifer 
resulting from unit pressure drop in oil field, and can 
be obtained from Andrew Gemant’s” constant pressure 
solution. Inserting n = 1 permits the evaluation of 
P,,(t,) and P»,(t,). This step entails solving four alge- 
braic equations in four unknowns. Similarly inserting 
n = 2, 3,4... gives values of P,,(t) and P,,(t) for 
all time instants. The production rate of oil field A is 
then predicted from P,,(t) and Equation 7. Similarly 
P(t) is used to predict the production rate of B. 

The above analyses are readily extended to cover 
multiple oil fields in a common aquifer, and are sub- 
ject to the following considerations: 

1. It was tacitly implied that the rate of water influx 
and the rate of production in reservoir barrels are equal. 
In case this assumption is not justified, more accurate 
predictions can be obtained by employing an iterative 
scheme. The procedure outlined earlier is used first to 


1x10 


yield preliminary predictions. Then, the rate of water- 
influx is calculated from the prekminary predicted pres- 
sure and material balance equations, taking into account 
the compressibility of the oilfield contents.* The water- 
influx rates thus obtained replace the production rates 
in the above analysis to secure more accurate predic- 
tions. 

2. The rate of production and the pressure drop may 
be represented by a stepwise variation as a function of 
time, provided that the step intervals are not taken to 
be too large. 


[ILLUSTRATIVE EXAMPLE 


In the present paper no attempt has been made to 
apply this method on an actual field case. However, 
an illustrative example is presented to demonstrate the 
methods and procedures to be followed in the solution 
of typical interference problems. It is planned to apply 
this method to actual field cases; these results will be 
presented in a later pape! 

Consider the case of a reservoir system comprised 
of two oil fields, A and B, which are located in an 
extensive aquifer and the distance between their cen- 
ters is 7.98 miles. A and B enclose 2,000 and 8,000 
acres of surface area, respectively, and are produced 
under complete water drive. The pressure-production 


This includes oil, gas, connat« te r ock compressibility 
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history of these fields, given in three-month increments, 
is shown in Figs. 4 and 5. Note that field A has been 
produced for the past four years and field B for the 
past two. The original pressure in the reservoir system 
was 4,000 psig, while the current pressure is 3,219 
and 3,220 psig in A and B, respectively. The system 
is further characterized by the following physical prop- 
erties; 

k = 212 md, f = 0.22 

h = 84 ft, B 1.15 bbI/STB 

C, = 5.0 X 10° per psi, = 0.30 cp 

Let it be required to predict the pressure in both 
fields if A were produced with the present rate of 75,000 
and B at 120,000 STB/D. 

If the contribution of the oil compressibility to the 
energy of the system is neglected, and if it is assumed 
that every reservoir barrel of crude oil removed from 
the oil fields is replaced by a reservoir barrel of water 
from the aquifer, then the problem may be solved 
directly by utilizing Equation 6. 

Since p, B, k, h, and q are known quantities, only 
the dimensionless pressure change terms must be deter- 
mined to solve the problem. This is done by first cal- 
culating the appropriate dimensionless time and radius 
and then locating the corresponding dimensionless pres- 
sure-change in the charts. Note that the numerical val- 
ues of the time, f, and consequently the dimensionless 
time, t. are always referred to the discovery date of the 
oil field under examination. The predicted performance 
is shown in Figs. 4 and 5. 

As a sample calculation of the former step the pres- 
sures eight years hence are determined as follows: 
Field A: 

t, = 12 years = 4,383 days 

kt, 

t, = (4.56) (10° Fa C.A 

(4.56) (10°) (212) (4,382) 


(.22) (.30) (5X10) (2,000) 


= 642.0 
r = 7.98 miles 
2,000 
lwo = eee .997 miles 
640 x 
7.98 g 
mn, = oy 2 
; .997 
TABLE 1 — CALCULATION OF PRESSURE CHANGES IN OILFIELDS A AND B 
DUE TO PRODUCTION FROM OILFIELD A 


P 642 — t.) P(8,642 t) 
642 — +t BS 
dimensionless 





t 
t O8.a0 A a 


days STB/D dimensionless time Pressure Change 
0 14,000 0.0 642.0 3.64 1.57 
90 10,000 13.2 628.8 3.62 1.56 
181 5,500 26.5 615.5 3.61 1.55 
273 10,500 40.0 602.0 3.60 1.54 
365 12,000 53.5 588.5 3.59 1.53 
455 — ,500 66.6 575.4 3.58 1.52 
546 15,500 80.0 562.0 3.57 1.51 
638 4,000 93.4 548.6 3.56 1.50 
730 4,000 106.9 535.1 3.55 1.49 

4383 642.0 


Substituting these values in = following equations, 


(0.30) (1. L9> | 
> 87.6 
P.., (887. ) 5 (312) (84) Galtier) q(t) 
r= 0 


P( 1,642-t;) 





VOL. 


1955 





204, 








8 

(0.30) aoe ag ve 

(887.6) 2m (212) (84) z laste galt; | 
i= 


P(8,642-1,) 
Yields P,, = 739.5 psi 


Py, = 315.1 psi 
Field B: 


II 


t, = 10 years = 3,653 days 


kt, 
t, = (4.56) (10%) juCe _ 
(4.56) (10°) (212) (3,653) 
(.22) (.30) (5X10) (8,000) 
= 133.8 
r = 7.98 miles 
fun = eed = 1.994 miles 
7.98 


ok: salle 





TABLE 2 — CALCULATION OF PRESSURE CHANGE IN OJLFIELDS B AND 4 
DUE TO PRODUCTION FROM OILFIELD B 


P(1,133.8 — ¢ #,) Pld, 133.8 








t 1 i Sgn 
Te ERE cod Seen 
days __$TB/D dimensionless time _ Pressure Change 
0 20,000 0.0 Sa So 1.49 
90 18,000 3.3 130.5 2.85 1.48 
181 14,000 6.6 127.2 2.83 1.46 
273 19,500 10.0 123.8 2.82 1.45 
365 18,500 13.4 1204 2.81 1.44 
456 — 6,000 16.7 117.1 2.80 1.42 
547 27,000 20.0 113.8 2.78 1.41 
639 8,000 23.4 110.4 2.77 1.40 
731 1,000 26.8 107.0 2.75 1.39 
3653 133.8 


Substituting these values in “¥§ following equations, 


(0.30) (1.15) a 
Pr. 887.6 B i*1 B 
(887.6) 5 (212) ZY Ga(tins) — Gn(ts) 


P(1,1335:6:— 8) 


8 
(0.30) (1. 2) 
Ps 887.6 ? 
( ) Se | (212) (84) = lan 1) Qn t, | 


P(4,133.8 — t,) 


Yields Py», = 925.7 psi 
Px, = 475.4 psi 
Thus, total pressure drop in A = P,, + Py, = 1,214.9 
psi, and 
total pressure drop in B = Py, + Py. = 
psi 


_— 


240.8 


NOMENCLATURE 


E, = Exponential integral. 
Modified Bessel function of the second kind 
and nth order. 


P,, = Pressure drop in oilfield A due to its produc 
tion. 
P,, = Interference pressure drop at oilfield A due 
to producing oilfield B. 
P, = Total pressure drop in oilfield A. 
P(...) = Pressure drop at radius, r, and time, /. 
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= Dimensionless pressure drop at dimensionless 
time ¢ and radius r. 

= Production rate of oilfield A at time t. 

= Radius of oilfield A. 


. : : ° r 
= Dimensionless radial distance, equals 
. r 


= Dimensionless distance between centers of oil- 
fields A and B in multiples of r,... 

= Dimensionless time, equals kt/fuC.r’.. 

= Dimensionless potential drop. 

= Reduced pressure, equals p/p,. 


APPENDIX 


LAPLACE TRANSFORM AND CONTINUOUS 
POINT SINK 


Laplace transform of P(r, t) is defined as, 
a0 


- 


ol oe BO ee! eo 


~ 


0 
When applied to Equation 5 of the text for the constant 
rate case, it gives 
K.(rvV's) 
P(r, s) =) eagles RD 
s” K,(Vs) 
Where K, and K, are modified Bessel functions of the 
second kind. For small and large t, Equation 2 is sim- 
plified and the ensuing inverse transforms are readily 
evaluated. 
If ¢ is small, the operator, s, is large and K,,(z) is 
reduced to”, 


a- 


~~ 


K,,(z) = q . e 


substituting for K, and K, in Equation 2, it becomes 


” ; 
Paley Pe pete) sl ans at a~..° (3) 


Vrs 
The inverse Laplace transform of Equation 3 is’, 
2ViIe r— | 


2\ 
P(r, t) = ierfc 
Vr 2V iI 


where ierfc x = 


; xerfe x 
Vr 


oo 


* ; 
and erfc x = af e Fag 
/ a 


x 
ierfe x is readily evaluated from “Tables of the error 
function and its derivative’”. The relation given by 
Equation 4 is accurate for r < .01. 
If ¢ is large, the operator s is small and K,,(z) is 
given by 


> 


K,(z) (1m -~+y},K gay 
where y “Euler’s constant” = 0.57722. Again substitut- 
ing for K, and K, in Equation 2, it reduces to 
Ins In*/r y 
P(r, s) t celeitinic: aivte Cee 


2s 5 


The inverse Laplace transform of Equation 6 is given 
by Carslaw and Jaeger’, 

l 4t ] 

~|In } oS ce 


r 


Pir, t) 


t 
This relation is accurate for == 2,000 


Equation 7 is also a solution of the continuous point 
sink for large ¢. The continuous point sink solution, 
originally developed by Lord Kelvin, is expressed as 


] . 
P(r, t) » eI 


where — E,( 


4t 


The “E,” function is known as the exponential integral 
and its values are given in “Tables of sine, cosine and 


4t 


exponential integral”. For 2,000 the exponential 
r 


4t , 
integral reduces to /n y, and the pressure gradient 


r 


oP 
becomes 
or r 2. ¢ S 30D 


Consequently, only for ¢ 500 the continuous point 
sink solution describes accurately the solution to a 
“ring” of production. 

It is evident that the continuous point sink solution 
does not apply for the determination of the pressure 
change when f is small. However, when the radius, r., 
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is small, such as a well radius, even small values of the 
absolute time, f, will give large values of the dimension- 
less time, ft, and the continuous point sink is applicable. 
On the other hand, in treating oil fields (in which case 
ry can have a large numerical value) the value of ¢ 
can be easily less than 500 even for large values of 
absolute time, f. 


In recapitulating, 


= J 
= 2V1 ote! man sort < 01 


Vr 2V1I 


1 ; wee 4 for t > 500 
-3[-«, a a 


Therefore, to bridge the gap between 1 Ol and 500, 
analog and digital computers are used. 


P(r, t) (9) 


ANALOG COMPUTATION 


Analog computation based on the fundamental equa- 
tion for flow of electric current in a dielectric medium’, 
was performed by the use of a resistance-capacitance 
type analyzer. The reservoir system under consideration 
was divided into segments by concentric circles and 
each segment was represented by a condenser and a 
resistance. Constant current, i, was drawn from the 
circuit and the voltage was recorded as a function of 
time for every segment. A sufficient number of segments 
was included in the circuit such that the last condenser 
showed no voltage drop during the time interval of 
interest, thus simulating an extensive aquifer. Because 
of practical considerations, 10 was the maximum num- 
ber of condensers used. The results obtained from the 
analog computer are approximate due to truncation 
error introduced by the limited number of segments 
available in the analog computer used, and the limita- 
tions imposed by the available recording instruments. 


DIGITAL COMPUTATION 


Digital computation was employed to achieve better 
accuracy. The diffusivity equation was translated to a 
set of difference equations which were solved simulta- 
neously with the aid of a digital computing machine 
“Datatron.” Difference equations: 


- 


Let u = Inr, and ¢ = (10) 


where ¢ is reduced pressure and p, is initial pressure 
in the aquifer. Substituting Equation 10 in Equation 5 
of the text, it becomes 
e° aa 
oe = on @ . 
ou Ct 
The step-by-step calculations of the pressure change 
were performed on the finite-difference form of this 
equation. Dividing the aquifer by a number, n, of con- 
centric circles and denoting the time stations by m, 
Equation 11 takes the following finite-difference form”” 
54 Pn-1 rid, 2 Pu + Pats + Pn a 2 dn t Puts ! 
(Au)* (Au)* ant 
sua ED m*1 ~ Bi, me} 
| At 
By rearranging ‘the terms and choosing equal Au incre- 
ments, the above equation can be rewritten as 


(11) 





; bus 2 ! “}- (eur ee da + dn — ge (12) 
= fp _ - (Au)? 2n \u F { 
) gait 2]1 At e |. Pats a 
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From the condition of constant rate and Equatior 10, 
it follows that 


0 : 
(:*) is Ne 13) 
Ou ao 


and the rate of water influx becomes 


2akh 
q=——p.H 
b 


The pressure distribution associated with this rate 1s 
converted directly to the corresponding pressure change, 
P(r, t), thus 


ae 
P(r, t) Sil diabglilies allen eet (14) 
Forn = 1 
da-1.m = Pom = Pw.m 
where ¢, is reduced field pressure. From Equation | 4 
and Taylor’s expansion we get, 
Pim = Gw.m + (Au)H + 


2 2 
(G8) eee 2 ee, (ae oe 6. | C15) 
2! Ou? J uo 


Expressing( 7) in difference form, Equation 15 


becomes 


_4 Pe Reg 
Dwim 2° 3 Pi, m 3 od» m cy Au 


By substituting the above relation in Equation 12 for 
n = 1, it takes the form 


a 3 Aur <j ee phates Phos, 116) 


Phat so v 
ib Bie 82 |e by 


The computations are based on Equations 12 and 16 
A technique was followed by which 30 difference 
equations were set to describe a limited radial system 
The pressure was computed for consecutive time sta 
tions ¢,, ft, ft, ... Up to f, at which the outer boundary 
showed the first sign of pressure decline. Pressures cor 
responding to f,., were used as boundary values for a 
larger radial system. By using this procedure it was 
possible to take smali radial increments with a limited 
number of difference-equations. Values of A? were 
chosen progressively smaller till the computed values of 
P(r, t) showed no change in the third decimal place 
The values of P(r, ft) thus obtained are accurate to 
three decimal places. This is adequate accuracy con 
sidering the degree of uncertainty associated with most 

properties of reservoir systems. 

It is self evident that Equations 12 and 16 are suit 
able to secure the pressure distribution associated with 
a limited radial system for all values of dimensionless 
time, ft. One might also add that the solution for the 


constant pressure case is readily computed from Equa 
tion 12 and an auxiliary equation stipulating constant 
pressure in the oil field. The auxiliary equation is 
given by 


,- 142 
\ 


1-4 ole. ~ if 

| ] f 

The case where the permeability varies with radia! 
distance is being examined by using a modified form of 
Equation 12 to account for the permeability variation 

It might be well to mention that several attempts 
were made before the solution was perfected to th« 


wes gin oh Seas (17 


met 


Pe 


desired degree of accuracy. Preliminary computations 
were performed neglecting the third term in the right 
hand side of Equation 15, and led to values of P(r, ¢) 
which included considerable error especially for small ¢. 


CONCLUSIONS 


Non-steady state pressure distribution in an extensive 
aquifer due to constant rate of production from a cir- 
cular sink has been developed partially from math- 
ematical considerations and partially from analog and 
digital computations. By virtue of the superposition 
theorem, the pressure distribution provides a practical 
analytical method for evaluating the interference effect 
in water-drive oil fields sharing a common extensive 
aquifer. 

The case of production near a fault is treated by 
applying the method of images. In addition, provisions 
are made for treating oil fields, which are produced 
by partial water-drive. 
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ABSTRACT 


A series of water floods was made 
on laboratory prepared unconsoli- 
dated sand colun:ns to study the ef- 
fects on oil recovery of the solid- 
water-oil contact angle, the oil-water 
interfacial tension, flood rate, and 
oil viscosity. 

A procedure was developed for 
treating silica sand with silicone pol- 
ymer to produce surfaces of varying 
degrees of wettability as evaluated 
by measuring the contact angle on 
a flat silica plate that had been simul- 
taneously treated with the sand. Fluid 
interfacial tensions were changed by 
the use of surface active agents in 
the flood water. Water-wet, oil-wet, 
and intermediate wettability systems 
were studied at high and at low val- 
ues of interfacial tension, at several 
flood rates, and with oils of different 
viscosities. 

For both oil-wet and water-wet 
systems and a low viscosity oil, re- 
coveries were functions of the oil- 
water interfacial tension; also, in- 
crease in flood rate resulted in in- 
creased oil recoveries. High inter- 
facial tension floods were more effi- 
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cient than low interfacial tension 
floods on water-wet systems, while 
low interfacial tension floods were 
more efficient on oil-wet systems. In- 
termediate or neutral wettability sys- 
tems were less sensitive to rate of 
flood advance and interfacial tension 
than either oil-wet or water-wet sys- 
tems. 

The effects of surface forces on oil 
recovery for high viscosity oils were 
not so well defined as for low vis- 
cosity oils. 


INTRODUCTION 


Reservoir rock surfaces vary in 
their wettability’, some being water- 
wet while others are apparently oil- 
wet. The factor of wettability enters 
into reservoir performance in both 
primary and secondary recovery, and 
into laboratory measurements of cap- 
illary pressure, connate water satura- 
tion, and relative permeabilities.’ In 
theory, the degree of wettability is 
measured by the contact angle for 
the system solid-water-oil. A water- 
wet system is defined as one in which 
the advancing water contact angle is 
less than 90°. An oil-wet system is 
one in which the advancing water 
contact angle is greater than 90°.’ 


‘References given at end of paper. 
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The surface forces in a solid-wate: 
oil system are related by the Young 
Dupree equation as follows’: 

hase ~ Big ™ See — Cae 

= g.~oos8 = —AF... (1) 

These surface forces are common 
ly expressed as the difference in ad 
hesion tension, or the product of the 
water-oil interfacial tension and the 
contact angle. If the contact angle is 
less than 90°, the adhesion tension 
difference is positive and oil displace 
ment by water from a solid surfac« 
will be spontaneous and favored by 
a high oil-water interfacial tension 
If the angle is greater than 90° the 
displacement will not be spontaneous 
but will require less energy the lower 
the oil-water interfacial tension. The 
contact angle and interfacial tension 
also appear in the capillary pressur: 
equation’: 

i ee CO a 5 (2) 

The quantity ocosO was used by 
Rapoport and Leas® as a factor in 
the generalized form of their “sca! 
ing coefficient.” Warren and Cal 
houn’ presented the results of «cos 
on oil recovery in oil-wet systems 
and correlated their data using a 
“scaling coefficient” suggested by 
Jones-Parra and Calhoun.” The “sca! 
ing coefficients” were advanced t 
help correlate the effects of the va 
riables influencing oil recovery by 


external drive and specifically to 
clarify the effect of injection rate 
upon oil recovery in laboratory flood- 
ing tests. The effect of increased in- 
jection rate has been reported by va- 
rious authors to increase, to decrease, 
and not to affect oil recovery. Levine’ 
as well as Rapoport and Leas’ cite 
several references as to the effect of 
injection rate on recovery. 

It was the object of this work to 
investigate separately the effects of 
contact angle and interfacial tension 
as measured external to the system 
on oil recovery in water flooding of 
unconsolidated sand columns at va- 
rious injection rates and at different 
oil-water viscosity ratios. The contact 
angle was controlled by treating the 
sand with silicone to render the sand 
surface to varying degrees of wet- 
tability. Interfacial tension was al- 
tered by the addition of surface ac- 
tive agents to the flood water. 


EXPERIMENTAL PROCEDURES 


MATERIALS 

The flow models were made of 
Lucite tubes 59 cm long, 2.52 cm in 
diameter, and 295 ml capacity. The 
tubes were packed with 80-120 mesh 
silica sand to an average permeabil- 
ity of 11 darcies to oil and average 
porosity of 33.5 per cent. The sur- 
face area of the sand was calculated 
from the porosity and permeability 
by the Kozeny equation” as 452 cm’*/ 
gem of sand. 

Oils of three different viscosities 
were used. A decane-dodecane par- 
affinic hydrocarbon, Soltrol C, was 
used in a majority of the experi- 
ments. Two oils of 32 cp viscosity 
but of different interfacial tensions 
against water, and an oil of 265 cp 
viscosity, were also used. Properties 
of the oils are listed in Table 1. 

Tap water of 0.895 cp viscosity 
at the temperature of the experi- 
ments, 25°C, or tap water contain- 
ing dissolved surfactant, were the 
displacing phases. The surfactants 
were Igepal CO-630 and Igepal CO- 
880 from Antara Chemicals, and 
Pluronic F-68 from Wyandotte 
Chemicals Corp. These are non-ionic 
type surfactants but with differing 
degrees of interfacial activity against 
water. 

The silicone used to treat the sand 
surfaces was Dow-Corning DC-1107, 
a 20-40 centistoke viscosity methyl- 
siloxane polymer. 


PROCEDURE 


The sand was treated with a 
known weight of silicone dissolved 
in isoproyl alcohol. Enough alcohcl 
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TABLE 1—PROPERTIES OF HYDROCARBONS USED IN WATERFLOOD EXPERIMENTS 
interfacial Tension 
Specific Viscosity, 25°C Vs. Water, 25°C 
Gravity, 25°C __Centipoises Dynes/cm 
0.754 1.46 = 
0.852 35.8 16.8 
0.836 35.8 34.5 
0.873 296 35.8 


Oil-Water 


Hydrocarbon Viscosity Ratio 


Soltrol C 
32 CP Neutral Oi! 
32 Cp Oil Blend 
265 Cp Neutral Oj! 


flooded horizontally at constant rates 
using a quadruplex Milton Roy Mini- 
pump. Fresh sand was used to pre- 
pare a new model for each flood. Oil 
production was observed visually at 
water breakthrough and periodically 
thereafter until eight pore volumes 
of water had been injected in the 
case of the Soltrol C systems and 15 
pore volumes in the cases of the 
higher viscosity oils. 


was used to give a thick slurry, about 
200 ml per 1,000 gm of sand. The 
sand slurry was mixed thoroughly 
and after standing overnight, the 
alcohol was removed by vacuum 
evaporation with frequent mixing. 
After removal of the alcohoi the 
silicone was set by heating at 150°C 
for two hours. A flat silica plate was 
processed along with the sand. A 
measure of the oi! wettability of the 
sand surface was then available by 
measuring the contact angle on the 
silica plate. 

The silica plate was placed in a 
rectangular optical cell and allowed 
to stand overnight under oil. The 
advancing silica-water-oil contact an- 
gle was formed by placing small 
drops, abovt 2 to 4 mm in diameter, 
of the ¢’,*'acing phase on the silica 
plate. Aa umage of the solid and drop 
was shown on a screen by projecting Table 2 give the average value of 
light through the cell and using a4 the contact angle and deviation from 
by 5 in. Graphic camera with the the average as a result of several 
back removed as a projector. The treatments on each system. Precision 
angle between the solid and the of measurement of the contact angle 
water-oil interface measured through was within 5°, although in some 
the water was obtained with a pro- cases the actual deviation from the 
tractor from the image. Angles were — contact angle may have been as high 
measured at two to five minutes, at as 15°, due possibly to slight differ- 
30 minutes, at one hour, and at 24 ences in treatment. The system re- 
hours after the formation of the ferred to as water-wet consisted of 
drop. untreated sand, the oil-wet system 

Interfacial tensions were measured was highly silicone treated, whereas 
with a DuNouy vernier tensiometer. the neutral or intermediate system 
The ring was cleaned by flaming and was given an approximately mono- 
placed in contact with the aqueous molecular coating of silicone. The 
phase prior to each measurement. quantity of silicone required for a 
The ring was pulled from the aque- monomolecular coating of the sand 
ous phase into the oil through inter- was calculated to be 0.0266 mg/gm 
faces at 15 to 30 minutes’ age. Read- of sand. This value is based on the 
ings were corrected from a curve ob- sand surface area calculated from 
tained by running the surface tension _porosity-permeability data and 1.7 
of liquids of known value on the sq m as the area occupied by | mg 
instrument. Interfacial tension as well _of silicone polymer.” Contact angles 
as contact angle determinations were ranged from 45 to 160°, as meas- 
made at approximately 25°C. ured through the water. 

Prior to the flooding experiments, Contact angle data for silica, sur- 
the Lucite tubes were packed with  factant solution, and Soltrol C are 
sand in a standard manner, evacu- listed in Table 3. 
ated, filled 100 ner cent with oil and The value of interfacial tensions 
allowed to stand for 18 to 24 hours. are listed with the individual runs 
The oil-filled columns were then in Figs. 1 through 7 for Soltrol C 


EXPERIMENTAL RESULTS 
AND DISCUSSION 

SURFACE PROPERTIES 

The wettability of the systems 
studied is designated by the value of 
the contact angle. The variation in 
contact angle with silicone treatment 
for the silica, water, Soltrol C sys- 
tem is given in Table 2. The data in 


TABLE 2—SILICA-WATER-SOLTROL C CONTACT ANGLES 
Average Contact Angles, Degrees 

2-5 Min. 30 Min 

45715 30 


62+10 60 
85+10 82 


Mg Silicone No. Silicone 


Per Gm Sand 24 Hrs. 
. 25+10 
46+ 6 
71+ 6 


Layers 


0 

0.0125 
0.025 ‘ 
0.05 . 108+10 100 90+ 5 


0.10 : 140+10 127 115415 
1.0 o aw 5 157 157 154> 5 
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TABLE 3—SILICA 





Mg Silicone Surfactant Solution 

Per Gm Sand Conc. PPM 
0 500 CO-630 
0 1000 CO-630 
0 1000 CO-880 
0 1000 F-68 
0.0125 500 CO-630 
0.025 500 CO-630 
0.025 1000 CO-630 
0.050 500 CO-630 
0.10 500 CO-630 
1.0 500 CO-630 
1.0 1000 CO-630 
1.0 500 CO-880 
1.0 1000 F-68 


and in Table 4 for high viscosity 
ous. Interfacial tensions varied from 
44 to 2.8 dynes/cm in the Soltrol C 
floods and from 35.8 to 1.2;dynes/cm 
for the high viscosity oil floods. 


SOLTROL C FLoops 
Results of oil displacement runs 











SURFACTANT SOLUTION 


—SOLTROL C CONTACT ANGLES 
Average Contact Angle, Degrees 





2-5 Min. 24 Hrs. 
58+10 25+10 
60+10 30+10 
35+10 25+10 
35+10 25+10 
65+ 5 S2z- § 
78+10 64+10 
80+10 65+10 
85+10 70+10 
140+10 125+10 
160+ 5 160+ 5 
160+ 5 160+ 5 
160+ 5 160+ 5 
160+ 5 160+ 5 


cent sat- 
given in 


on sand columns 100 per 
urated with Soltrol C are 
Figs. 1 through 7. 

The effect of wettability for high 
interfacial tension floods at a water 
injection rate of 50 ml/hour (ap- 
proximate linear advance of | ft/hour 
is shown in Fig. |. These data sug- 





SUMMARY OF HIGH VISCOSITY C 





gest that oil recovery at water break 
through js not affected by the ad- 
vancing contact angle if it is less 
than 90°. At the flood rate used, 
however, oil recovery after break 
through was greater for lower values 
of the contact angie. Oil recovery 
decreased at breakthrough for con 
tact angles above 90° as the angles 
increased. 

Water floods on a water-wet sys 
tem at four interfacial tensions and 
a 50 ml/hour injection rate are given 
in Fig. 2. These results indicate that 
high interfacial tension floods 
more efficient than low interfacial 
tension floods on water-wet or 
treated system. Interfacial 
has an opposite effect on oil-wet sys 
tems as shown in Fig. 3 where oi! 


are 


un 


tension 





DISPLACEMENTS 








Oil Water Silicone Interfacia 
Plood Viscosity Treatment Tensio 
No. Ratio Mg /Gm Dynes /cw 
$0 35.8 0 34.5 
$1 35.8 0 $4.5 
$2 35.8 0 34.5 
$3 35.8 0 16.8 
54 35.8 8] 16.8 
55 35.8 0 16.8 
56 35.8 0 16.8 
$7 35.8 0 16.8 
58 35.8 0 1.2 
59 35.8 0 1.2 
60 35.8 0 1.2 
61 35.8 0 1.2 
62 35.8 0 1.2 
63 35.8 0.025 34.5 
64 35.8 0.025 34.5 
65 35.8 0.025 34.5 
66 35.8 0.025 16.8 
67 35.8 0.025 16.8 
68 $5.8 0.025 16.8 
69 35.8 0,025 16.8 
70 35.8 0.025 1.2 
71 35.8 0.025 1.2 
72 35.8 0.025 1.2 
73 35.8 0.025 1.2 
: 74 35.8 1.0 34.5 
: 75 35.8 1.0 34.5 
7¢ 35.8 1.0 34.5 
77 35.8 1.0 16.8 
78 35.8 1.0 16.8 
79 35.8 1.0 16.8 
80 35.8 1.0 1.2 
81 35.8 1.0 1.2 
82 35.8 1.0 1.2 
63 296 0 35.8 
84 296 0 35.8 
85 296 0 1.8 
86 296 0 1.8 
87 296 0.025 35.8 
88 296 0.025 35.8 
89 296 0.025 1.8 
30 296 0.025 1.8 
91 296 1.0 35.8 
j 92 296 1.0 35.8 
93 296 1,0 35.8 
94 296 1.0 1.8 
95 296 1.0 1.8 
96 296 1.0 1.8 
97 296 1.0 1.8 


* The scaling factor enon yr is for contact angles below 90°. 
9 


For contact angles above 90° the scaling factor is LViqy/j,ctose. 
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PORE VOLUMES FLOOD WATER INJECTED 


Fic. 1— HIGH INTERFACIAL TEN- 

SION FLOODS ON VARYING WETTA- 

BILITY SOLTROL C SySTEMS AT 50 
ML/HR INJECTION RATE. 


recovery increased as the interfacial 
tension was lowered. These results 
are as predicted by Equation 1. In 
the case of the water-wet system the 
spontaneity of oil displacement by 
water increases with ocos®. For the 
oil-wet system, where cos is nega- 
tive, less energy is required to dis- 
place oil as acosQ becomes smaller 
numerically. 

The rate of water injection has an 
appreciable effect on oil recovery, 
especially ai very low rates. This is 
shown for high and low interfacial 
tension floods on water-wet systems 
in Fig. 4. Oil recovery increased with 
rate up to a certain point, beyond 
which further rate increases had less 
effect on oil recovery. This stabilized 
recovery was reached at 15-25 ml/ 
hour for high interfacial tension 
floods and 100 to 200 ml/hour for 
low interfacial tension floods. A max- 
imum in oil recovery at water break- 
through was observed at 25 ml/hour 
for the high interfacial tension floods. 
This maximum in oil recovery for 
high interfacial tension floods as a 
function of rate may be explained by 
the tendency of water to imbibe 
into the smaller or tighter channels. 
At low flood rates flow channels are 
established by imbibition through the 
smaller channels. At slightly higher 
rates more water pushes into the 
larger channels as the flow capacity 
of the imbibition channels is ex- 
ceeded. At still higher rates more 
of the injected water is carried by 
the larger channels, with the result 
that some oil is trapped and by- 
passed in the tight channels. The 
maximum oil recovery is observed, 
then, at the injection rate where 
water flows at the same rate of linear 
advance through large and small 
channels. In the case of low inter- 
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facial tension floods, channels once 
formed have less tendency to branch 
or expand due to lower capillary 
pressure difference between phases. 
This means a higher injection rate is 
required to create a flow pattern or 
system of channels equivalent to the 
high interfacial tension floods. 

Fig. 5 shows that oil recovery in- 
creased with flood rates for both 
high and low interfacial tension 
floods on intermediate wettability 
systems. Stabilized recovery was 
reached at a rate of from 1 to 5 
ml/hour for high interfacial tension 
floods and at 5 to 25 ml/hour for 
low interfacial tension floods. Max- 
imum oil recovery was obtained at 5 
ml/hour for the high and at 2° 
ml/hour for the low interfacial ten- 
sion floods. For the neutral or inter- 
mediate wettability system there is 
no tendency for water to imbibe into 
the tighter channels and by-pass oil 
in the larger channels. As a result, 
water moves through the larger chan- 
nels as well as the small channels in 
a piston-like displacement. This re- 
sults in a low value of injection rate 
for stabilized recovery as compared 
to the water-wet system or the oil- 
wet system. The maximum oil recov- 
ery observed at the optimum injec- 
tion rate, as well as the lower rate 
required for stabilized oil recovery 
with high interfacial tension floods, 
is believed due to the explanation 
offered in the preceding paragraph. 

The effect of rate on oil recovery 
for oil-wet systems is presented in 
Fig. 6. Oil recovery became stable 
ct about 100 ml/hour for both high 
and low interfacial floods. For oil- 
wet systems water does not imbibe 
but must be forced into all flow 
channels. Less pressure is required 
to penetrate large channels than small 
channels. At low rates water flows 
through the larger channels, while 
higher rates result in the penetration 
of smaller channels. As a result, oil 
recovery increases with flood rate. 
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Fic. 2—VARYING INTERFACIAL TEN- 

SION FLOODS ON WATER-WET SOL- 

TROL C SystTeEMs AT 50 ML/Hour 
INJECTION RATE 
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Fic. 3—VARYING INTERFACIAL TEN- 
SION FLOOLS ON OIL-WET SOLTROI 
C Systems aT 50 MiL/Hour 
INJECTION RATE. 


rhe effects of contact angle, inter- 
faciai tension, and rate on oil recov- 
ery at water breakthrough are sum- 
marized in Fig. 7. These data show 
that oil recovery for intermediate 
wettability systems is less sensitive 
to rate of water injection and to the 
oil-water interfacial tension than for 
either the water-wet or oil-wet sys- 
tems. At rates above stabilized val- 
ues, oil recovery was highest for low 
interfacial tension floods on _ inter- 
mediate wettability systems. 

Several different combinations of 
variables were grouped in an effort 
to reduce experimental data to a 
common curve. One attempt is shown 
in Fig. 8 as LVy.o/p,cosO for con- 
tact below 90° and LVyu./p,0cosO 
for contact angle above 90°. These 
values are plotted against residual oil 
saturation at water breakthrough. 


HiGH Viscosity O1t FLoops 


Water floods were made on three 
using 32 cp oil. 
tensions and sev- 
red. Results are 


wettability system: 
rhree 
eral rates were 


interfacial 


presented in Table 4. 

The 32 cp oil systems are charac- 
terized by higher oil recoveries at 
water breakthrough at low rather 
than high flood rates. This result 
appeared visually to be due to chan- 
neling. Engelberts and K!tinkenberg”, 
who observed similar results, called 
the phenomenon “viscous fingering.” 
Jones-Parra™ reported this to be due 
to the fact that the inlet saturation 
distribution does not reach equilib- 
rium at high rates. This channeling 
effect was not so pronounced on oil- 
wet systems. It was especially appar- 
ent for low interfacial tension floods 
on untreated and lightly silicone 
treated sands at higher injection 
rates. The severe channeling for low 
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interfacial tension floods is probably 
due to the low capillary pressure dif- 
ference between the oil and wate 
phases. This results in less channel 
branching than in the case of high 
interfacial tension floods. 

For the untreated systems, a com- 
parison of the average of residual 
saturations at 19, {00, and 340 ml! 
hour rates shows ihe high interfacial 
tension floods to be slightly less effi- 
cient than low interfacial tension 
floods at water breakthrough but 
more efficient at 15 pore volumes of 
injected water. Such a comparison 
for the 0.025 mg/gm treated systems 
;or neutral wettability) does not re- 
veal any significant effect of inter- 
facial tension on oil recovery. For 
oil-wet systems, oil recovery was 
higher at water breakthrough and 
at 15 pore volumes of injected water 
the lower the interfacial tension. 

The results of water floods on 
three wettability systems saturated 
with 265 cp oil are given in Table 4 
Oil recovery at water breakthrough 
decreased with rate for the untreated 
and 0.025 mg/gm treated systems 
whereas oil recovery increased with 
rate for the oil-wet system. 

Attempts to correlate data for the 
high viscosity oils with Soltrol C 
were not successful except for oil- 
wet systems where the general group- 
ing shown in Fig. 9 was obtained. 
This grouping shows the residual oil 
saturation at eight pore volumes of 
injected water for Soltrol C and 15 
pore volumes of injected water for 
the high viscosity oils plotted against 
LV p./ p,ecosO. 


SUMMARY AND CONCLUSIONS 


Water floods have been made on 
laboratory-prepared unconsolidated 
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sand models to study the effects of 
the solid-water-oil contact angle, the 
oil-water interfacial tension, the 
water injection rate, and oil-water 
viscosity ratio on oil recovery. Ad- 
vancing contact angles ranged from 
35° to 160°, interfacial tensions from 
1.2 to 44 dynes/cm, water injection 
rates from 1.2 to 400 ml/hour, and 
oil-water viscosity ratios from 1.46 
to 296. The sand models of 59 cm 
length and 2.52 cm diameter had av- 
erage permeabilities of 11 darcies 
and average porosities uf 33.5 per 
cent. 

With the limitation that the con- 
tact angles and interfacial tensions 
were static values measured external 
to the flow models, the foilowing 
conclusions have been drawn: 


1. Wettability in the displacement 
of oil by water is important in that 
it determines the manner in which 
interfacial tension and flood rate af- 
fect oil recovery. 

2. Oil recovery increases with in- 
terfacial tension for water-wet sys- 
tems and decreases with interfa-ial 
tension for oil-wet systems. Neuiral 
wettability systems are les® sensit‘ve 
to interfacial tension than either o:- 
wet or water-wet Cases. 

3. For low oil-water viscosity ra- 
tios, oil recovery increases wi’h i- 
jection rate up to a point where fur- 
ther increase in rate results in no 
additional oil recovery. Neutrai wet- 
tability systems stabilize at Icwer in- 
jection rates than either oil-wet or 
water-wet systems. 

4. High viscosity oil systems are 
characterized by higher oil recoveries 
at low rather than high flood rates. 

5. A procedure has been developed 
for treating silica sand to produce 
surfaces of measurable wettability. 
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NOMENCLATURE 


o = Oil-water interfacial ten 
sion, dynes ‘cm. 

o,, = Solié-oil interfacial ten 
sige, Cynes/cm. 

o.~ = Solid-water interfacial ten 
sion, dynes/cm. 

© = Contact angle between 
solid and the water-oil 
interface measured 
» through the water phase 
P. = Capillary pressure, dynes 
cm’. 
to = Oil viscosity, cp. 
tty = Water viscosity, cp. 
W, = Cumulative water injected 
pore volumes. 
V = Water injection rate, cm 
min/cm’. 
L = Length of system, cm. 
r = Pore radius, cm. 

B.T. = Water breakthrough. 

A, = Adhesion tension between 
solid and water, ergs 
cm’, 

A... = Adhesion tension between 
solid and oil, ergs/cm 

AF = Free energy of displace 
ment of oil by wate: 
from a solid surface 
ergs/cm’. 
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ABSTRACT 


A new electric logging method, 
called displacement logging, often 
gives a direct indication of the pres- 
ence of mobile hydrocarbons in 
hydrocarbon-bearing formations. 
This method is based on the detec- 
tion fluids in these cases. The exist- 
which is formed between the invad- 
ing mud filtrate and virgin forma- 
tion fluids in these cases. The exist- 
ence of this phenomenon has been 
verified by laboratory studies and 
by actual field logs taken in known 
oil sands. 


INTRODUCTION 


During the past decade increas- 
ing use has been made of conven- 
tional electric logs to determine the 
amount of oil or gas in formations 
traversed by a borehole. This de- 
termination requires knowledge of 
several factors, such as true re- 
sistivity of the formation, formation 
factor, and water resistivity. Some- 
times all of the factors necessary to 
make a reliable saturation calcula- 
tion are not available. Also, the 
computations and corrections which 
are necessary in such a study often 
require the services of an expe- 
rienced log analyst. It has long been 
the desire of exploration people to 
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have a logging device which will 
locate hydrocarbon-bearing forma- 
tions directly without requiring com- 
plex interpretation methods. The 
displacement log, while not com- 
pletely fulfilling this desire, appears 
to be a step in that direction. 


BASIS FOR THE DEVICE 

The basis for this device stems 
from a phenomenon which occurs 
when mud filtrate invades a forma- 
tion which contains displaceable hy- 
drocarbons. As the mud filtrate in- 
vades the formation, most of the oil 
or gas is pushed out of the invaded 
zone, leaving residual oil. However, 
in addition to displacing most of the 
oil or gas, the connate water is also 
swept out, since the originally ir- 


reducible water is mobile in the 
presence of a water phase.”** Con- 
sequently, this irreducible water 


piles up as a bank between the 
virgin oil zone and the invading fil- 
trate. 

If the mud filtrate is of higher re- 
sistivity than the formation water— 
as is commonly the case—the re- 
sulting annular zone, containing the 
bank of formation water, will be of 
lower resistivity than either the 
virgin oil-bearing formation resistiv- 
ity or the invaded zone resistivity. 
This resistivity configuration, going 
from high to low to high again as 





‘References given at end of paper. 


the distance from the borehole in 
creases, is found only in those zones 
which contain displaceable hydro 
carbons. A device which would re 
liably detect this low resistivity zone 
would give a direct indication of 
mobile hydrocarbons. 


ILLUSTRATIONS 

This phenomenon is illustrated in 
Fig. 1 and 2. Fig. 1 shows a sketch 
of a zone containing salty connate 
water plus displaceable hydrocar 
bons, along with graphs showing 
fluid content and electrical resistivity 
vs distance from the borehole, be 
fore invasion occurs. Fig. 2 shows 
this zone, together with fluid con 
tent and resistivity graphs, after in- 
vasion by a fresh water mud. The 
resistivity pattern on moving away 
from the borehole is high: low: high 

The absence of the low zone 
phenomenon in water zones is illus 
trated in Fig. 3 and 4. Fig. 3 
shows a sketch of a zone containing 
salty formation water and residual 
oil, along with graphs showing fluid 
content and electrical resistivity vs 
distance from the borehole, before 
invasion occurs. Fig. 4 shows this 
water zone, together with fluid con- 
tent and resistivity graphs, after in- 
vasion by a fresh water mud. In 
this case the resistivity pattern on 
moving away from the borehole is 
high: low since there is no change 
in water saturation at the head of 









the invading front. The same re- 
sistivity pattern holds for a 100 per 
cent watur saturated zone. 

This low resistivity zone, or “low 
zone,” between the invaded zone 
and virgin formation then is found 
only in zones containing displace- 
able hydrocarbons. It has been 
studied in detail in the laboratory by 
core and model studies and has been 
found to exist in the field in zones 
which later proved to be productive. 


CORE STUDIES 

Initial laboratory work was per- 
formed to determine the efficiency 
of fluid displacement in the invaded 
zone and to determine the resistivity 
contrast that would exist under 
various conditions. Fig. 5 shows the 
resistivity variations along a sand 
core which had been reduced to con- 
nate water by flushing with oil and 
then flooded with fresh water. The 
“low zone” phenomenon has been 
observed in unconsolidated sands, 
consolidated sand cores, and in both 
homogeneous and vugular limestone 
cores. 


RESISTIVITY VARIATION 

The resistivity variation in any 
given case will depend on the water 
saturation in the virgin formation, 
the resistivity of the formation 
water, the resistivity of the invading 
filtrate, and the residual oil satura- 
tion in the invaded and “low zones.” 
The displacement of the connate 
water by the filtrate is not 100 per 
cent efficient,. so the contrast be- 
tween the “low zone” and the in- 
vaded zone is reduced because of 
some mixing of the saline and fresh 
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waters. It is believed that in sand- 
stone reservoirs typical resistivity 
ratios may range between 2:1:10 and 
20:3:10 for the invaded zone: “low 
zone”: virgin formation resistivities. 

The width of the low resistivity 
zone depends on the connate water 
saturation and the depth of inva- 
sion. Larger values of connate water 
saturation and greater depths of in- 
vasion lead to wider low resistivity 
zones. For 20 per cent connate 
water and five hole diameters in- 
vasion, “low zone” widths of about 
one-half of a hole diameter are be- 
lieved to exist in sandstone reser- 
voirs. 


STABILITY OF RESULTING 
RESISTIVITY CONFIGURATION 

Another factor that was consider- 
ed in initial laboratory studies was 
the stability of the resulting resistiv- 
ity configuration. Three processes 
that tend to eliminate the resistivity 
contrast are: (1) diffusion of the 
bank of salty formation into the in- 
vaded zone, (2) imbibition of the 
bank of salty formation water into 
the virgin, oil-bearing formation, 
and (3) gravitational separation of 
the three zones due to their different 
densities. 

Limited investigation indicates 
that unless formation permeabilities 
are of the order of 1 darcy or higher, 
the effects due to diffusion and 
gravitational separation can usually 
be neglected for periods of several 
weeks. One test of the stability of 
the resistivity contrast was made on 
a 400 md Torpedo sand section, ap- 
proximately '%4 in, thick, over a 


period of 13 days. 


sand was 


This particular 
since other data 
indicated thit the imbibition effect 
would be Fig. 5 also shows 
the change in the resistivity profile 
due to all causes in the Torpedo 
sand section. From this we expect a 
Significant resistivity contrast to exist 
Over a period of several weeks after 
an oil zone is penetrated. 


chosen 


large. 


MODEL STUDIED 


WEDGE MopeE! 

After being satisfied that the effect 
noted in the laboratory appeared 
promising as a possible means for 
mobile hydrocarbons di- 
rectly, attention was directed toward 
a method of recording the presence 
of the “low zone” by some type of 
logging device to be run in the bore- 
hole. We have used the term “dis- 
placement log” to indicate any log- 
ging method which is designed to 
give an indication of the “low zone” 
formed due to displacement of oil 
mud filtrate. Since, gen- 
erally, the horizontal depth of in- 
vestigation of conventional logging 
arrangements increases with spacing, 
recording a number of conventional 
spacings appeared feasible. It was 
hoped that the resulting profile of 
apparent resistivity vs spacing would 
detect the low resistivity zone by 
showing a minimum. A unique effect 
in the apparent resistivity profile is 
since it appears 
impossible, or at least impractical, 
to apply reliable corrections to con- 
vert the apparent resistivity values 
profile. 
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An electrolytic model was used to 
simulate a formation with a “low 
zone.” Different salt water solutions 
and agar mixtures were used to set 
up the various resistivity configura- 
tions desired. A Lucite section was 
constructed to simulate the logging 
sonde. 

Initial electrolytic model studies 
conducted with a series of two- 
electrode spacings did not result in 
profiles capable of detecting the 
presence of the “low zones” set up 
in the model. However, three-elec- 
trode spacings did appear to have 
promise. It was found that in order 
to detect wide ranges of the re- 
sistivity contrast expected in nature, 
effects due to the flow of current in 
the mud-filled borehole would have 
to be reduced or eliminated, if pos- 
sible. Insertion of a series of insulat- 
ing discs or bumpers between the log- 
ging electrodes accomplished a sig- 
nificant reduction in borehole effects 
when the bumpers approached 80 
to 90 per cent of the borehole 
diameter. 

Much of the early work was done 
with this electrolytic model. This 
work indicated that many of the 
resistivity contrasts expected in na- 
ture could be detected with a system 
of three-electrode spacings employ- 
ing bumpers. Problems with dif- 
fusion between zones of different 
resistivity in the model and the 
relative cost of producing a great 
many model runs led to the use of 
a resistivity analog model as a more 
reliable and cheaper method to 
evaluate the usefulness of various 
systems for displacement logging. 


BORE HOLE 





FORMATION 
SALT WATER 


DISTANCE FROM HOLE 


Pa ee 


EERE cent 
DISTANCE FROM HOLE 
Fic. 4—WaATER ZONE AFTER 
INVASION. 





ELECTRICAL 
RESISTIVITY 








VOL. 204, 1955 


—— INITIAL INVASION 
~--1!3 DAYS AFTER INITIAL | 
INVASION 


OM 


APPARENT RESISTIVITY— 
3 











DISTANCE - INCHES 
Fic. 5 — Resistivity vs DISTANCE 
ALONG OIL-BEARING SANDSTONE 
CoRE INVADED BY FRESH WATER. 


RESISTIVITY ANALOG MODEL 

The analog model used was de- 
signed by Hubert Guyod, Houston 
logging consultant, to simulate earth 
formations and mud-filled boreholes 
by means of an ciectrical resistance 
network.* 

Resistance networks were built 
that would permit various types of 
displacement logs (apparent resistiv- 
ity for a given electrode arrange- 
ment vs electrode spacing) to be run 
at different depths in a simulated 
formation. It was possible to make 
the runs with insulating bumpers 
which were 100, 90 and 84 per cent 
of the borehole diameter. 

In the work covered by this paper, 
the depth reference point is taken 
as the position of the current or sin- 
gle electrode in a three-electrode ar- 
rangement. It was found that this 
method of presentation resulted in 
an unique indication of the “low 
zone.” The use of the midpoint of 
the two potential probes as_ the 
depth reference point, as is done 
in conventional logging, may also 
be a satisfactory method. The latter 
method is probably easier to do 
instrumentally. However, the basic 
object of displacement logging, de- 
tecting the “low zone,” remains the 
same no matter how the data are 
presented for log interpretation. 

In view of all model work to 
date, we define an indication of the 
“low zone” to be either (1) a min- 
imum in the graph of apparent re- 
sistivity vs electrode spacing which 
is stationary, i.e., it continues to 
occur at the same spacing for re- 
sistivity profiles run at different 
depths in the zone being logged; or 
(2) a stationary inflection in the 
graph of apparent resistivity vs elec- 
trode spacing when the apparent re- 
sistivity continues to increase with 
spacing. 

A case is shown in Fig. 6 of a 


displacement log taken with a three- 
electrode configuration, at 1-ft inter 
vals opposite a 10-ft thick bed (as 
suming a hole diameter of 8 in.) 
The bed of interest has a resistivity 
contrast of 1:10:1:30 for the mud, 
invaded zone, “low zone,” and virgin 
formation, respectively. The su: 
rounding bed resistivity, R,, is equal 
to three times the mud resistivity 
Each profile is plotted with the posi- 
tion of the current electrode as the 
depth reference point. 

It is important to note that the 
indication due to the “low zone” has 
the unique characteristic of remain 
ing On the same spacings even while 
the sonde is moving through the bed 
of interest. Using the current refer- 
ence method of plotting, a minimum 
or inflection due to bed effects would 
move across the profile as the sonde 
moved. It is this characteristic that 
distinguishes the “low zone” indica 
tion from any inflection or minimum 
due to bed boundaries or from laye! 
ing occurring within a bed. 

Over 200 model runs have been 
made to study the effect of changing 
bed thicknesses and resistivity con 
trasts in the invaded zone, “low 
zone,” virgin formation, and adja- 
cent beds using both three-electrode 
and differential three-electrode con 
figurations. 


RESULTS OF MopeL WorkK 

The results of the model 
may be summarized as follows: 

1. If an indication of the “low 
zone” is obtained on a displacement 
log, it is unique. Bed boundary and 
layer effects move through the re 
sistivity profile when the sonde moves 
past the bed. 

2. It has not proved possible to 
date to design a logging system that 
will detect the “low zone” in all! 
cases. The ability to detect the “low 
zone” is dependent on the resistivity 
values of the invaded zone, the “low 
zone,” the virgin formation, and the 
adjacent formations. It also is de- 
pendent on the bumper diameter of 
the sonde, the thickness of the bed 


work 


"ZON OmMOP<z- 


APPARENT RESISTIVITY 








SPACING - - HOLE DIAMETER «® 
(SPACINGS EXTEND FROM 6” TO 82°) 


Fic. 6—Mope. RESISTIVITY 
PROFILES. 





of interest, and the type of electrode 
configurations used. 

3. It appears that at least 60 to 70 
per cent of producing sand reservoirs 
can be detected using 90 per cent 
bumpers. In certain areas the suc- 
cess ratio is expected to be higher, 
and in others, lower. For limestone 
reservoirs, model work indicates that 
detection of the “low zone” will be 
much more difficult due to the heter- 
ogeneous nature of these beds. 


FIELD TESTS 


A field program undertaken after 
core studies and early model work 
proved the existence of the “low 
zone” in actual wells. 


SPECIAL INSTRUMENTATION 

The main problems in special in- 
strumentation were with the physical 
dimensions of the sonde and the nec- 
essary electrical switching that per- 
mitted readings to be taken on a 
plurality of electrode spacings. The 
electrode sections of the sondes, built 
for use in 8% in. and 9% in. bore- 
holes, were like the sketch shown in 
Fig. 7. Note that the electrode sec- 
tion of the sonde provides a mud by- 
pass through the center to reduce 
swabbing effects. 

Leaf spring centralizers were used 
to prevent flexing and excessive wear 
on the thin rubber bumpers and to 
prevent excessive scraping of filter- 
cake from the wall of the hole. In- 
side a pressure tight housing above 
the electrode arrangement, means 
were provided to switch the recorder 
to the appropriate spacing. 


LOGGING PROCEDURES 

The initial field tests, run in shal- 
low holes in East Texas, were made 
by a number of runs past the zone 
of interest. Each successive run was 
made with a different spaced three- 
electrode arrangement. 

Later logs were made by stopping 
the sonde opposite the formation of 
interest and successively switching 
the recording device through the vari- 
ous spacings. After all readings were 
taken at one depth, the sonde was 
usually dropped about 1 ft, and an- 
other series of readings was taken. 
Displacement logs were often taken 
opposite uniform shale formations 
to check instrument performance. 


EXPERIMENTAL DATA 
AND Discussion 

Fig. 8 is a portion of a conven- 
tional electric log survey from 5,150 
to 5,440 ft in an Atlantic Refining 
Co. well, Sholem Alechem field, Step- 
hens County, Okla. Fig. 9 shows dis- 
placement logs taken opposite a zone 


236 


which later produced pipe line oil. 
An indication like those obtained in 
model runs due to the “low zone” is 
easily seen in the region of spacing 
positions 8, 9, and 10. 

Displacement logs were taken in 
the same well opposite a water sand 
at depths of 5,219, 5,220 and 5,221 
ft. These logs are shown in Fig. 10. 
The displacement logs show clearly 
that R, is greater than R,, and that 
no indication of hydrocarbons is 
present. 


LoGcs IN SALT, LAKE FIELD 

Fig. 11 shows a portion of the 
conventional electric log and micro- 
device in an Atlantic Refining Co. 
well, Salt Lake field, Aransas County, 
Tex. Some displacement logs taken 
in this section are shown in Fig. 12. 

The high resistive streak shown at 
7,795 ft on the micro-device moves 
through the spacings as the depth of 
the sonde is changed. Any effect due 
to streaks of different resistivity oc- 
curring in the sand will not remain 
on the same probe readings. It is 
probable that the displacement logs 
indicate an inflection in the region of 
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Fic. 8 — ELECTRICAL SURVEY, AT- 
LANTIC REFINING Co. WELL, SHOLEM 
ALECHEM FIELD, STEPHENS COUNTY, 
OKLA. 
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Fic. 9 DISPLACEMENT LOGGING 

RESULTS. WELL LOCATED IN SHOLEM 

ALECHEM FIELD, STEPHENS COUNTY. 
OKLA. 


Remarks: Dips from spacings 8 to 11 in- 
dicated movable hydrocarbons. This zone pro- 
duces clean oi! from perforations 5,316 ft to 
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Fic. 10 DISPLACEMENT LOGGING 

RESULTS. WELL LOCATED IN SHOLEM 

ALECHEM FIELD, STEPHENS COUNTY, 
OKLA. 


Remarks: Indication of a water zone 
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Current probe at 7,794 ft 

Fic. 12 — DisPLACEMENT LOGGING 
RESULTS. WELL LOCATED IN SAL1 
LAKE FIELD, ARANSAS COUNTY, TEX. 


Remarks: An indication of oil is noted 
around spacings 8, 9, and 10. A resistivity 
streak masks the indication at 7,792 ft. This 
sand produces in this field. 


probe positions 8, 9 and 10. Though 
it is not the producing sand in this 
well, this sand produces oil in this 
field and in this well the sand is 
above the oil-water contact. Shows 
in sidewall cores were obtained in 
this interval. 


FIELD RESULTS 

To date, displacement logs have 
been successfully recorded in 10 dif- 
ferent Atlantic wells. A total of 120 
intervals were logged, 25 of which 
were in known or probable hydro- 
carbon-bearing sands. In these hy- 
drocarbon-bearing sands, 19 inter- 
vals were interpreted as having an 
indication due to the “low zone.” 
(Most of the remaining ones were 
distorted by resistivity inhomogenei- 
ties which made interpretation im- 
possible.) Of a total of 88 water for- 
mations, or formations of unknown 
fluid content, which were logged, 
only four were interpreted as hav- 
ing an indication due to the “low 


Most of the information supplied 
in the Martin-Campbell paper is 
based upon work of The Atlantic 
Refining group prior to the time a 
logging system capable of record- 
ing simultaneously a multitude of 
three-electrode resistivity curves was 
available. 

Atlantic’s field work, their own 
model study and that with the 
Guyod model, indicated that the 
logical log should consist of a fam- 
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TABLE 1 — FIELD DATA ON DISPLACEMENT LOGS RECORDED SUCCESSFULLY IN TEN ATLANTIC WELLS 


Known or Known or Shale Presumably 
Probable Probable Intervals Water or Non-hydrocorbon 
Number Oil or Gas Oil or Gas logged, No Unknown Intervals in 
Intervals Intervals Intervals Indications of Intervals Which Oi! Wo 
Weill Logged Logged Detected Oil or Gas Logged Indicated 
1. East Texas 2 i 1 0 1 0 
2. Southwest Texas 16 0 0 2 14 0 
3. North Lovisiana 22 0 0 1 20 1 
4. East Texas 12 0 0 0 12 1 
5. Oklahoma 20 7 7 tt) 13 0 
6. Southwest Texas 13 6 3 0 7 0 
7. Southwest Texas 11 7 5 2 3 1 
8. East Texas 6 1 1 1 4 1 
9. Southwest Texos 7 3 2 0 4 0 
10. Mississippi 1 0 0 1 10 0 
TOTALS 120 25 19 7 88 4 


zone.” None of these four zones has 
been tested. No stationary minimum 
or inflection was found in any shale 
formation. Table 1 lists these field 
data. 

It is often difficult to interpret the 
displacement logs when resistivity 
inhomogeneities are present. This 
means that in beds which are not 
homogeneous over a sufficient thick- 
ness to permit reliable interpretation, 
the displacement logging method is 
limited. 

The field tests have shown that: 

1. It is possible to obtain a dis- 
placement log under field conditions. 

2. Indications are found in hydro- 
carbon-bearing sands that are simi- 
lar to those recorded in model set- 
ups. 

3. Resistivity inhomogeneities 
make interpretation difficult in many 
cases. 


CONCLUSIONS 

The main conclusions reached are 
as follows: 

1. The setting up of the “low 
zone” is a real phenomenon and the 
resistivity configuration remains sig- 
nificant over a period of several 
weeks. 

2. The ‘“‘low zone’ causes an 
unique effect, which can often be de- 
tected, on displacement logs. 

3. The “low zone” has been de- 
tected under field conditions. 


DISCUSSION 


WALTON J. GREER 
MEMBER AIME 


HOUSTON, TEX. 
ily of about eight three-electrode 
curves with spacings ranging from 
10 in. to 60 in. An instrument to 
provide such a log has been built. 

Fig. | shows a biock diagram of 
the system. Two electric fields are 
established in the formation by con- 
stant currents of two different fre- 
quencies flowing between the two 
pairs of current electrodes: C—C, 
and C —C.. Each potential electrode 
disposed along the tool measures 


HALLIBURTON OIL WELL 






We believe that displacement log 
ging is a method which will often 
indicate the presence of mobile hy 
drocarbons directly. Such a method 
should materially aid in the searc! 
for oil. 

Certain improvements in the dis 
placement logging instrumentation 
and equipment need to be made. 
however, before such a method wil! 
be practical. An improvement that 
would permit better discrimination 
between moving and stationary in 
flections and minima is the simulta 
neous recording of the apparent re 
sistivities from different spaced elec 
trode arrangements. A service com 
pany has entered an agreement with 
The Atlantic Refining Co. to develop 
such a system. 
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potential signals from each of these 
fields according to the spacing be 
tween it and the current electrodes 
Thus P, receives a signal at 133 
cycles per second corresponding to 
that of a three-electrode configura 
tion of 10 in., and a signal at 200 
cycles per second corresponding to 
that of a three-electrode configura- 
tion of 16 in. In a similar manner 
P, receives a signal at 133 cycles per 
second corresponding to that of a 
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three-electrode configuration of 22 
in. spacing and a signal at 200 cycles 
per second corresponding to that of 
a three-electrode configuration of 28- 
in. spacing. Similarily, each of the 
other potential electrodes receives 
signals corresponding to two three- 
electrode configurations. The signals 
picked up by these potential elec- 
trodes are used to modulate in fre- 
quency various carriers. The carrier 
signals then are injected into the 
single conductor line circuit for 
transmission to the surface. The in- 
dividual signals are separated there 
and the log recorded in the conven- 
tional fashion. 

Fig. 2 shows a portion of a log 
made with this system recorded on 
conventional scales. For clarity of 
the record each of the curves is 
based at a different vertical rule line. 
Thus, the 10-in. curve is recorded on 
a scale of 20 ohm-m with the zero 
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Fic. 2—PorTION OF LoG MADE WITH 
NeW SYSTEM RECORDED ON 
CONVENTIONAL SCALES. 











Fic. 3—PoRTION OF LOG MADE WITH 
NEW SYSTEM RECORDED ON 
CONVENTIONAL SCALES. 








at the extreme left vertical rule line. 
The 16-in. curves, 22-in. curves, etc., 
are zeroed along the vertical rule 
lines indicated at the upper part of 
the figure. 

To fully evaluate the displacement 
concept it is necessary that a plot be 
made of the resistivity profile as illus- 
trated in Martin’s presentation. It is 
interesting to note, however, how 
clearly the invasion picture can be 
shown by this family of three-elec- 
trode curves. At the depth of 7,326- 
51 there is a water sand. The shorter 
spaced curves are influenced consid- 
erably by the invaded zone, while 
the longer spaced curves are influ- 
enced less until at about 48-in. spac- 
ing there is no apparent influence. 

Fig. 3 shows another section of 
this same log. Here are included 
some thin streaks, and a very good 
presentation is given of the responses 
of the three-electrode resistivity meas- 
uring configurations of different spac- 
ings to thin sections. At a depth of 
8,365 there is a very thin section 
and for each spacing the character- 
istic lag or shadow zone is apparent. 
At a depth of 8,322 there is a streak 
of about 4 ft thickness. Here can 
be noted the characteristic response 
of the configuration to a bed of crit- 
ical thickness. 

For all evaluations of the displace- 
ment logging operations to date it 
has been necessary to plot the re- 
sistivity profile at a given depth for 
the various spacings. Until recently 
this plot was made manually. A 
mechanical reader has been devel- 
oped which facilitates rapid plotting 
of the profile directly from the log 
of the eight three-electrode curves. 
Fig. 4 and Fig. 5 show such profiles 
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of a water sand and an oil sand re- 
spectively. Even with the mechanical 
reader considerable labor is required 
to put the data into proper form for 
final interpretation. It is anticipated 
that there eventually will be provided 
an instrument which will automat- 
ically plot the profiles as the log is 
being made so that better well-site 
interpretations can be made. 

To date 13 wells in South Texas, 
Louisiana and Oklahoma have been 
surveyed with the fully automatic 
and simultaneously recording equip- 
ment. 

Future 
curves of 
more logical sequence as 
by past field experience. 


instruments will provide 
spacings to 72 in. in a 
indicated 
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I would like, at the outset, to em- 
phasize that the effect which the au- 
thors call the low zone is indeed 
real. I note this because when the 
paper was presented verbally, it 
seemed that many listeners doubted 
its validity. Data were obtained on 
water flooding at the Pennsylvania 
State University just after the war 
which show the low resistivity zone 
very strikingly. In these experiments 
distilled water was used to displace 
oil from water-wet sands containing 
both oil and brine. The floods were 
linear and the resistivity along the 
direction of flooding was monitored 
continuously. The build-up of a bank 
of brine in which the oil saturation 
was small and, in consequence, the 
resistivity low, is clearly revealed. 

Campbell and Martin are to be 
congratulated for recognizing that 
this phenomenon can he applied to 
the problem of log interpretation. 
The practical method they suggest 


DISCUSSION 


M. R. J. WYLLIE 
MEMBER AIME PITTSBURGH, PA. 
of recognizing the low zone is in- 
genious and, I believe, workable in 
many instances. Nevertheless, it is 
my view (necessarily subjective) that 
the basic approach they propose is 
anachronistic. The recognition of the 
low zone, apart from such practical 
difficulties as stem from equipment 
and formation heterogeneity (espe- 
cially interstitial clay), depends upon 
a reasonable depth of filtrate inva- 
sion, a high ratio of R,, to R,, and 
adequate hydrocarbon mobility. In 
many formations, at least with prop- 
erly tailored mud, these conditions 
will be met. But this cannot be re- 
lied upon. If the method is not uni- 
versally applicable, its simplicity, in- 
stead of being advantageous, becomes 
the reverse. All sorts of well-meaning 
log interpretations will inevitably be 
made by semi-skilled interpreters who 
are impressed by the simplicity of 
the low zone technique but unaware 
of the physical processes on which 
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its validity rests. Mistaken interpreta- 
tions will be made and these will 
not necessarily err on the side of 
optimism. Management will not be 
impressed. 

It may be argued that the low zone 
technique will supplement other in- 
terpretation techniques. This would 
be admirable but there are already 
too many logs run. With all due 
deference to the authors, ! subscribe 
to the following view. The only basis 
on which log interpretation can rest 
satisfactorily is one which is uni- 
versally applicable. This means that 
directly or indirectly R,, R,, and the 
SP must be found (for aqueous-base 
muds). The resistivity parameters are 
best measured by focusing logs run 
in muds whose resistivities are specif- 
ically designed to the problems at 
hand, i.e., maximum R, likely to be 
encountered and values of R,,,/R 
suited to the logs. 


AUTHORS’ REPLY to M. R. J. WYLLIE 


It seems that M. R. J. Wyllie, 
though agreeing with the concept of 
the formation of the low zone, feels 
that displacement logging is already 
outmoded. He mentions what he con- 
siders to be factors limiting its appli- 
cability. 

Heterogeneities of the formation, 
in the form of d:fferent resistivity 
beds, will effect the recognition of 
the low zone. However, recent field 
tests show that the displacement log 
has apptication in beds of thickness 
down to 3 or 4 ft. Though we have 
no field data to prove it, a sand with 
heterogeneity due to interstitial clay 
(a shaly or “dirty” sand) may be 
subject to evaluation with a displace- 
ment log. If R,, is greater than R,, in 
an invaded shaly sand, the resistivity 
pattern on moving away from the 
borehole is high:low:high. This is 
precisely the condition for displace- 
ment logging. 

It is true that displacement logging 
must have a reasonable depth of fil- 
trate invasion. Probably just a few 
inches will do. It is also true that any 
interpretation method which uses R, 
must depend on a few inches filtrate 
invasion. 

A high ratio of R,, to R, is not 
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necessary for displacement logging. 
When &,, is greater, by any amount, 
than R,,, the necessary resistivity pat- 
tern is set up by invasion of a sand 
bearing mobile hydrocarbons. If the 
invaded zone resistivity is less than 
the true resistivity (with R,, > Ry), 
then we are assured that mobile hy- 
drocarbons are present. Model data 
indicate that such a resistivity config- 
uration gives the unique inflection 
necessary for interpretation. This is 
true even with no low zone as long 
as the invaded zone resistivity is less 
than the true resistivity. The fact 
that the low zone is present accentu- 
ates the indication. We do not think 
that mud resistivity control will be 
generally necessary. When necessary, 
it should be no more of a problem 
than the tailoring of the mud pro- 
posed in Wyllie’s “universal” method. 

As for adequate hydrocarbon mo- 
bility, that’s what we are looking for. 
If the oil saturation changes upon 
invasion from 70 per cent to 40 per 
cent, we would expect the resistivity 
of the virgin formation to be about 
three times that of the low Zone. 
Greater oil saturation changes would 
give greater contrasts. If the change 
in oil saturation is much less than 30 


per cent, the zone may not be com- 
mercial. Zones of this type would 
also be hard to find with Wyllie’s 
proposed method. 

We feel that one of the main ad 
vantages of the displacement log 
method is that interpretations can be 
made by semi-skilled interpreters. If 
the unique indication is present, the 
interpreter is certain of movable hy- 
drocarbons in the formation whethe: 
or not there is any other supporting 
evidence, such as “shows.” A logging 
expert can examine zones of possible 
interest which do not have displace- 
ment log indications by conventional! 
interpretation methods. We may say 
here that the displacement log would 
permit calculation of values for R, 
and R,. An SP will eventually be run 
with this device. The bumper sonde 
serves to an extent the same purpose 
as the pad on contact logs and the 
focusing systems in the focused de- 
vices used to determine R,. 

Like Wyllie, we too would like to 
use a method that is “universally ap- 
plicable.” Until such time as one is 
found, however, we will continue to 
use any log interpretation method 
that will aid in the search for mov- 
able hydrocarbons. wee 
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The method of comparing laboratory and field per- 
formance as set forth by the authors is certainly an 
interesting and thought provoking approach. The anal- 
ysis is, however, based on a hypothetical reservoir and 
data, and we believe one of the basic premises to be in 
error. 

The authors assumed a typical laboratory K,/K 
curve which was characterized by a critical gas satura- 
tion of 10 per cent of pore space. In order to further 
substantiate the validity of their assumption, the au- 
thors cited references which indicate that the critical 
gas saturation of typical laboratory K,/K, curves fall 
in the general range of 5 to 20 per cent of pore space. 


CONTINENTAL Oil CO 
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This is in disagreement with results obtained in our 
laboratory. 

Results of all our laboratory K,/K.,, tests would lead 
us to believe that the critical gas saturation is extremely 
small—tess than 2 per cent of pore pace in the majority 
of cores. The high critical gas saturations reported in 
the literature are due, we believe, to the inability of 
many laboratory K,/K, procedures to measure gas flow 
and saturations accurately in the lower ranges. If K,/K 
data originating at low gas saturations had been used 
by the authors, it is believed that the conformance of 
actual and theoretical reservoir performances would 
have been improved. kkk 
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(Dec. JPT) 
CANADA: crude oil, economic behavior 8 years after 
discovery and predictions of future growth, 
(Art.) July JPT, 7 
CARMAN FORMULA: validity for consolidated sands, 
(Tech. Note) Nov. JPT, 53 
CASE HISTORIES. See Field 
CASING: corrosion, external: extent, protective cathodic 
current required, and instrument for study, 
92 (June JPT) 
caused by electrochemical potentials, 66 (Apr. 
JPT) 
current from electric log SP curves: quantitative 
derivation of profiles, 66 (Apr. JPT) 
CAUDLE, B. H., et al.: Encroachment of Injected Fluids 
beyond the Normal Well Pattern, 79 (May 
JPT); discussion, 84, 85 (May JPT) 
CEMENT SLURRIES: high strength granular sealing ma- 
terial (walnut shells) beneficial in cement- 
ing casing and squeeze cementing, ( Art.) 
Aug. JPT, 15 
CEMENTING: radioactive tracers for location of cement 
top, 5 (Jan. JPT) 
CHAPMAN, A. H., and Duce, J. T.: Economics of the 
International Petroleum Industry, { Art.) 
March JPT, 9 
CHEEK, R. E., and MENziE, D. E.: Fluid Mapper Model 
Studies of Mobility Ratio, (Tech. Note) 
Nov. JPT, 49 
CoLLINGwoop, D. M., and KERN, C. P.: Discussion on 
Encroachment of Injected Fluids beyond 
the Normal Well Pattern, 83 (May JPT) 
COMPUTATION: presenting data in usable form for 
machines: tables for approximating smooth 
curves by polynomials of degree five or 
lower, (Tech. Note) Dec. JPT, 43 
CoMPUTER: IBM, (Art.) May JPT 9 
electrical analog: design for solving unsteady-state 
heat conduction problems of radial system 
containing a moving cylindrical heat source, 
205 (Dec. JPT) 
reservoir analyzer, Carter type, 73 (May JPT} 
CoNnLey, F. R., and FELSENTHAL, M.: discussion of 
Analysis of Reservoir Performance k,/k 
Curves and a Laboratory k,/k, Curve Meas- 
ured on a Core Sample, 240 (Dec. JPT) 
CONSTRUCTION: of filled island for offshore drilling, 
modification program, (Art.) Dec. JPT, 13 
CoRNELL, Davin: Calculations of Stabilized Gas Well 
Performance Curves from Back Pressure 
Test Data, (Tech. Note) May JPT, 53 
CorNELL, Davin, and SUKKAR, Y. K.: Direct Calcula- 
tion of Bottom-hole Pressures in Natural 
Gas Wells, 43 (March JPT) 
CORROSION: casing. See Casing 
Cores: fresh and aged, East Texas Woodbine, differ- 
ences, laboratory study, 86 (June JPT) 
CPC. See Calculator 










CraliG, F. F., Jr., et al.: Oil Recovery Performance of 
Pattern Gas or Water Injection Operations 
from Model Tests, 7 (Jan. JPT); discus- 
sion, 15 (Jan. JPT) 

CRAWFORD, P. B., and LANpRUM, B. L.: Effect of 
Drainhole Drilling on Production Capacity, 
(Tech. Note) Feb. JPT 45 

Effect of Unsymmetrical Vertical Fractures on Pro- 
duction Capacity, (Tech. Note) Apr. JPT, 
47 

CRAWFORD, P. B., et al: Estimating Size and Shape of 
Vertical and Horizontal Fractures, 175 
(Oct. JPT) 

Crockett, J. J., et al.: Studies of the Mechanism of 
Paraffin Deposition and Its Control, ( Art.) 
Sept. JPT, 23 

Croes, G. A., and SCHWARZ, N.: Dimensionally Scaled 
Experiments and the Theories on the Water- 
drive process, 35 (March JPT) 

CULLENDER, M. H.: /sochronal Performance Method 
of Determining Flow Characteristics of Gas 
Wells, 137 (Sept. JPT) 

CUNNINGHAM, R. A., and Murray, A. S.: Effect of 
Mud Column Pressure on Drilling Rates, 
196 (Nov. JPT); discussion, 203 (Nov. 
JPT) 

CuRVE FITTING. See Computation 


D 


De ANDRADE, P. V.: Result of SP Log Interpretation, 
(Tech. Note) Nov. JPT, 59 

DEPLETION EQUATION: derivation and solution, 121 
(Aug. JPT) 

De Witte, LEENDERT: Study of Electric Log Interpre- 
tation Methods in Shaly Formations, 103 
(July JPT) 

De WITTE, LEENDERT, and RApb, F. J.: Corrosion of 
Oil Well Casing by Earth Currents, 66 
(Apr. JPT) 

DIETZ THEORY On water-drive process: study by di- 
mensionally scaled experiments, 35 (March 
JPT) 

DISPLACEMENT LOGGING: detection of presence of mo- 
bile hydrocarbons in formations, laboratory 
and field studies, 233 (Dec. JPT) 

DouG.as, JIM, JR., et al.: Calculation of Unsteady-state 
Gas Flow within a Square Drainage Area, 
190 (Nov. JPT) 

DRAINHOLEs: effect on well productivity, model studies, 
30 (Feb. JPT) 

DRAINAGE BOUNDARIES: reservoirs at steady state, model 
studies, 182 (Nov. JPT) 

DRILLING: drainholes effect on production capacity, 
model studies, (Tech. Note) Feb. JPT 45 

drainholes, lateral holes, slant holes, definitions, 
(Art.) Feb. JPT, 11 

effect on well productivity of slant holes, drain- 
holes, and lateral holes, (Art.) Feb. JPT, 11 

offshore, from filled island, 70 wells, (Art.) Dec. 
JPT, 13 

radioactive tracers for discovery of lost circulation 
zones and casing leaks, 4 (Jan. JPT) 

rate: effect of mud column pressure, 196 (Nov. 
JPT) 

DRILLING FLUIDS: lignosulfonate stabilized emulsions, 
(Art.) June JPT, 9 

DRILLING MUD: effect of column pressure on drilling 
rates, 196 (Nov. JPT) 

Duce, J. T. and CHapMan, A. H.: Economics of the 
International Petroleum Industry, (Art.) 
March JPT, 9 

Dyes, A. B.: Discussion on Oil Recovery Performance 
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of Pattern Gas or Water Injection, 15 (Jan. 


JPT) 


K 

4 

EARLOUGHER, R. C. and AMSsTUTZ, R. W.: Water 
Source and Requirements, ( Art.) Jan. JPT, 
18 

Epwarops, J. M., and Simpson, A. L.: Method for Neu- 
tron Derived Porosity Determination for 
Thin Beds, 132 (Aug. JPT) 

EINARSEN, C. A.: High Strength Granular Sealing Ma- 
terial Increases Efficiency of Primary Ce- 
menting and Squeeze Cementing, (Art.) 
Aug. JPT, 15 

ELECTRIC LOGGING: displacement. See Disyiacement 

ELECTRIC LOGs: interpretation for low resistivity sands, 
(Art.) Aug. JPT, 10 

interpretation for shaly oil-bearing formations. 
quantitative, 103 (July JPT) 

ELECTROLYTES: lignosulfonates in solution, (Art.) June 
yea, 9 

ELK BASIN Tensleep reservoir: pressure maintenance by 
inert gas injection, calculation confirms the- 
oretical recovery, 49 (March JPT) 

ELK Hits: Naval Petroleum Reserve No. 1, report of 
operations 1945 through 1954, Nov. JPT, 
13 

ELKINS, LINCOLN | 
nance by Inert Gas Injection in the High 
Relief Elk Basin Field, 56 (March JPT) 

EMULSIONS. See Drilling Fluids 

ENERGY: world input, 1860-2050, relation of petro- 
leum industry, (Art.) March JPT, 9 

ENGLAND, A. C., JR., et al.: Estimating Size and Shape 
of Vertical and Horizontal Fractures, 175 
(Oct. JPT) 

ERICKSON, R. A., et al Encroachment of Injected 
Fluids beyond the Normal Well Pattern, 79 
(May JPT); discussion, 84, 85 (May JPT) 

EXPLORATION FOR OIL: 
placement 


Discussion on Pressure Mainte- 


displacement logging. See Dis- 


F 


FELSENTHAL, M., and CONLEy, F. R.: discussion of 
Analysis of Reservoir Performance k,/k 
Curves and a Laboratory k,/k, Curve Meas- 
ured on a Core Sample, 240 (Dec. JPT) 
FIELD CASE HISTORIES: gas injection project, Canal field, 
California, (Art.) Aug. JPT, 25 
highly volatile oil reservoir performance, 156 (Oct. 
JPT) 
secondary recovery operations in Tijernia No. 3 
Unit reservoir, May JPT, 23 
Shannon reservoir, Cole Creek field, Wyoming, 
(Art.) July JPT, 19 
FLASH CALCULATOR on International Business Machine 
card programmed calculator (CPC), ( Art.) 
May JPT, 9 
FiacG, A. H., et al.: Radioactive Tracers in Oil Pro- 
duction Problems, 1 (Jan. JPT) 
FLOw: gas within square drainage area, unsteady state, 
calculation, 190 (Nov. JPT) 
heat conduction in laboratory sand column filled 
with a moving hydrocarbon, (Tech. Note) 
Dec. JPT, 45 
water-oil: design and operation of models, scaling 
laws, 143 (Sept. JPT) 
FLOWMETER: subsurface, for measuring low rates of 
flow, (Tech. Note) March JPT, 53 
FLUID FLOW: equation applied in direct calculation of 
problems of vertical gas flow, 43 (March 
IPT) 


PETROLEUM TRANSACTIONS, AIME 











FLUID MAPPER: mobility ratio studies, (Tech. Note) 
Nov. JPT, 49 
FORMATION, OIL-WELL: solubility of samples in treating 
acid, (Tech. Note) Nov. JPT, 57 
FORMATION EVALUATION: coordinated, program, advan- 
tage, (Art.) July JPT, 11 
earliest recorded and subsequent developments, 
(Art.) Dec. JPT, 10 
FORMATION FRACTURES: development and application of 
frac treatments in Permian Basin, 58 (Apr. 
JPT) 
effect on waterflooding in North Texas Strawn 
sand, 16 (Feb. JPT) 
fracture treatment. See Fracture Treatment 
pressure induced: study of lost circulation with test 
cell developed at Gulf Research and De- 
velopment Co., (Tech. Note) Sept. JPT, 77 
size and shape of vertical and horizontal, model 
studies, 175 (Oct. JPT) 
unsymmetrical vertical, effect on productivity, 
model studies, (Tech. Note) Apr. JPT, 47 
FORMATION FRACTURING: fluids employed; refined oil, 
desirable characteristics, (Tech. Note) Dec. 
JPT, 39 
Foster, G. L., et al.: Investigation of Various Refined 
Oils for Formation Fracturing, (Tech. Note) 
Dec. JPT, 39 
FRAC METHOD of well stimulation (see also Formation 
Fractures): development and application in 
Permian Basin, 58 (Apr. JPT) 
FRACTURE TREATMENT: high injection rates, Spraberry 
formation, effect on performance and eco- 
nomics, (Art.) May JPT, 14 
FRACTURING FLUIDS: refined oils; desirable character- 
istics, (Tech. Note) Dec. JPT, 39 
FRAZIER, B. G., and Biron, J. N.: Field Welding on Oil- 
field Tubular Goods, ( Art.) Sept. JPT, 29 


G 


GARTHWAITE, D. L., et al.: Pressure Maintenance by 
Inert Gas Injection in the High Relief Elk 
Basin Field, 49 (March JPT): discussion, 
57 (March JPT) 

GAS DRIVE (see also Gas Lift): five-spot pattern; pre- 
dicting oil recovery in uniform sands, 7 
(Jan. JPT) 

Gas FLOW through porous media: unsteady state within 
square drainage area, calculation, 190 
(Nov. JPT) 

Gas INJECTION: Canal field, California, field case his- 
tory, (Art.) Aug. JPT, 25 

pattern injection: oil recovery, model tests, 7 (Jan. 
JPT) 

Gas LIFT: (see also Gas Drive): solution gas-drive res- 
ervoirs, free piston operations, theoretical 
and practical aspects, (Art.) Sept. JPT, 33 

GAS-OIL RATIO: solution: effect on primary recovery 
from depletion type reservoir, 121 (Aug. 
JPT) 

Gas WELLS: back-pressure method of testing; isochronal 
determination of characteristic or true slope 
of curves, 137 (Sept. JPT) 

stabilized performance curves, calculation from 
back-pressure data, (Tech. Note) May 
JPT, 53 

temperature logs using special high-speed thermom- 
eter and locating gas producing zones, 111 
(July JPT) 

Gates, C. F., and BRown, G. G.: Production Engineer- 
ing Aspects of Shannon Reservoir, Cole 
Creek Field, Wyoming, ( Art.) July JPT, 19 

GauGeEs. See Subsurface 

GEFFEN, T. M., et al.: Oil Recovery Performance of 
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Pattern Gas or Water Injection Operations 
from Model Tests, 7 (Jan. JPT); discus- 
sion, 15 (Jan. JPT) 

Goocn, F. W., JR., and ADAMs, R. H.: Report of Op- 
erations, 1945 through 1954, Naval Petro- 
leum Reserve No. 1, Elk Hills, Nov. IPT, 
13 

GOULDY, ROLAND: Discussion on Effects of Extensive 
Well Fracturing on Waterflooding in North 
Texas Strawn Sand, 21 (Feb. JPT) 

GRAVITY OF OIL: effect on primary recovery from de- 
pletion type reservoir, 121 (Aug. JPT) 

Guyop, HUBERT: analog model for study of displace- 
ment logging, 234 (Dec. JPT) 


H 


HALLER, R. W.: Big Muddy Water Flood, field case his- 
tory, Nov. JPT, 9 

HANEY, W. F., JR., et al: Estimating Size and Shape 
of Vertical and Horizontal Fractures, 175 
(Oct. JPT) 

HAWKINS, M. F., Jr.: Material Balances in Expansion 
Type Reservoirs above Bubble Point, (Tech. 
Note) Oct. JPT, 49 

HARMSEN, G. J.: Concept “Hydraulic Radius” in Por- 
ous Media, (Tech. Note) Nov. JPT, 53 

HAWKINS, M. F., Jr., and AINSWorRTH, W. J.: Two Er- 
rors in Pressure Measurement Using Sub 
surface Gauges, (Tech. Note) Sept. JPT, 80 

HAWKINS, M. F., JR., and ROEMERSHAUSER, A. E.: Ef- 
fect of Slant Hole, Drainhole, and Lateral 
Hole Drilling on Well Productivity, ( Art.) 
Feb. JPT, 11 

HENDRICKSON, A. R., et al.: Jnvestigation of Various 
Refined Oils for Formation Fracturing, 
(Tech. Note) Dec. JPT, 39 

Huitt, J. L.: Device for Determining Concentration of 
Formation Treating Acids, (Tech. Note) 
Nov. JPT, 57 

Hurst, R. E., et al.: Development and Application of 
“Frac” Treatments in the Permian Basin 
58 (Apr. JPT) 

HUTCHINSON, T. S.: Discussion on Analysis of Reser- 
voir Performance k,/k, Curves and Lab 
oratory Curve Measured on a Core Sample, 
131 (Aug. JPT) 

HYDRAULIC RADIUS: concept of use in relation to porous 
media, (Tech. Note) Nov. JPT, 53 


I 


IMAGES: method for analyzing pressure interference 
near a fault, 227 (Dec. JPT) 
INTERNATIONAL BUSINESS MACHINES CORPORATION 
See Calculator 

INTERNATIONAL PETROLEUM INDUSTRY, economics, 
(Art.) March JPT, 9 

IsLAND: filled for drilling, construction of first on Pacific 
Coast, (Art.) Dec. JPT, 13 


J 


JENKINS, R., and RAMEY, H. J., Jr.: Discussion on 
Analog Computer for Studying Heat Trans- 
fer During Thermal Recovery Process, 2\\ 
(Dec. JPT) 

JESSEN, F. W., and MILLER, J. C.: Investigation of Use 
of Spectrograph for Correlation in Lime- 
stone Rock, (Tech. Note) Feb. JPT, 41 


K 


KANSAS: salt-mud district: casing corrosion, relation be- 
tween SP and casing currents, 66 (Apr 
JPT) 

Katz, D. L., and KOBAYASHI, RIKI: Metastable Equi- 

librium in the Dew Point Determination of 







Natural Gases in the Hydrate Region, 
(Tech. Note) Aug. JPT, 51 

KERN, C. P., and CoLLiIncwoop, D. M.: Discussion on 
Encroachmeit of Injection Fluids beyond 
the Normal Well Pattern, 83 (May JPT) 

KOBAYASHI, RIKI, and Katz, D. L.: Metastable Equi- 
librium in the Dew Point Determination of 
Natural Gases in the Hydrate Region, 
(Tech. Note) Aug. JPT, 51 

KOHL, J., et al.: Locating Casing Shoe Leaks With 
Radioactive Argon, 213 (Dec. JPT) 

KREBILL, F. K., et al.: Pressure Maintenance by Inert 
Gas Injection in the High Relief Elk Basin 
Field, 49 (March JPT); discussion, 57 
(March JPT) 

KOZENY FORMULA: validity for consolidated sands, 
(Tech. Note) Nov. JPT, 53 

KRUEGER, R. F., and VoGEL, L. C.: Analog Computer 
for Studying Heat Transfer During Thermal 
Recovery Process, 205 (Dec. JPT); discus- 
sion, 211 (Dec. JPT) 

KUNZ, K. S., and Trxrer, M. P.: Temperature Surveys 
in Gas Producing Wells, 111 (July JPT) 


L 


LANDERS, J. E., et al.: External Casing Corrosion Con- 
trol, 92 (June JPT) 

LANDRUM, B. L., and CRAwForp, P. B.: Effect of 
Drainhole Drilling on Production Capacity, 
(Tech Note) Feb. JPT, 45 

Effect of Unsymmetrical Vertical Fractures on Pro- 

duction Capacity, (Tech. Note) Apr. JPT, 
47 

LANDRUM, B. L., et al.: Estimating Size and Shape of 
Vertical and Horizontal Fractures, 175 
(Oct. JPT) 

LATIMER, J. R., JR., and SHEARIN, H. M., JR.: Advan- 
tage of a Coordinated Formation Evaluation 
Program, (Art.) July JPT, 11 

LEAKS: underground gas storage reservoirs; locating by 
use of radioactive argon tracer, 213 (Dec. 
JPT) 

LeEBEAUX, J. M., and Supputn, L. F.: Theoretical and 
Practical Aspects of Free Piston Operation, 
(Art.) Sept. JPT, 33 

Lerkovits, H. C., and MATTHEWS, C. S.: Studies on 
Pressure Distribution in Bounded Reservoirs 
at Steady State, 182 (Nov. JPT) 

Lester, G. W., and Nowak, T. J.: Analysis of Pressure 
Fall-off Curves Obtained in Water Injection 
Wells to Determine Injective Capacity and 
Formation Damage, 96 (June JPT) 

LETTER SYMBOLS. See Symbols 

Lire, W. E., et al.: Fracturing at High Injection Rates, 
Spraberry Trend Area Field, (Art.) May 
JPT, 14 

LIMESTONE: geologic correlation of cores by spectro- 
graphic analysis, (Tech. Note) Feb. JPT, 41 

LINDNER, J. D., and Roark, G. E.: Effects of Extensive 
Well Fracturing on a Waterflood Operation 
in North Texas Strawn Sand, 16 (Feb. 
JPT); discussion, 21 (Feb. JPT) 

LIGNOSULFONATES as emulsifying agents: character- 
istics of emulsions and application to oil 
drilling, (Art.) June JPT, 9. 

chemistry and surface active properties, (Art.) 

June JPT, 9 

LOST CIRCULATION: pressure induced: study of fractur- 
ing with test cell developed at Gulf Re- 
search and Development Co., (Tech. Note) 
Sept. JPT, 77 


M 


Marpock, E. S., et al.: Radioactive Tracers in Oil Pro- 


214 


duction Problems, 1 (Jan. JPT) 
MarTIN, J. L., and BLUM, H. A.: Log Interpretation 
Problem in Low Resistivity Sands, ( Art.) 
Aug. JPT, 10 
MartTIN, J. L., and CAMPBELL, W. M.: Displacement 
Logginge—a New Exploratory Tool, 233 
(Dec. JPT) 
MATERIAL BALANCE EQUATION: applied to expansion 
type reservoirs above bubble point, (Tech. 
Note) Oct. JPT, 49 
simplified Schilthuis form, (Tech. Note) Jan. 
JPY, 33 
MATTHEws, C. S., and Lerxovits, H. C.: Studies on 
Pressure Distribution in Bounded Reser- 
voirs at Steady State, 182 (Nov. JPT) 
MATTHEWS, C. S., et al.: Single-fluid Five-spot Floods 
in Dipping Reservoirs, 160 (Oct. JPT) 
McDowELL, J. M.: Performance of Water Drive Reser- 
voirs, Including Pressure Maintenance, as 
Determined by the Reservoir Analyzer, 73 
(May JPT) 
McGHEE, JOHN, et al.: Wettability versus Displacement 
in Waterflooding in Unconsolidated Sand 
Columns, 227 (Dec. JPT) 
MENzIE, D. E., and CHEEK, R. E.: Fluid Mapper Model 
Studies of Mobility Ratio, (Tech. Note) 
Nov. JPT, 49 
METHANE-WATER SYSTEMS: pressure-temperature dia- 
gram, (Tech. Note) Aug. JPT, 51 
Meyer, H. I., and OrGAnNickK, E. I.: Flash Calculations 
as Performed on the IBM Card Pro- 
grammed Calculator, (Art.) May JPT, 9 
MILLER, F. G., and SEBAN, R. A.: Conduction of Heat 
Incident to Flow of Vaporizing Fluids in 
Porous Media, (Tech. Note) Dec. JPT, 45 
MILLER, J. C., and JESSEN, F. W.: Investigations of 
Use of Spectrograph for Correlation in 
Limestone Rock, (Tech. Note) Feb. JPT, 
41 
MOBILITY RATIO: effect on areal sweepout efficiency, 
fluid mapper studies, (Tech. Note) Nov. 
JPT, 49 
MODELS: electrical: displacement. See Displacement. 
effect of slant hole, drainhole, and lateral hole 
drilling on well productivity, (Art.) Feb. 
JPT) 
electrolytic tank of reservoir, study of well 
productivity increase from drainholes, 30 
(Feb. JPT) 
studies of effect of drainhole drilling on pro- 
duction capacity, (Tech. Note) Feb. JPT, 
45 
study of effect of unsymmetrical vertical frac- 
tures on production capacity, (Tech. Note} 
Apr. JPT, 4 
study of size and shape of vertical and hori- 
zontal fractures, 175 (Oct. JPT) 
fluid mapper: mobility ratio studies, (Tech. Note) 
Nov. JPT, 49 
for study of water-oil flow: scaling laws for use 
in design and operation, 143 (Sept. JPT) 


Pyrex cores: study of waterflood efficiency in oil- 
wet systems, 22 (Feb. JPT) 


sandstone: tests of oil recovery from pattern gas 
or water injections, (Jan. JPT) 
tank: reservoir flow: study of pressure distribution 
in bounded reservoirs at steady state, 182 
(Nov. JPT) 
Moore, J. M., et al.: Development and Application of 
“Frac” Treatments in the Permian Basin, 
58 (Apr. JPT) 
Moore, P. L.: Discussion on Effect of Mud Column 
Pressure on Drilling Rates, 203 (Nov. JPT) 


PETROLEUM TRANSACTIONS, AIME 










MorsB, R. A., et al.: Oil Recovery Performance of Pat- 
tern Gas or Water Injection Operations 
from Model Tests, 7 (Jan. JPT); discus- 
sion, 15 (Jan. JPT) 

MortTaDA, MOHAMED: Practical Method for Treating 
Oilfield Interference in Water-drive Reser- 
voirs, 217 (Dec. JPT) 

MuELLER, T. D., et al.: Analysis of Reservoir Perform- 
ance k,/k, Curves and a Laboratory k,/k, 
Curve Measured on a Core Sample, 128 
(Aug. JPT; discussion, 131 (Aug. JPT); 
240 (Dec. JPT) 

Murray, A. S., and CUNNINGHAM, R. A.: Effect of 
Mud Column Pressure on Drilling Rates, 
196 (Nov. JPT); discussion, 203 (Nov. 
JPT) 

Myers, J. L. P., et al.: Radioactive Tracers in Oil Pro- 
duction Problems, 1 (Jan. JPT) 


N 


NATHAN, C. C.: Solubility Studies on High Molecular 
Weight Paraffin Hydrocarbons Obtained 
from Petroleum Rod Waxes, 151 (Sept. 
JPT) 

NATURAL Gas: direct calculation of bottom-hole pres- 
sures in wells, 43 (March JPT) 

in equilibrium with gas hydrates: dew point and 
equilibrium temperatures, differences, 
(Tech. Note) Aug. JPT, 5! 

NATURAL GAS WELLS. See Gas Wells 

NAVAL PETROLEUM RESERVE No. 1, Elk Hills, report of 
operations 1945 through 1954, Nov. JPT, 
13 

NEUTRON LOGs. See Porosity 

NEWACHECK, R. L., et al.: Locating Casing Shoe Leaks 
With Radioactive Argon, 213 (Dec. JPT) 

NEWCOMBE, JACK, et al.: Wettability versus Displace- 
ment Waterflooding in Unconsolidated Sand 
Columns, 227 (Dec. JPT) 

Nowak, T. J., and Lester, G. W.: Analysis of Pres- 
sure Fall-off Curves Obtained in Water In- 
jection Wells to Determine Injective Ca- 
pacity and Formation Damage, 96 (June 


IPT) 
O 


OFFSHORE DRILLING: See Drilling 
OIL PRODUCTION: summary of technological d elop- 
ments from first well to present time ( Art.) 
Dec. JPT, 10 
OIL RECOVERY: depletion equation, derivat 
lution, 121 (Aug. JPT) 
from depletion type reservoirs: c relation with 
relative permeability ratio, \PI gravity, and 
solution gas-oil ratio, 120 Aug. JPT) 
thermal: electrical analog for solv: \g unsteady-state 
heat conduction problems of radial system 
containing a moving, cylindrical heat source, 
205 (Dec. JPT) 
heat conduction in laboratory sand column 
filled with a moving hydrocarbon, (Tech. 
Note) Dec. JPT, 45 
OLANSEN, N. A.: Secondary Recovery Operations in the 
Tijerina No. 3 Unit Reservoir, May JPT, 23 
O’QuINN, P. C., et al.: Fracturing at High Injection 
Rates, Spraberry Trerid Area Field, (Art.) 
May JPT, 14 
ORGANICK, E. I., and Meyer, H. I.: Flash Calculations 
as Performed on the IBM Card Pro- 
grammed Calculator, (Art.) May JPT, 9 
Ososa, J. S., et al.: Differences in Behavior of Fresh 
and Aged East Texas Woodbine Cores, 86 
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(June JPT); discussion, 91 (June JPT) 


P 


PARAFFIN: high meiecular weight: solubilities in hydro- 
carbon solvents and in chlorinated and 
oxygenated solvents, 151 (Sept. JPT) 

PARAFFIN DEPOSITION: mechanism and prevalence, 
(Art.) Sept. 23 

solvent response, (Art.) Sept. JPT, 23 
viscosity-temperature curves, value, (Art.) Sept. 
JPT, 23 

PEACEMAN, D. W., et al.: Calculation of Unsteady-state 
Gas Flow within a Square Drainage Area, 
190 (Nov. JPT) 

PERMIAN BasINs frac treatments, development and ap- 
plication, 58 (Apr. JPT) 

West Texas: geologic correlation of limestone 
cores by spectrographic analysis, (Tech. 
Note) Feb. JPT, 41 

PERRINE, R. L.: Well Productivity Increase from Drain- 
holes as Measured by Model Studies, 30 
(Feb. JPT) 

PETROLEUM INDUSTRY: by-product developments and 
trends as seen by an editor, (Art.) March 
JPT, 14 

international, economics, (Art.) March JPT, 9 

PERKINS, F. M., Jr., et al: Differences in Behavior of 
Fresh and Aged East Texas Woodbine 
Cores, 86 (June JPT); discussion, 91 (June 
JPT 

Pincs, C. J., JR., and TEMPELAAR-LIETZ, W.: Canal 
Field Gas Injection Project, field case his- 
tory, (Art.) Aug. JPT, 25 

Piston, free. See Gas Lift 

POROUS MEDIA: concept “hvdraulic radius,” (Tech. 
Note) Nov. JPT ; 

Prats, M., et al.: Single}. *’. -spot Floods in Dip- 

ping R«; voirs, 160 \Uct. »PT) 

.iNTENANCE: by inert gas injection, Elk 

Basin Tensleep reservoir, 49 (March JPT) 

PR’ .SURE MEASUREMENT with subsurface gauges, pre- 
vention of two errors in sealing gauge, 
(Tech. Note) Sept. JPT, 80 

PoHoRILes, E. M., et al.: New Apparatus for Studying 
Pressure Induced Fracturing in Reference 
to Lost Circulation, (Tech. Note) Sept. 
JPT, 77 

Propuctivity: drainholes for increasing, model studies, 
30 (Feb. JPT) 

effect of slant hole, drainhole, and lateral hole 
drilling, (Art.) Feb. JPT, 11 

effect of unsymmetrical vertical formation frac- 
tures, model studies, (Tech. Note) Apr. 
JPT, 47 

Porosity, thin beds: determination from neutron logs: 
concept of drag as related to logging speed 
and time constant, effects of neutron source 
to detector spacing, and their relationships 
to deflection of neutron curve, 132 (Aug. 
JPT) 

PyLes, E. E.: First Pacific Coast Permanent Drilling 
Island Modified to Allow Drilling of 70 
Wells, (Art.) Dec. JPT, 13 


R 


RACHFORD, H. H., JR., and SCHULTZ, W. P.: Some Use- 
ful Tables for Approximating Smooth 
Curves by Fifth and Lower Degree Poly- 
nomials, (Tech. Note) Dec. JPT, 43 

RACHFORD, H. H., Jr., et al.: Calculation of Unsteady- 
state Gas Flow within a Square Drainage 
Area, 190 (Nov. JPT) 

Rapp, F. J., and DeWitte, L.: Corrosion of Oil Well 


PRESSURP 
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Casing by Earth Currents, 66 (Apr. JPT) 

RADIOACTIVE ARGON: tracer in locating casing shoe 
leaks, 213 (Dec. JPT) 

RADIOACTIVE TRACERS: locating zones of permeability, 
1 (Jan. JPT) 

problems and successes, | (Jan. JPT) 
safe handling, | (Jan. JPT) 

RaMeyY, H. J., JR., and JENKINS, R.: Discussion on 
Analog Computer for Studying Heat Trans- 
fer During Thermal Recovery Process, 211 
(Dec. JPT) 

RaMsEy, D. E., et al.: Development and Application 
of “Frac” Treatments in the Permian Basin, 
58 (Apr. JPT) 

Rapoport, L. A.: Scaling Laws for Use in Design and 
Operation of Water-oil Flow Models, 143 
(Sept. JPT) 

RELATIVE PERMEABILITY: effect on primary recovery 
from depletion type reservoir, 121 (Aug. 
JPT) 

performance type and laboratory k,/k, curves for 
hypothetical reservoir, differences, 128 
(Aug. JPT) 

RESERVOIR ANALYZER: determination of performance of 
water-drive reservoirs, including pressure 
maintenance, 73 (May JPT) 

RESERVOIR ENGINEERING: letter symbols, report of com- 
mittee, (Art.) Jan. JPT, 37 

RESERVOIR PERFORMANCE: k,/k, curves for hypothetical 
reservoir contrasted with laboratory curve 
measured on core sample, 128 (Aug. JPT) 

RESERVOIR PRESSURE: study for bounded reservoirs at 
steady state, 182 (Nov. JPT) 

RESERVOIR ROCK: coordinated formation evaluation 
program, advantage, (Art.) July JPT, 11 

RESERVOIRS: bounded: drainage area, estimating shape 
for each well, 182 (Nov. JPT) 

case history of performance of highly volatile type, 
156 (Oct. JPT) 

expansion type, newly discovered, estimation of 
size by material balance equation, (Tech. 
Note) Oct. JPT, 49 

RICHARDSON, J. G., et al.: Differences in Behavior of 
Fresh and Aged East Texas Woodbine 
Cores, 86 (June JPT); discussion, 91 (June 
JPT) 

RosBerts, T. G., and Arps, J. J.: Effect of Relative 
Permeability Ratio, Oil Gravity, and Solu- 
tion Gas-oil Ratio on Primary Recovery 
from Depletion Type Reservoir, 120 (Aug. 
JPT) 

Roark, G. E., and LinpNerR, J. D.: Effects of Extensive 
Well Fracturing on a Waterflood Operation 
in North Texas Strawn Sand, 16 (Feb. 
JPT): discussion, 21 (Feb. JPT) 

ROEMERSHAUSER, A. E., and HAWKINS, M. F., Jr.: 
Effect of Slant Hole, Drainhole, and Lateral 
Hole Drilling on Well Productivity, ( Art.) 
Feb. JPT, 11 

Ropcers, J. K., et al.: New Apparatus for Studying 
Pressure Induced Fracturing in Reference 
to Lost Circulation, (Tech. Note) Sept. 
JPT, 77 

RoseENE, R. B., et al.: /nvestigation of Various Refined 
Oils for Formation Fracturing, (Tech. Note) 
Dec. JPT, 39 

RUMBLE, R. C.: A Subsurface Flowmeter, (Tech. Note) 
March JPT, 53 

Rzasa, M. J., et al.: Wettability versus Displacement in 
Waterflooding in Unconsolidated Sand Col- 
umns, 227 (Dec. JPT) 


S 


SCHWARZ, N., and Croes, G. A.: Dimensionally Scaled 
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Experiments and the Theories on the Water- 
drive Process, 35 (March JPT) 

SCHULTZ, W. P., and RACHFoRD, H. H., Jr.: Some 
Useful Tables for Approximating Smooth 
Curves by Fifth and Lower Degree Poly- 
nomials, (Tech. Note) Dec. JPT, 43 

SCINTILLOMETER: determination of thickness of thin 
beds, 132 (Aug. JPT) 

SEBAN, R. A., and MILLER, F. G.: Conduction of Heat 
Incident to Flow of Vaporizing Fluids in 
Porous Media, (Tech. Note) Dec. JPT, 45 

SECONDARY RECOVERY: encroachment of injected fluids 
beyond normal well pattern, estimation, 79 
(May JPT) 

sweepout pattern efficiencies in area between last 
row of wells and reservoir boundary, lab- 
oratory study, 79 (May JPT) 
Tijerina No. 3 unit reservoir, Texas, case history, 

May JPT, 23 

SHANNON RESERVOIR, Cole Creek field, Wyoming, field 
case history, (Art.) July JPT, 19 

SHEARIN, H. M., JR., and LATIMER, J. R., JR.: Advan- 
tage of a Coordinated Formation Evaluation 
Program, (Art.) July JPT, 11 

SHock, D. A., et al.: External Casing Corrosion Con- 
trol, 92 (June JPT) 

SHOCK, D. A., et al.: Studies of the Mechanism of Par- 
affin Deposition and Its Control, (Art.) 
Sept. JPT, 23 

SILICA SAND: treating with silicone polymer to proauce 
surfaces of varying degrees of wettability, 
227 (Dec. JPT) 

SIMPSON, A. L., and Epwarps, J. M.: Method for Neu- 
tron Derived Porosity Determination for 
Thin Beds, 132 (Aug. JPT) 

SLoBop, R. L., et al.: Encroachment of Injected Fluids 
beyond the Normal Well Pattern, 79 (May 
JPT; discussion, 84, 85 (May JPT) 

SP LoG: water wells in Brazil, interpretation, (Tech. 
Note) Nov. JPT, 59 

SP PROFILE: quantitative derivation, 66 (Apr. JPT) 

SPECTROGRAPHIC ANALYSIS: geologic correlation of 
limestone cores from wells in Permian Basin 
of West Texas, (Tech. Note) Feb. JPT, 41 


value in studying various phases of sedimentation, 
(Tech. Note) Feb. JPT, 41 
SPRABERRY TREND AREA FIELD: fracture treatment at 
high injection rates, (Art.) May JPT, 14 


STEWART, F. M., et al.: Pressure Maintenance by Inert 
Gas Injection in the High Relief Elk Basin 
Field, 49 (March JPT): discussion, 57 
(March JPT) 

STRAWN LIME FIELD: coordinated formation evaluation 
program, (Art.) July JPT, 11 

STRICKLER, W. R.., et al.: Single-fluid Five-spot Floods 
in Dripping Reservoirs, 160 (Oct. JPT) 

SUBSURFACE GAUGES for pressure measurement: pre- 
vention of two errors in sealing gauge, 
(Tech. Note) Sept. JPT, 80 

Supsury, J. D., et al.: External Casing Corrosion Con- 
trol, 92, (June JPT) 

Studies of the Mechanism of Paraffin Deposition 

and Its Control, (Art.) Sept. JPT, 23 

SuppuTH, L. F., and LesBeaux, J. M.: Theoretical and 
Practical Aspects of Free Piston Operation, 
(Art.) Sept. JPT, 33 

SUKKAR, Y. K., and CoRNELL, Davin: Direct Calcula- 
tion of Bottom-hole Pressures in Natural 
Gas Wells, 43 (March JPT) 

SUMAN, J. R.: Oil Production—a Continuing Challenge, 
(Art.) Dec. JPT, 10 


PETROLEUM TRANSACTIONS, AIME 











committee, (Art.) Jan. JPT, 37 


T 


TEMPELAAR-LiETZ, W., and PinGs, C. J., Jr.: Canal 
Field Gas Injection Project, field case his- 
tory, (Art.) Aug. JPT, 25 

Tepuitz, A. J., et al.: New Apparatus for Studying 
Pressure Induced Fracturing in Reference 
to Lost Circulation, (Tech. Note) Sept. 
JPT, 77 

TERRY, J. M., et al: Radioactive Tracers in Oil Produc- 
tion Problems, 1 (Jan. JPT) 

Texas: North: Strawn sand, effect of formation frac- 
turing on waterflooding, 16 (Feb. JPT) 

Tijerina. See Tijerina. 
West: Permian Basin. See Permian Basin 
well completions for waterflooding, Permian 

sand, (Art.) Jan. JPT, 21 

TIJERINA No. 3 RESERVOIR, TEXAS: secondary recovery 
operations, case history, May JPT, 23 

TixieR, M. P., and KUNZ, K. S.: Temperature Surveys 
in Gas Producing Wells, 111 (July JPT) 

Torrey, P. D.: Preparation of Water for Injection into 
Oil Reservoirs, (Art.) Apr. JPT, 9 

TRACERS: radioactive. See Radioactive 

Tracy, G. W.: Simplified Form of the Material Balance 
Equation, (Tech. Note) Jan. JPT, 53 

Twalits, W. O.: Canadian Oil at the Eighth Milestone, 
(Art.) July JPT, 7 


V 


VARJAO, PAULO DE ANDRADE: Results of SP Log In- 
terpretation, ‘Tech. Note) Nov. JPT, 59 

VoGEL, L. C., and KRUEGER, R. F.: Analog Computer 
for Studying Heat Transfer During Thermal 
Recovery Process, 205 (Dec. JPT); discus- 
sion, 211 (Dec. JPT) 


Ww 


WALNUT SHELLS: beneficial addition to cement slurries, 
(Art.) Aug. JPT, 15 

Warp, T. G., JR., et al.: Well Completions for Water- 
flood Operations, (Art.) Jan. JPT, 21 

WARREN, JOSEPH E., and CALHOoUN, J. C., Jr.: Study of 
Waterflood Efficiency in Oil-wet Systems, 
22 (Feb. JPT) 

WarREN, J. E., et al.: Analysis of Reservoir Perform- 
ance k,/k, Curves and a Laboratory k,/k.,, 
Curve Measured on a Core Sample, 128 
(Aug. JPT); discussion, 131 (Aug. JPT); 
249 (Dec. JPT) 

WATER DRIVE: efficiency: model displacement experi- 
ments, continuation of Engelberts and 
Klinkenberg series, 35 (March JPT) 

WATER-DRIVE RESERVOIRS: change in pressure decline 
and water influx as function of rate of fluid 
withdrawals, as function of permeability of 
formation, and as function of size of aqui- 
fer, 73 (May JPT) 

common aquifer: pressure interference; analog and 
digital computations, 217 (Dec. JPT) 
analysis by charts based on principle of super- 
position, 217 (Dec. JPT) 
near a fault, use of method of images, 217 
(Dec. JPT) 
Woodbine formation, 217 (Dec. JPT) 
position of water injection wells with reference to 
oil-water contact, 73 (May JPT) 
pressure maintenance, effect on rate of natural 
water influx, 73 (May JPT) 

WATERFLOODING: Big Muddy field, Wyoming, field case 

history, Nov. JPT, 9 

1955 


VOL, 204, 


SYMBOLs: letter, for reservoir engineering, report of 





economics, Expected Gross Revenue less Summa- 
tion of Computed Acquisition Costs, Costs 
of Development, and Operating Costs (over 
life of flood) equals Estimated Gross Profit, 
(Art.) Jan. JPT, 12 
effect of formation fracturing in North Texas 
Strawn sand, 16 (Feb. JPT) 
efficiency in oil-wet systems: study with Pyrex 
cores, 22 (Feb. JPT) 
essentials for success, engineering, experience, eco- 
nomics, (Art.) Jan. JPT, 12 
five-spot pattern: predicting oil recovery in uniform 
sands, 7 (Jan. JPT) 
single fluid in dipping reservoirs, 160 (Oct. 
JPT) 
mobility ratio effect on areal sweepout efficiency, 
fluid mapper studies, (Tech. Note) Nov. 
JPT, 49 
pattern water injection: oil recovery, model tests, 
7 (Jan. JPT) 
preparation of water: design and operation of 
plants, (Art.) Apr. JPT, 17 
methods, objects, and results, (Art.) Apr. 
JPT, 9 
radioactive tracers, problems and successes, | (Jan. 
JPT) 
theories evaluated by dimensionally scaled experi- 
ments, 35 (March JPT) 
unconsolidated sand: wettability versus displace- 
ment, laboratory study, 227 (Dec. JPT) 
water: calculation of requirements, selection and 
development of source, (Art.) Jan. JPT, 18 
well completion for injection and producing wells, 
Permian sand, West Texas, (Art.) Jan. 
JPT, 21 
WATER INJECTION WELLS: pressure fall-off curves: gen- 
eral characteristics, 96 (June JPT) 
theoretical quantitative analysis and application of 
theory to field data, 96 (June JPT) 
WATER WELLS: Brazil, SP log interpretation, (Tech. 
Note) Nov. JPT, 59 
WATKINS, J. W.: Design and Operation of Plants for 
Preparation of Water for Injection into Oil 
Reservoirs, (Art.) Apr. JPT, 17 
Watson, M. S., JR., et al.: Well Completions for Water- 
flood Operations, (Art.) Jan. JPT, 21 
WELDING: at rig site: five categories discussed, ( Art.) 
Sept. JPT, 29 
field: technique all important on steels used for 
oilfield tubular goods, (Art.) Sept. JPT, 29 
WELL COMPLETION: for injection and producing wells, 
Permian sand, West Texas, (Art.) Jan. 
JPT, 21 
radioactive tracers for cement top location, 5 (Jan. 
JPT) 
West, W. J., et al: Analysis of Reservoir Performance 
k,/k, Curves and a Laboratory k,/k, Curve 
Measured on a Core Sample, 128( Aug. 
JPT); discussion, 131 (Aug. JPT); 240 
(Dec. JPT) 
WILLSON, C. O.: Some Sidelights in Petroleum’s Eco- 
nomic Progress, (Art.) March JPT, 14 
WoLF, ALEXANDER: Discussion on Encroachment of 
Injection Fluids beyond the Normal Well 
Pattern, 84 (May JPT) 
WoopDBINE FORMATION, East Texas: pressure interfer- 
ence, 217 (Dec. JPT) 
Woops, R. W.: Case History of Reservoir Performance 
of Highly Volatile Type Oil Reservoir, 156 
(Oct. JPT) 
WYOMING: Big Muddy field, waterflooding, field case 
history, Nov. JPT, 9 
Shannon reservoir. See Shannon wae 












AIME, PETROLEUM BRANCH 
MAJOR MEETING SCHEDULE 


1956—Feb. 19-23, Annual AIME Meeting, Hotel 
Statler, New York City 


1956—March 29-30, “Petroleum Conference, Eco- 
nomics and Valuation,” Statler-Hilton Hotel, 
Dallas. 

1956—May 14-15, ‘‘Formation Fracturing Sympo- 
sium,” Wichita, Kans. 

1956—May 17-18, Joint Meeting of Rocky Mountain 
Petroleum Sections of AIME, Gladstone Hotel, 
Casper, Wyo. 

1956—Nov. 19-20, “North Texas Secondary Recov- 
ery Symposium,” Wichita Falls, Tex. 

1956—Oct. 14-17, Petroleum Branch Fall Meeting, 

Biltmore Hotel, Los Angeles. 
1957—Feb. 24-28, Annual AIME Meeting, New 
leans. 

1957—Oct. 6-9, Petroleum Branch Fall Meeting, 
Dallas. 

1958—Feb. 16-20, Annual AIME Meeting, Hotel 
Statler, New York City. 

1958—Oct. 12-15, Petroleum Branch Fall Meeting, 
Houston. 

1959—Feb. 15-19, Annual AIME Meeting, San 
Francisco. 

1959—Oct. 11-14, Petroleum Branch Fall Meeting, 
Dallas (tentative). 

1960—Feb. 14-18, Annual AIME Meeting, Hotel 
Statler, New York City. 

1961—Feb. 19-23, Annual AIME Meeting, St. Louis. 

1962—Feb., Annual AIME Meeting, New York City. 

1963—Feb. 17-22, Annual AIME Meeting, Dallas. 
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® continuous Drilling Mud Research 


In Magcobar’s modern new headquarters at Houston, three research 
laboratories and dozens of chemists and engineers are constantly searching 
for ways to make your drilling faster, more economical, and more efficient. They are 
concerned with the problems of even deeper wells, heavier muds, higher 
temperatures and the variety of other complex drilling 
problems that will give producers trouble tomorrow, 
if not solved today. The tangible contributions to 
drilling mud progress by Magcobar in the past give concrete 
evidence of the continued progress you can expect 
from Magcobar in the future. Continuous research is 
just another one of the many services that Magcobar gladly 
accepts as its individual responsibility to the oil industry. 
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